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Apresentagao 


Atualmente, as empresas estao muito preocupadas com a redugao de custos. Tornar 
suas operagoes mais eficientes e um caminho potencialmente capaz de atingir esse 
objetivo. Para que a eficiencia operacional seja alcangada, no entanto, muitas decisoes 
precisam ser tomadas pelo nfvel gerencial da empresa. Algumas dessas decisoes sao 
criticas e suas consequencias podem ser formidaveis ou desastrosas no que diz res- 
peito a eficiencia da operagao. 

A busca pelas decisoes otimas nem sempre e um processo simples, mas muitas 
vezes e encarado como tal, carecendo, em boa parte do emprego de um cientificismo 
maior. O processo de tomada de decisao e, frequentemente, apoiado apenas pela 
experiencia profissional e intuigao do decisor. Desta forma, e de suma importancia 
pesquisar, analisar, criticar, desenvolver e aplicar ao cotidiano operacional das em¬ 
presas metodos quantitativos de previsao e de apoio a busca pelas decisoes otimas, 
que objetivam complementar (e nao substituir!) a experiencia profissional e a intuigao 
do tomador de decisao. 

Em outras palavras, e vital explorar ferramentas capazes de embasar as decisoes 
operacionais a serem tomadas pelo administrador, assim como justifica-las e comunica- 
las para toda a empresa. Nesse sentido, o avango das tecnologias computacionais 
pode facilitar ainda mais a aplicagao de metodos quantitativos de previsao e de apoio 
a tomada de decisao. 

Pois a proposta deste livro e justamente explorar essas ferramentas. Para tal, a obra e 
dividida em duas partes: metodos de apoio a tomada de decisao e metodos de previsao. 
Em ambas, cada capftulo apresenta um metodo diferente, com uma revisao conceitual 
e o passo a passo de sua aplicagao a algum caso pratico, alem da apresentagao e dis- 
cussao de resultados, com as respectivas planilhas de implementagao. 

Na Parte I, a Programagao Linear, as Arvores de Decisao e a Simulagao de Mon¬ 
te Carlo sao apresentadas teoricamente e aplicadas a problemas reais na area de 
Administragao. O ultimo capftulo dessa parte realiza uma analise comparativa destas 
metodologias em relagao a outro problema real. 

A Parte II consiste em dois capftulos que aplicam a Regressao Multipla (com 
variaveis dummy) a problemas reais de previsao na area de negocios, e um tutorial (com 
texto e planilha) sobre um complexo metodo de series temporais, a Decomposigao 
Classica. 

Este trabalho e destinado a alunos, professores e pesquisadores dos cursos de 
graduagao, pos-graduagao e politecnico em Administragao, Engenharia de Produgao e 
Logfstica, inserindo-se perfeitamente como livro-texto para disciplinas como Metodos 
Quantitativos, Metodos Quantitativos Aplicados e Metodos de Apoio a Decisao. Pode 
ser tambem utilizado como bibliografia complementar para disciplinas como Previsao 
de Vendas, Analise de Decisoes, Pesquisa Operacional e Otimizagao em Sistemas 
Logfsticos. 
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Os profissionais de mercado tambem podem usufruir dessa obra para consultas 
visando aplicagoes proprias. Analistas de Loglstica, Inteligencia de Marketing e de 
Planejamento podem encontrar grande utilidade na leitura desta obra para a realizagao 
de suas atividades de transportes, previsao de vendas e dimensionamento de equipes. 

Desta forma, este trabalho procura preencher uma lacuna existente na interface 
entre a academia e o mundo dos negocios, apresentando, alem do arcabougo teorico, 
pesquisas academicas (e seus resultados!) recentes na area, totalmente aplicadas e em 
linha com as metodologias abordadas. 


Marco Aurelio Caririo Bouzada 
Organizador 
Fevereiro de 2013 




Parte I 


Metodos de apoio 
a tomada de decisao 


Os metodos de apoio a tomada de decisao podem ajudar as empresas a, por 
exemplo: 

• decidir que tipo de tecnologia utilizar para a obtengao de um novo produto, 
considerando os riscos envolvidos e as consequencias potenciais de cada 
alternativa; 

• designar elementos de uma equipe para a realizagao de diversas tarefas de 
forma otimizada; 

• optar pela modalidade de seguro mais interessante em face dos riscos ineren- 
tes ao negocio e aos custos e das recompensas referentes a cada alternativa 
considerada; 

• alocar recursos em diferentes opgoes de rmdia de forma eficiente; 

• dimensionar corretamente a produgao de determinado item, tendo em mente 
a flutuagao natural da demanda e eventuais custos de sobra e de falta in- 
corridos; 

• decidir onde instalar uma antena de transmissao, levando em conta o que 
pode acontecer se cada um dos possiveis locais for escolhido. 

Obviamente, tais decisoes podem ser tomadas sem o auxilio de um metodo de 
apoio e, lamentavelmente, o sao com alguma frequencia. Quando isso acontece, 
no entanto, muitas vezes as consequencias de tais decisoes tomadas de forma 
intuitiva ou baseada apenas na experiencia profissional do gestor podem ser 
desastrosas: 

• a tecnologia escolhida pode nao ser a mais indicada para o caso, resultando 
em custos associados muito maiores que os ganhos; 

• funcionarios podem ser alocados a tarefas de forma que elas sejam realizadas 
em um periodo muito maior do que o necessario; 

• o seguro escolhido, por ser o mais barato, pode nao cobrir satisfatoriamente 
os riscos da operagao de uma empresa; 





2 ) PARTE | 1 Metodos de apoio a tomada de decisao 


ELSEVIER 


• as mfdias escolhidas para a alocagao de recursos podem alcangar uma 
quantidade insatisfatoria de potenciais consumidores; 

• o mal dimensionamento da produgao pode fazer com que um item fique 
"encalhado" no estoque; 

• o local escolhido para a instalagao de uma antena pode ser, tecnicamente, 
o mais interessante, mas proibitivamente caro. 

Justamente por isso, acaba se tornando um gestor diferenciado aquele que 
procura minimizar a ocorrencia de consequencias como as descritas anterior- 
mente, com a utilizagao de metodos quantitativos para subsidiar a tomada de 
decisao empresarial. Tais metodos nao tern como objetivo substituir o feeling 
e o bom senso do tomador de decisao; apenas complementa-los, dando apoio 
ao processo. 

Os metodos quantitativos de apoio a tomada de decisao (que sao o objeto 
desta Parte I) podem ser divididos em dois tipos: os determinfsticos e os nao 
determinfsticos (estes ultimos podendo ser subdivididos em probabilfsticos e 
nao probabilisticos). 

Os metodos determinfsticos so podem ser usados satisfatoriamente nas chamadas 
decisSes em condigoes de certeza (quando nao ha duvida a respeito dos valores 
das variaveis relevantes). O representante classico e a Programagao Linear, cujo 
caso particular (Programagao Inteira) e explorado no Capftulo 1, que apresenta 
uma tentativa de solucionar um problema caracterfstico de logfstica da Marinha 
do Brasil: a alocagao, o transporte e a distribuigao de generos frigorificados para 
as Organizagoes Militares. 

Ja os metodos nao determinfsticos nao probabilfsticos sao adequados para as de¬ 
cisSes em condigoes de incerteza, que ocorrem quando nao se tern praticamente 
nenhuma informagao a respeito de algumas variaveis relevantes, a nao ser os 
seus possfveis valores de ocorrencia (mas sem as respectivas probabilidades). As 
alternativas mais utilizadas em tais situagoes sao os criterios Maximax, Maximin 
e Minimax. 

Por fim, os metodos probabilfsticos podem ser aplicados em situagoes consi- 
deradas no meio do caminho entre as condigoes de certeza e as condigoes de 
incerteza: as decisSes em condigoes de risco - quando nao se sabe ao certo os 
valores de algumas variaveis relevantes, mas sao conhecidas as suas distribuigoes 
de probabilidades. Em tais situagSes, dois dos metodos mais utilizados sao a 
Arvore de Decisao e a Simulagao. 

Quando ha diferentes momentos em que e necessaria uma tomada de decisao, 
com ocorrencia de eventos probabilfsticos entre eles, a Arvore de Decisao e 
mais recomendada. O Capftulo 2 usa esta metodologia para avaliar a escolha 
de um fornecedor sob condigoes de risco. 
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Se, no entanto, a situagao compreende muitas complexidades a serem modeladas 
e/ou nao se esta preocupado apenas com a decisao otima, mas tambem com o 
desempenho do sistema como um todo em face de cada uma das altemativas 
de decisao, entao a Simulagao e mais adequada. 

Tal ferramenta e utilizada no Capitulo 3 para verificar se a intermodalidade de 
transporte e uma alternativa viavel no processo de escoamento de arroz. 

O Capitulo 4, o ultimo da Parte I, utiliza como pano de fundo um problema real 
de determinagao do mix de projetos de consultoria para comparar a adequagao 
de todas estas metodologias citadas: Programagao Linear (de fato, o seu caso 
particular, a Programagao Binaria), criterios Maximax, Maximin e Minimax, 
Arvore de Decisao e Simulagao. 






Capftulo 1 


A aplica^ao da Programagao Inteira 
na solu^ao logfstica do transporte 
de carga na Marinha do Brasil 


Ricardo Franga Santos 

MADE/UNESA 

Eugenio Correa de Souza Junior 

MADE/UNESA 

Marco Aurelio Carino Bouzada 

MADE/UNESA 


Objetivos de aprendizagem do capftulo 







Saber modelar um problema real em termos determimsticos. 

Entender o conceito de variaveis de decisao, fungao-objetivo e restrigoes 
em um modelo de Programagao Linear. 

Aprender a usar o Solver do Excel para otimizar modelos de Programagao Linear. 
Entender as peculiaridades da Programagao Inteira e da Programagao Binaria. 
Conhecer algumas das limitagbes do Solver do Excel na busca pela solugao otima. 


INTRODUgAO 

A Marinha do Brasil realiza suas atividades de abastecimento por meio de um sistema 
constituido por orgaos de distribuigao, processos e recursos de qualquer natureza, 
interligados e interdependentes, estruturado com a finalidade de prover generos neces- 
sarios a manutengao das Organizagoes Militares (OMs). Para isso conta com depositos 
primarios, depositos secundarios e depositos regionais. 

Dentre os depositos regionais existe o Deposito de Subsistencia da Marinha no 
Rio de Janeiro (DepSubMRJ), criado em 1955, cuja missao e executar as tarefas de 
receber, periciar, estocar, controlar e fornecer generos alimentfcios, a fim de contribuir 
para a prontidao operativa dos Meios Navais, Aeronavais, Fuzileiros Navais e demais 
Organizagoes Militares da Marinha. 


5 
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Em relagao a fungao logfstica de transporte, diversos requisites devem ser atendidos 
pelo DepSubMRJ, principalmente em relagao a temperatura da carga e do vefculo, as- 
sim como em relagao a outras caracterfsticas deste ultimo, como estado de conservagao, 
de higiene e de seguranga. 

O DepSubMRJ dispoe de um numero limitado de vefculos de transporte de cargas 
que atendem a esses requisitos, que nem sempre estao disponfveis devido a problemas 
mecanicos ou ausencia de motoristas. 

A frota de vefculos e constitufda de caminhoes leves e ageis que carregam cargas 
menores, e de caminhoes pesados para cargas maiores e percursos mais longos. A dis- 
tribuigao da frota e realizada de modo empfrico, baseada na experiencia e competencia 
dos profissionais, a partir das quais e produzida uma programagao semanal, sustentada 
em fatores como demanda, localizagao do destinatario, condigoes de acesso, distancia, 
disponibilidade do motorista e consumo de combustfvel, entre os principals fatores. 
Por outro lado, existe uma demanda, por parte das OMs, no que tange aos generos 
frigorificados, que necessitam ser atendidos por essa frota. A programagao semanal 
atende tanto as OMs na Cidade do Rio de Janeiro como as das cidades proximas. As 
entregas sao realizadas de segunda-feira a sexta-feira, inclusive feriados. Na Tabela 1.1 
estao listados os custos fixos e variaveis da operagao. 

Na Tabela 1.2 estao registrados os valores mfnimo e maximo de entrega diaria. 

Na Tabela 1.3 podem ser observadas a distancia e a demanda semanal de cada OM. 

O transporte representa um significativo custo logfstico por parte das empresas, 
absorvendo de um a dois tergos do total. Quando realizado por frota e equipamentos 
proprios, pode propiciar melhor desempenho operacional, maior disponibilidade e 


TABELA 1.1 Custos fixos e variaveis 


Custos fixos (em R$/mes) 


Custos variaveis (em R$/km) 


Depreciagao do vefculo - 13.500,00 


Manutengao - 2,25 


Licenciamento e seguro - 2.025,00 


Combustfvel - 1,00 


Reserva para seguro - 1.600,00 


Lubrificagao - 0,05 


Salarios (motorista e ajudante) - 22.500,00 


Pneus - 0,20 


TABELA 1.2 Entrega minima e maxima por dia (em toneladas) 


Segunda 


Terga 


Quarta 


Quinta 


Sexta 


Entrega minima 


31 


25 


37 


29 


13 


Entrega 

maxima 


45 


32 


43 


33 


38 
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TABELA 1.3 Distancia e demanda semanal de cada OM 

"N 

OM 

Demanda semanal (em t) 

Distancia (em km) 

Alfa 

8 

12 

Bravo 

15 

5 

Charlie 

16 

140 

Delta 

18 

160 

Echo 

20 

15 

Foxtrot 

26 

10 

Golf 

19 

180 

Hotel 

14 

06 

India 

16 

25 

Juliet 

13 

20 

V 


_^ 


capacidade de transporte e gerar custos menores. Isto geralmente e possfvel quando o 
volume de carga e elevado, sendo mais economico possuir um servigo de transporte 
proprio do que contrata-lo. Por outro lado, nem sempre isso e possfvel devido as 
exigencias de requisitos especiais, tais como entrega rapida com confiabilidade muito 
elevada, indisponibilidade de algum equipamento especial, manuseio especial da carga 
ou um servigo que deva estar disponfvel quando necessario. 1 

Na analise realizada neste capftulo, sera considerado apenas o servigo de transporte 
por frota propria, em que os custos fixos e variaveis sao perfeitamente conhecidos. 

O gerente de trafego de carga do DepSubMRJ e responsavel pela decisao de utiliza- 
gao da frota, em que menores custos e melhor desempenho na entrega sao perseguidos. 

A montagem do piano de rota e um problema encontrado para direcionar os 
vefculos atraves das diversas vias. O calculo da rota pode ser feito pela minima dis¬ 
tancia, pelo mfnimo tempo ou por uma combinagao de ambos. Alem disso, o custo 
total e dependente nao so do custo do transporte, mas tambem do custo de estocagem. 2 

Por outro lado, as empresas estocam produtos visando a melhoria da compensagao 
entre oferta e demanda, alem de diminuir seus custos totais. Neste sentido, os custos de 
armazenagem e manuseio de materiais podem ser justificados pela redugao dos custos 
de transporte e de produgao, pois os estoques podem reduzir os custos de transporte 
ao permitir que quantidades maiores e mais economicas de carga sejam transportadas. 

A questao que surge, entao, e como obter, de forma determinfstica, a configuragao 
otima de carga e de trajeto da frota de caminhoes, de modo a se obter o menor custo 
na distribuigao de suprimentos da Marinha? 


1. Ballou (2006). 

2. Ballou (2008). 
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Para atacar o problema, foi usada a abordagem da Programacao Linear Inteira, 
cujo arcabouco teorico esta brevemente resumido a seguir. 

PROGRAMAgAO LINEAR 

A Programacao Linear (PL) tem sido bastante utilizada desde a decada de 1940 pelos 
militares e, ate hoje, muitas empresas se beneficiam do seu uso. Segundo pesquisa da 
revista Fortune, 85% das empresas ja utilizaram ou utilizam a Programacao Linear. A 
PL teve como um de seus inventores George B. Dantzig, que descobriu o algoritmo 
Simplex, mecanismo matematico utilizado para resolver os problemas da PL, a qual 
serviu de base para outras ferramentas gerenciais tais como a Programacao Nao Linear, 
a Programacao Inteira e outras tecnicas de otimizacao. Outro matematico, Leonid G. 
Khanchian, descobriu o algoritmo de tempo polinomial, que teoricamente e superior 
a algoritmos polinomiais, como o algoritmo Simplex. Outros algoritmos mais rapidos 
foram descobertos mais tarde, como o algoritmo de Karmarkar. 3 

Todo modelo de Programacao Linear tem duas caracterfsticas importantes: uma 
funcao-objetivo, a ser maximizada ou minimizada, e as restricoes. Esse modelo de 
otimizacao restrito fornece meios de solucionar problemas em que os recursos sao 
escassos de modo a otimizar um objetivo de interesse. 4 

O algoritmo Simplex utiliza uma serie de definicoes e conceitos associados, entre 
os quais a definicao de solucao basica, o conceito de variavel basica e nao basica e 
a definicao de solucao basica viavel. Uma solucao viavel e aquela que pode ser im- 
plementada considerando-se que os valores das variaveis de decisao estao de acordo 
com as restricoes. 5 

O Simplex e o mesmo algoritmo utilizado nos softwares especlficos para a solucao 
de problemas de PL e caminha de uma solucao viavel para outra, diminuindo o valor da 
funcao-objetivo ate o ponto otimo ser alcancado. Esse algoritmo e simples e resolve uma 
gama de problemas e modelos do mundo real, gerando uma grande economia de recursos. 6 

Para m equacoes e n variaveis existem no maximo n\/ml(n-m)l solucoes basicas. 
Ainda assim, o metodo Simplex geralmente converge para a solucao otima em apro- 
ximadamente 3m/2 iteracoes. 7 

Em relacao as restricoes da utilizacao do algoritmo Simplex, sao conhecidos os 
problemas de escala quando existem muitos calculos numericos e a variabilidade entre 
o maior e o menor numero e muito grande, causando problemas de calculo numerico. 
Outros problemas dizem respeito as: (i) solucoes inexistentes - quando nao ha solucao 
que satisfaca as restricoes; (ii) solucoes ilimitadas - quando o algoritmo nao consegue 
encontrar um ponto de maximo em um problema de maximizacao ou de minimo num 
problema de minimizacao; e (iii) solucoes multiplas - quando mais de uma solucao 
satisfaz todas as restricoes e alcanca um valor extremo da funcao-objetivo. 8 


3. Moore; Weatherford (2005). 

4. Moore; Weatherford (2005). 

5. Colin (2007). 

6. Colin (2007). 

7. Colin (2007). 

8. Colin (2007). 
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O algoritmo Simplex nao resolve qualquer problema de PL. Quando todas as 
variaveis sao estritamente positivas nas solugdes basicas viaveis, estamos diante de 
um problema conhecido como nao degenerado. Em outras palavras, se uma ou mais 
das variaveis basicas numa solugao tem o valor zero, a solugao e denominada solugao 
basica degenerada. A degeneragao traz consigo problemas de convergencia do algo¬ 
ritmo, sendo diffcil identificar a solugao otima. Geralmente, os softwares utilizados 
no mercado, tais como o Solver, ja tratam dessas particularidades, de forma invisfvel 
para o usuario. 9 

A Programagao Inteira (PI) surgiu como um caso particular da PL, quando na 
resolugao dos problemas havia a necessidade de uso de variaveis inteiras. O nfvel de 
dificuldade para resolugao de problemas de PI e similar ao da resolugao de problemas 
com PL, porem, enquanto estas possuem um espago infmito de solugoes, aquelas pos- 
suem um espago fmito de solugoes, o que implica que a solugao de PL seja mais facil 
que a solugao de PI. 10 

A PROGRAMAgAO DEVEICULOS SEM A PROGRAMACIAO INTEIRA 

Atualmente, a distribuigao da frota e realizada de modo empfrico, baseada na expe¬ 
rience e competencia dos profissionais envolvidos nessa tarefa. O objetivo da analise 
apresentada neste capftulo consiste em programar, de forma eficaz, a alocagao de 
vefculos frigorfficos do DepSubMRJ, atendendo as demandas semanais e minimizando 
os custos de distribuigao. Existe uma programagao semanal que atende tanto as OMs 
na Cidade do Rio de Janeiro quanto as cidades proximas. As entregas sao realizadas de 
segunda-feira a sexta-feira, inclusive feriados. Como restrigoes, ha uma programagao 
de demanda e percursos em distancia a percorrer para atender cada OM. A frota de 
vefculos frigorfficos e constitufda por caminhoes de capacidades diferenciadas (14t, 

8t e 3t), o que caracteriza a capacidade e a agilidade do caminhao no percurso. Os 
custos fixos e variaveis do transporte, a demanda semanal de cada OM e o valor 
maximo e mfnimo de entrega diaria levantados junto a instituigao foram apresentados 
nas Tabelas 1.1 a 1.3). 

A PROGRAMAgAO DEVEICULOS COM A PROGRAMAgAO 
INTEIRA 

Para resolver a questao proposta foi utilizada a ferramenta Solver do Excel, bus- 
cando-se obter a programagao que importasse menor custo possfvel. Como primeiro 
passo foi necessario identificar os fatores/variaveis que influenciavam diretamente no 
custo envolvido na distribuigao de suprimentos da Marinha. A partir dessa primeira 
analise, chegou-se, entao, a algumas premissas. 

A primeira diz respeito a composigao do custo total, que engloba o custo fixo, que 
nao se altera, e o custo variavel, que sofre influencia de acordo com a programagao 
executada. Neste caso, para diminuir o custo total, o que importa e fazer a distribuigao 



9. Colin (2007). 

10. Colin (2007). 
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nas OMs com a menor quilometragem possfvel, pois e isso que influencia diretamente 
o custo variavel. 

O segundo pressuposto indica que, para se obter menor quilometragem, deve haver 
o mmimo possfvel de viaturas deslocando-se para as OMs, a fim de satisfazer suas 
necessidades. No caso em que haja a opgao de suprir determinada OM, empregando 
uma ou duas viaturas, deve-se escolher a primeira alternativa, pois a segunda implica 
o dobro da quilometragem (um percurso por cada viatura). 

A terceira premissa preconiza que e desejavel suprir as OMs de uma so vez, 
pois, quanto mais viagens forem realizadas, maior sera a quilometragem. As me- 
lhores escolhas sao, portanto, aquelas em que se utiliza apenas um dia para suprir 
determinada OM. 

Finalmente, deve-se observar que, embora nao haja diferenga no calculo dos custos 
variaveis, deve haver penalizagao para a alocagao de viaturas com pouca carga. Afinal, 
quando isto ocorre, esta se deixando de aproveitar espago util nas viaturas, o que nao 
deixa de ser uma perda, ainda que implfcita. 

Em seguida, foi realizada a programagao no Excel, procurando-se encontrar uma 
resposta determinfstica para a questao de estudo, como alternativa para a programagao 
da frota, executada apenas com base na experiencia dos profissionais da Marinha. 

Primeira tentativa de solugao 

Baseando-se nas premissas listadas anteriormente, buscou-se a princfpio modelar 
o problema no Excel, criando uma tabela geral de variaveis ajustaveis, utilizando 
numeros inteiros, correspondendo a tonelagem alocada em cada viatura. Estas variaveis 
indicariam, diariamente, quais viaturas seriam empregadas, qual a capacidade utilizada 
das viaturas e que OMs seriam abastecidas. 

Para um custo mmimo, ou seja, menor quilometragem, foram estabelecidas algumas 
restricoes. Em primeiro lugar, diariamente, cada viatura nao poderia ser abastecida com 
uma tonelagem acima de sua capacidade. Era imperativo tambem que as demandas 
de cada OM deveriam ser atendidas, ao termino da semana. Por fim, as variaveis de 
decisao deveriam ser binarias. A Tabela 1.4 traz a configuragao do problema para a 
segunda-feira, considerando-se que para cada dia da semana o numero de celulas ajus¬ 
taveis era a mesma. Dessa forma, o valor 0 significaria a nao utilizagao da viatura, ao 
passo que o valor 1 indicaria o uso do vefculo. 

Apesar de a configuragao do problema obedecer a uma logica correta, ao se utilizar 
a ferramenta Solver do Excel nao foi possfvel se chegar a uma solugao, sendo exposta 
a seguinte mensagem: “Excesso de celulas ajustaveis”. Sem duvida, haveria grande 
necessidade de processamento dos calculos para se testar todas as variaveis na tabela 
geral. Por este motivo, esta que seria a solugao mais direta no Solver teve que ser 
descartada. 

Em diversas outras tentativas de resolugao do problema, as limitagoes de calculo 
do Solver ficaram evidentes. Em alguns casos a execugao do programa demorava horas 
sem que apresentasse uma solugao. Como resultado dessas experiences, chegou-se 
a conclusao de que a modelagem do problema deveria considerar fases de resolugao 
que nao excedessem a capacidade da ferramenta empregada. 
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TABELA 1.4 Extrato das celulas 

variaveis 

da tabela geral 




Segunda 


OM A 

OM B 

OM C 

OM D 

OM E 

OM F 

OM G 

Vtr 1 (14t) 

0 

0 

0 

0 


0 

0 

0 

Vtr 2 (14t) 

0 

0 

0 

0 


0 

0 

0 

Vtr 3 (14t) 

0 

0 

0 

0 


0 

0 

0 

Vtr 4 (14t) 

0 

0 

0 

0 


0 

0 

0 

Vtr 5 (8t) 

0 

0 

0 

0 


0 

0 

0 

Vtr 6 (8t) 

0 

0 

0 

0 


0 

0 

0 

Vtr 7 (8t) 

0 

0 

0 

0 


0 

0 

0 

Vtr 8 (3t) 

0 

0 

0 

0 


0 

0 

0 

Vtr 9 (3t) 

0 

0 

0 

0 


0 

0 

0 

Totais 

0 

0 

0 

0 


0 

0 

0 

V 








y 


Primeira solugao 

Levando-se em consideragao as premissas levantadas e a conclusao da necessidade 
de atacar o problema por fases, conforme anteriormente mencionado, chegou-se a 
primeira solugao determimstica com o uso da ferramenta Solver no Excel. Para essa 
primeira modelagem valida, alguns procedimentos foram adotados. Em primeiro 
lugar, foi criada uma planilha com parametros, para servir como base para as demais 
planilhas que seriam criadas. Alem disso, a solugao foi dividida em duas etapas, sendo 
empregada a ferramenta Solver em cada uma delas. 

Na primeira etapa, para cada OM foi realizado um modelo Solver individual, 
com o proposito de verificar qual seria a melhor alocagao de viaturas para satisfazer 
a demanda de uma unica vez. A celula objetivo priorizou a menor quilometragem 
possfvel, com a menor ociosidade possfvel. Por este motivo, sua composigao foi de 
viaturas empregadas, multiplicada pela distancia, mais o espago ocioso. Com isto, 
procurou-se tentar atender as restrigoes preestabelecidas. 

Na segunda etapa, foi executado um Solver geral para decidir quais OMs seriam 
atendidas em cada dia da semana. Para tal, tomou-se como base o resultado da etapa 
anterior, que alocou a melhor quantidade de viaturas para satisfazer as OMs de uma 
unica vez. Neste Solver foram impostas as restrigoes de quantidade de viatura e de que 
cada OM nao poderia ser atendida mais de uma vez na semana. O mais importante na 
segunda etapa foi satisfazer as restrigoes. Por isso, a media foi um artiffcio utilizado 
para criar a variavel objetivo; afinal, qualquer solugao apresentaria o mesmo custo, pois 
a quilometragem seria sempre a mesma. Inicialmente, a variavel-padrao foi definida. 
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tomando-se como base o desvio-padrao, porem, isso dificultou o trabalho do Excel 
para apresentar uma solu§ao valida. 

Cabe ressaltar que nao foram utilizados numeros necessariamente binarios para 
as variaveis de decisao na segunda etapa. No entanto, na verdade os numeros sempre 
eram 0 ou 1, devido a restr^ao de se atender a OM apenas uma vez na semana. As 
Tabelas 1.5 a 1.7 trazem o detalhamento de uma das configurates validas encontradas 
nesta primeira solu 9 ao determinlstica. 

Segunda solugao 

A primeira solu£ao nao atendia o requisito de entregas maxima e minima impostas pelo 
problema. Embora isso nao fosse um prejulzo consideravel, por nao implicar maiores 
custos ou indisponibilidade de viatura, foi necessario elaborar uma segunda solutjao, 
que satisfizesse a todos os requisitos. 


TABELA 1.5 OMs atendidas 

por dia 




N 

OM 

Seg 

Ter 

Qua 

Qui 

Sex 


Alfa 

0 

1 

0 

0 

0 


Bravo 

0 

0 

1 

0 

0 


Charlie 

0 

1 

0 

0 

0 


Delta 

1 

0 

0 

0 

0 


Echo 

1 

0 

0 

0 

0 


Foxtrot 

0 

1 

0 

0 

0 


Golf 

1 

0 

0 

0 

0 


Hotel 

1 

0 

0 

0 

0 


India 

0 

0 

0 

0 

1 


Juliet 

0 

0 

0 

1 

0 


Entrega 

71 

50 

15 

13 

16 
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TABELA 1.7 Distribuigao de viaturas e quilometragens envolvidas 


OM 

Vtr 14t 

Vtr 8t 

Vtr 3t 

Demanda 

Distancia 

Km por OM 

Alfa 

0 

1 

0 

8 

12 

12 

Bravo 

0 

2 

0 

15 

5 

10 

Charlie 

0 

2 

0 

16 

140 

280 

Delta 

1 

1 

0 

18 

160 

320 

Echo 

1 

1 

0 

20 

15 

30 

Foxtrot 

2 

0 

0 

26 

10 

20 

Golf 

1 

1 

0 

19 

180 

360 

Hotel 

1 

0 

0 

14 

6 

6 

India 

0 

2 

0 

16 

25 

50 

Juliet 

1 

0 

0 

13 

20 

20 

Total 






1.108 

V._ 






y 


Cabe ressaltar que a segunda solugao e mais trabalhosa que a primeira, embora 
seja a mais completa e a unica que realmente satisfaz o problema proposto. Foram, 
portanto, adotados alguns procedimentos. Em primeiro lugar, criou-se uma planilha 
de parametros, de forma semelhante a anterior. Em seguida, dividiu-se a resolugao 
em tres fases. 

Para a primeira fase partiu-se da premissa de que a melhor coisa a ser feita e suprir 
o maximo de OMs o mais rapido possfvel, obedecendo, e claro, as restrigoes de entrega 
maxima e minima. Alem disso, a OM deve ser suprida em apenas um unico dia, para 
minimizar a quilometragem. Por este motivo, foi criado um Solver diario, para decidir 
quais OMs seriam supridas diariamente. 

A variavel-objetivo foi a sobra de toneladas a serem supridas, devendo ser a minima 
possivel. Cabe ressaltar que as planilhas estavam relacionadas de forma que a demanda 
considerada em determinado dia fosse igual a demanda que deixou de ser atendida no 
dia anterior, de modo a nao se atender mais de uma vez a determinada localidade. A 
Tabela 1.8 traz o resultado do Solver para segunda-feira. 

Com os dados da primeira fase chegou-se a uma tabela que demonstrava como as 
OMs seriam atendidas, ao longo da semana, conforme exposto na Tabela 1.9. Faltava, 
no entanto, determinar como fazer a distribuigao das viaturas para que fosse atendido 
o objetivo maior do problema: determinar a melhor programagao para reduzir custos. 
Para tal, foi realizada a segunda fase. 

Como base para a segunda fase, alem da Tabela 1.9, foram adotados na varia¬ 
vel-objetivo os principios de menor quilometragem possivel e tambem o de menor 
capacidade ociosa possivel (espago vazio em viatura utilizada). Estes sao os mesmos 
parametros utilizados na solugao anterior. 
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TABELA 1.8 Primeira fase (segunda-feira) 

Entrega (max) 

45 

45 




Entrega (min) 

31 

45 

Efetuada 

Resta 


Alfa 

8 

0 

0 

8 


Bravo 

15 

0 

0 

15 


Charlie 

16 

0 

0 

16 


Delta 

18 

0 

0 

18 


Echo 

20 

0 

0 

20 


Foxtrot 

26 

1 

26 

0 


Golf 

19 

1 

19 

0 


Hotel 

14 

0 

0 

14 


India 

16 

0 

0 

16 


Juliet 

13 

0 

0 

13 


Total 

165 



120 


s -- 





_' 


TABELA 1.9 Consolidate) da primeira fase 



"N 

OM 

Seg 

Ter 

Qua 

Qui 

Sex 

Dist 

Alfa 

0 

0 

1 

0 

0 

12 

Bravo 

0 

0 

1 

0 

0 

5 

Charlie 

0 

0 

0 

1 

0 

140 

Delta 

0 

1 

0 

0 

0 

160 

Echo 

0 

0 

1 

0 

0 

15 

Foxtrot 

1 

0 

0 

0 

0 

10 

Golf 

1 

0 

0 

0 

0 

180 

Hotel 

0 

1 

0 

0 

0 

6 

India 

0 

0 

0 

1 

0 

25 

Juliet 

0 

0 

0 

0 

1 

20 

Totais 

45 

32 

43 

32 

13 



V- 
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TABELA 1.10 Segunda fase (segunda-feira) 



"N 


Vtr (t) 



Dist envo 





14 

8 

3 


Demanda 

Dem at 

Dem ref 


0 

0 

0 

0 

8 

0 

0 


0 

0 

0 

0 

15 

0 

0 


0 

0 

0 

0 

16 

0 

0 


0 

0 

0 

0 

18 

0 

0 


0 

0 

0 

0 

20 

0 

0 


2 

0 

0 

20 

26 

28 

26 


1 

1 

0 

360 

19 

22 

19 


0 

0 

0 

0 

14 

0 

0 


0 

0 

0 

0 

16 

0 

0 


0 

0 

0 

0 

13 

0 

0 

Utl 

3 

1 

0 

< = 

Cap total 

50 

45 

Disp 

4 

3 

2 





s '- 







_^ 


As restricdes adotadas em cada Solver desta segunda fase referiam-se, sobretudo, 
a limitagao de viatura. A Tabela 1.10 traz o resultado do Solver de segunda-feira, nesta 
segunda fase: “Dem at” refere-se a capacidade das viaturas empregadas, enquanto 
“Dem ref’ diz respeito a necessidade das OMs. 

Na ultima fase dessa resolugao, foi feita apenas a elaboragao de uma planilha de 
consolidagao, com o resultado final da primeira e da segunda fases da resolugao. A 
terceira fase detalha, portanto, diariamente, o percurso das viaturas que serao em¬ 
pregadas na distribuigao de suprimentos pelas OMs. 

Compara^ao entre as solu^oes 

A unica vantagem da primeira solugao em relagao a segunda e ser relativamente mais 
simples, exigindo menos tabelas Solver. Esta aparente vantagem, no entanto, talvez 
nem mesmo possa ser considerada, tendo em vista que somente na segunda solugao 
foram atendidos os requisitos propostos pelo problema de entrega maxima e minima. 

Alem disso, a segunda solugao utiliza de modo mais eficaz o espago das viaturas. 
Isto acontece porque ela nao considera as OMs de forma isolada, e sim a entrega que 
deve ser efetuada no dia. Evita-se, com isso, a sobreposigao de sobras de espagos vazios 
nas viaturas quando estas sao alocadas, considerando-se as OMs de forma isolada. 
Como resultado, pode-se, inclusive, utilizar viaturas de tres toneladas, o que nao 
ocorria na primeira solugao. 
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Finalmente, e possfvel concluir que a segunda solugao e, com seguran 9 a, a otima, 
porque esta em consonancia com o objetivo maior da ferramenta Solver, que e apontar 
a melhor solugao. No caso da primeira solugao, podem existir alternativas melhores, 
e se assim for o instrumento de apoio a decisao esta comprometido. Por isso, embora 
mais trabalhosa, a segunda solugao e a que deve ser considerada na solugao do pro- 
blema proposto. 

CONSIDERAgOES FINAIS 

A resolugao do problema real com o uso do Excel permitiu que se chegasse a uma 
visao acerca das potencialidades e das limitagoes da ferramenta Solver. Verificou-se, 
nesse sentido, a necessidade de executar a modelagem do problema, considerando o 
processamento eletronico que possui limitagoes. No entanto, apos tratar o problema 
em fases, a ferramenta Solver pode resolver questoes complexas, como a considerada 
neste trabalho. 

Em relagao as limitagoes da ferramenta Solver reportadas na literatura, tomou-se 
cuidado em relagao ao problema mencionado de escala para que nao houvesse muita 
variabilidade entre o maior e o menor numero nos calculos efetuados, o que poderia 
gerar problemas de convergencia. 11 

Outra limitagao da ferramenta Solver que deve ser mencionada diz respeito ao 
problema da degeneragao. Conforme reportado na literatura, 12 se uma ou mais das 
variaveis basicas numa solugao tern o valor zero, a solugao e chamada solugao basica 
degenerada, e isto implica problemas de convergencia do algoritmo, sendo diflcil 
identificar a solugao otima. Ainda que o algoritmo Simplex trate esses problemas, este 
fator deve ser considerado na tentativa de resolugao na primeira solugao. 

Entretanto, na primeira solugao, para atender os requisites preestabelecidos e apre- 
sentar um valor para a celula alvo que fosse otimo, o Solver apresentou, inicialmente, 
em alguns casos, uma demora excessiva. Conforme reportado na literatura, 13 para m 
equagoes e n variaveis deveria haver, no maximo, n\/m\(n-m)\ solu 9 oes basicas e o 
algoritmo deveria convergir para a solu 9 ao otima em aproximadamente 3m/2 itera 9 oes. 
Trazendo essas formulas para o problema em questao, verifica-se que, no primeiro 
caso, existiam 362 equa 9 oes relacionadas as restri 9 oes e 315 variaveis, o que poderia 
levar a 543 itera 9 oes, explicando a demora na sohnjao e a nao convergencia para um 
valor otimo. Por outro lado, com a solu 9 ao vislumbrada, utilizando-se a defasagem 
dos dias da semana, conseguiu-se uma convergencia com 18 itera 9 oes. 

Durante a resolu 9 ao do problema, constatou-se que a formula 9 ao de um bom 
modelo e crucial para a sua solu 9 ao e que e importante a analise crftica detalhada para 
avalia 9 ao da real necessidade de se modelar uma variavel como inteira. Tal constata 9 ao 
foi baseada na falta de convergencia do resultado na primeira solu 9 ao e no sucesso 
alcan 9 ado com a nova modelagem na segunda solu 9 ao, considerando-se a defasagem 
proposta. Cabe considerar que, embora trabalhosa a montagem do problema, apos os 


11. Colin (2007). 

12. Colin (2007). 

13. Colin (2007). 
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dados serem correlacionados na planilha nao existem dificuldades para se efetuar novas 
solucdes a partir de modificagoes que venham a ocorrer. No caso, por exemplo, de se 
acrescerem viaturas ou haver aumento de demanda, a partir da mudanga na planilha de 
parametros todos os demais calculos serao efetuados de forma automatica. A solugao 
apresentada nao serve, portanto, apenas para responder ao momento atual, podendo 
ser utilizada tambem em ocasioes futuras. 

Outro aspecto importante que deve ser salientado e que nas solugoes de problemas 
reais desta natureza nao devem ser considerados apenas os custos e parametros ex- 
plfcitos. Trafegar com viaturas com espago ocioso, matematicamente, nao apresentaria 
aumento de custos. Conforme mencionado na literatura, 14 na logfstica o custo total 
e dependente nao so da variavel custo do transporte, mas tambem do custo de es- 
tocagem. No entanto, sabe-se que, ao menos indiretamente, isto representa perda, pelo 
espago ocioso. De forma semelhante, diante de um mesmo custo, e melhor atender 
o mais rapido possfvel a demanda do usuario. Cabe, portanto, a quem vai modelar o 
problema, a tarefa de considerar estes e outros aspectos implfcitos na resolugao de 
questoes reais. 

Como indicagao para pesquisas futuras, fica a sugestao de se tentar modelar tais 
problemas de modo ainda mais realista, considerando todos os aspectos mais impor- 
tantes em relagao a situagao em foco. Por exemplo, a demanda de cada OM poderia 
ser considerada de forma probabilfstica, atraves de uma distribuigao de probabilidades 
condizente com o seu historico de comportamento, o que estaria mais de acordo com 
a realidade. Isso seria possfvel por meio do uso de uma abordagem probabilfstica, 
como a Simulagao de Monte Carlo, 15 em vez da Programagao Inteira, utilizada neste 
trabalho, que considerou a demanda de forma puramente determinfstica. 

Cabe ressaltar, no entanto, que, no caso de modelagens ainda mais realistas (e, 
portanto, potencialmente mais complexas), tornam-se ainda mais importantes os 
cuidados em relagao a modelagem mencionados nestas consideragoes finais, para 
que o esforgo computacional na busca da solugao otima nao ultrapasse os limites da 
viabilidade operacional. 

EXERCICIOS 

1. Uma empresa tern que planejar seus gastos em Pesquisa e Desenvolvimento para 
os proximos cinco anos. A empresa pre-selecionou quatro projetos e deve escolher 
dentre estes quais priorizar. Os dados relevantes ao problema encontram-se na 
planilha Exercfcio 1 . Nela tambem se encontra a disponibilidade de capital a 
ser alocado em cada um dos anos, bem como o valor presente lfquido (VPL) de 
cada projeto. Como todos os projetos apresentam VPL positivos, todos seriam 
candidatos. Vale notar que existe uma limitagao no valor a ser investido anualmente. 

a. Quais projetos devem ser priorizados? 

b. Agora considere que os projetos 1, 3 e 4 sao concorrentes e apenas um deles 
pode ser escolhido. Quais projetos devem ser priorizados? 



14. Ballou (2008). 

15. Hertz (1980); Paula et al. (2007). 
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c. Sem considerar a letra b, dessa vez considere que o projeto 1 depende de uma 
tecnologia que deve ser desenvolvida pelo projeto 3, isto e, so pode ser aprovado 
se o projeto 3 tambem o for. Quais projetos devem ser priorizados? 

2. Um cliente procurou um analista de investimentos com o objetivo de investir em 
agoes. O analista o entrevista para identificar o seu perfil de risco e apresentar as 
opcbes de investimento. O cliente esta decidido a investir em papeis de bancos. O 
analista deve apresentar uma proposta de investimento que atenda as expectativas 
do cliente. O retorno e o risco de cada investimento estao apresentados na planilha. 
Exerclcio 2. Durante a entrevista, o cliente fez as seguintes exigencias: 

• pelo menos um quarto da carteira deve ser composta por afoes do Banco A, 
devido ao seu baixo risco; 

• a participate das a§oes do Banco B deve ser de, no maximo, 12,5%, devido 
ao seu alto risco; 

• a participato de Renda Fixa na carteira deve ser de, no maximo, 10%, devido 
ao seu baixo retorno; 

• a taxa maxima suportada de risco e de 5,5%. 

3. Uma corretora de valores gerencia recursos de terceiros atraves da composhjao de 
carteiras de investimento para clientes, com diferentes perfis de risco. Um deles 
exige que: 

• nao mais de 30% do total aplicado seja aplicado em um unico investimento; 

• um valor nao inferior a 40% do total aplicado seja investido em tltulos de 
maturidade maior do que 10 anos; 

• o total aplicado em tltulos de alto risco seja, no maximo, de 40% do total inves¬ 
tido. 

A planilha Exerclcio 3 mostra os dados dos tltulos selecionados. Determine qual per- 
centual do total deve ser aplicado em cada tltulo. 

4. Um indivlduo dispoe de R$ 6.000 para investir em ate dois ativos. Ele tem tres 
op 5 oes de ativos, que vao fornecer o retorno apos cinco anos. Cada um dos inves¬ 
timentos pode ser feito em valores multiplos de R$ 1.000 e ele deseja investir, no 
maximo, R$ 4.000 em cada ativo. Os retornos esperados apresentam ganhos de 
escala (ou seja, sempre e melhor fazer uma aplicado de R$ 3.000 em um ativo do 
que uma aplica§ao de R$ 1.000 e outra de R$ 2.000 no mesmo ativo, por exemplo) 
e podem ser visualizados na planilha Exerclcio 4. Qual deve ser a polltica de 
investimento? 

5. Uma empresa possui as seguintes op 9 oes de investimento em projetos para os 
proximos quatro anos: investimento mac^o em P&D, polltica agressiva de Vendas, 
abertura de uma nova loja e recapacita 9 ao dos Recursos Humanos. Enfrentando 
limita 9 oes anuais de capital, a empresa deve escolher em quais projetos investir 
para obter o maior retorno posslvel, a valor presente. O VPL dos projetos, o capital 
necessario para cada um deles e a proje 9 ao do capital disponlvel estao ilustrados 
na planilha Exerclcio 5. 
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a. Determine que projetos a companhia deve selecionar para maximizar o valor 
presente lfquido do capital investido. 

b. Considere que, devido a questoes estrategicas, a empresa decida que a recapaci- 
tagao dos Recursos Humanos e um projeto que nao pode deixar de ser realizado, 
independentemente de seu VPL estimado. Que mudangas esta informagao gera 
na formulagao do problema e no resultado? 

c. Sem considerar o que foi mencionado na letra b, o que voce recomendaria se a 
companhia conseguir um adicional de R$ 10.000 para cada um dos quatro anos? 

d. E considerando o que foi mencionado na letra b? 
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Objetivos de aprendizagem do capi'tulo 

• Saber modelar um problema real em situagao de risco atraves de uma Arvore de 
Decisao. 

• Entender o conceito e as diferengas entre no de decisao e de evento. 

• Compreender o conceito do Valor Monetario Esperado. 

• Aprender a estimar o Valor da Informagao Imperfeita e o Valor da Informagao Perfeita. 

• Saber conduzir uma Analise de SensibiIidade acerca dos parametros mais importantes. 

• Conseguir considerar, na analise, o perfil de aversao a risco do tomador de decisao. 


INTRODUgAO 

Uma decisao precisa ser tomada quando se esta diante de um problema que dispoe de 
mais de uma alternativa para sua solugao. A Teoria da Decisao parte do pressuposto 
de que os individuos sao capazes de expressar suas preferencias basicas, e sao racio- 
nais, quando enfrentam situagoes de decisoes simples. Com base nessa proposigao, a 
metodologia desenvolvida pela Teoria da Decisao possibilita a resolugao de problemas 


1. Gomes; Gomes; Almeida (2002). 
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de decisao mais complexos. Estes problemas sao comuns em uma infinidade de areas, 
tanto publicas quanto privadas, e desde tempos remotos o homem tenta resolve-los 
apoiando-se em abstragoes, heurfstica e raciocfnios dedutivos, a fim de validar e guiar 
escolhas. Para resolver tais problemas, ate a primeira metade do seculo XX, utilizava-se 
basicamente a esperanga matematica para tomada de decisao em condigoes aleatorias; 
porem, em muitas situagoes, o risco associado a tal procedimento era inaceitavel. 

No processo de modelagem da realidade, o pesquisador troca a riqueza e a abran- 
gencia da realidade pelo poder de analise e pela capacidade de experimentagao. Neste 
caso, os elementos do modelo a ser formulado podem ser vistos como variaveis de 
controle ou de decisao e ser manipulados para se atingir os objetivos. O modelo com- 
preende ainda variaveis de estado ou de natureza, sobre as quais nao se tern controle, 
mas que afetam as consequencias ou resultados da decisao. 

O ambiente atual dentro e fora das organizagoes e dinamico e cheio de incerteza. Os 
gerentes estao sobrecarregados devido a complexidade dos processos organizational. As 
exigencias diarias em relagao a decisoes rapidas e assertivas demandam o uso de ferramen- 
tas de tomada de decisao que possam aumentar a qualidade e a velocidade nesse processo. 

O processo de tomada de decisao dentro das organizagoes modernas demanda um 
profissional que conhega o ambiente no qual a empresa esta inserida e que tambem 
esteja ciente das suas constantes e frequentes mudangas. 

Em constante cobranga por eficiencia, o profissional deve dispor de conheci- 
mentos tecnicos sobre o assunto em relagao ao qual vai opinar e tambem ser agil no 
levantamento de todas as alternativas possfveis. Para atender essa necessidade, o uso 
dos metodos quantitativos evoluiu para o conceito da Teoria da Decisao, dentro da 
qual se encontram diversas ferramentas e tecnicas que contribuem para maior eficacia 
gerencial. Os metodos quantitativos modernos desenvolveram a abordagem logica 
para situagoes de incerteza e risco, o que se aplica ao ambiente atual das organizagoes. 

E diante de uma dessas situagoes de risco que se pode recorrer a uma ferramenta 
que permita, a partir de probabilidades conhecidas ou subjetivamente atribuidas e 
com o auxflio de heurfsticas, visualizar e interpretar facilmente algumas opgoes para 
a tomada de decisao mais acertada. 

ARVORE DE DECISAO 

Na decada de 1970, comegaram a surgir os primeiros metodos de apoio com o intuito de 
enfrentar situagoes espetificas, nas quais um tomador de decisao, atuando com racionali- 
dade, deveria resolver um problema em que varios eram os objetivos a serem alcangados 
de forma simultanea. Dal surgiu a Arvore de Decisao como um dos mais importantes 
metodos de analise de decisao, pautando-se em conceitos como o de risco e retorno. 

As Arvores de Decisao utilizam a estrategia clividir-e-conquistar ( divide-and- 
conquer ), em que arvores sao construldas utilizando-se apenas alguns atributos. 
As Arvores de Decisao sao uma das tecnicas de aprendizado de maquina (machine 
learning), na qual um problema complexo e decomposto em subproblemas mais sim¬ 
ples. Recursivamente, a mesma estrategia e aplicada a cada subproblema. 2 


2. Gama (2004). 
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Quinlan, professor da Universidade de Sidney, considerado o “pai das Arvores 
de Decisao”, elaborou um novo algoritmo, chamado ID3, desenvolvido em 1983. O 
ID3 e suas evolugoes (11)4. ID6, C4.5 , See 5) sao algoritmos muito utilizados para 
gerar Arvores de Decisao. O atributo mais importante e apresentado na arvore como 
o primeiro no, e os atributos menos importantes, segundo o criterio utilizado, sao 
mostrados nos nos subsequentes. As principals vantagens das Arvores de Decisao 
sao o fato de que elas “tomam decisoes” levando em consideragao os atributos que 
sao considerados mais relevantes, segundo a metrica escolhida, alem de serem com- 
preensfveis para as pessoas. Ao escolher e apresentar os atributos em ordem de 
importancia, as Arvores de Decisao permitem aos usuarios conhecer quais fatores 
mais influenciam os seus trabalhos. 

Dentro do conceito do metodo de Arvore de Decisao esta o modelo de arvore 
binomial, com o qual se pode descrever uma multiplicidade de alternativas numa 
escolha de opgao otimizada. A tarefa deste metodo e estabelecer uma relagao entre 
a situagao de tomada de decisao e uma expressao numerica das opgoes disponfveis. 

Um problema representa o infcio de uma arvore que se compoe de ramos, nos e 
folhas. Para cada no deve ser definido o proximo criterio de ramificagao, uma decisao 
ou cenario. Os ramos representam decisoes e cenarios, e as folhas sao alternativas 
finais para a escolha de uma unica e melhor opgao. Ao receber novos dados, a arvore 
modifica-se e gera uma nova decisao otimizada. 

O problema em construir uma Arvore de Decisao pode ser expresso recursivamente: 
primeiro deve-se selecionar um atributo para colocar no no-raiz e fazer um ramo para 
cada possfvel valor. Isso divide o problema em subconjuntos, um para cada valor do 
atributo. Agora o processo pode ser repetido recursivamente para cada ramo. Se a 
qualquer instante todos os exemplos em um no tern a mesma classificagao, interrompe-se 
o desenvolvimento daquela parte da arvore. Como determinar, no entanto, qual atributo 
dividir? Escolhe-se o atributo que gere uma arvore menor e que tenha chances de clas- 
sificar melhor, ou seja, precisa-se mensurar o grau de pureza de cada no. Com isto, e 
possfvel escolher o atributo que produz os nos filhos mais puros. 3 

A Arvore de Decisao esta diretamente ligada com o que se denomina flexibilidade 
gerencial, pois o metodo de analise permite rever a estrategia inicial e alterar pianos de 
acordo com novas condigoes economicas. As flexibilidades gerenciais possibilitam tanto 
capitalizar futuras oportunidades favoraveis ao negocio quanto diminuir perdas, isto e, 
melhoram o potencial de ganhos e limitam o potencial de perdas. 4 A flexibilidade geren¬ 
cial consiste em um conjunto de opgoes reais tais como postergar um projeto, expandir 
ou contrair escala de produgao ou avaliar a alteragao das materias-primas de um projeto. 

Considera-se a flexibilidade gerencial como uma possibilidade, nao como uma 
obrigagao. Muitos autores acreditam que as Opgoes Reais no futuro nao constituirao 
meramente uma ferramenta analftica mas, em contextos mais rigorosos e intensivos, 
um processo organizacional. Contudo, a Arvore de Decisao e um metodo que permite 
indicar, de forma grafica e cronologica, um caminho a ser seguido em um processo de 
decisao, explicitando etapas a serem cumpridas para alcangar o objetivo pretendido. 



3. Carvalho (2002). 

4. Minardi (2004). 
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A arvore representa um processo de decisao em que os nos sao os momentos no 
tempo em que o decisor deve efetuar uma decisao. As etapas da tecnica para o auxflio da 
tomada de decisao sao definidas como: a) defmigao do tema; b) defmigao do objetivo, 
metas e submetas; c) construgao da Arvore de Decisao; d) revisao da Arvore de Decisao; 
e) encerramento. 5 

Algumas vantagens da utilizagao desta tecnica podem ser relacionadas: 6 

• permite a subdivisao do objetivo em metas e submetas, indicando como alcanga-las; 

• orienta o tomador de decisao a medida que este responde a pergunta: O que e 
necessario fazer para alcangar a meta pretendida?', 

• permite o exame, pelo tomador de decisao, de todas as possibilidades; 

• permite a criagao de algoritmos facilmente implementados em computadores. 

Uma decisao tomada com base em todos os dados e informagoes disponiveis e alter- 
nativas possfveis e que, principalmente, utilize ferramentas de metodos quantitativos 
apropriados, deve ser considerada como uma boa decisao, mesmo que apresente algo 
que um administrador nao vise como meta. Complementando essa convicgao, uma 
decisao que nao considerou todos os procedimentos prescritos pela metodologia e 
considerada ma decisao, mesmo que o resultado obtido seja favoravel. Isso porque, 
como nao foram exploradas todas as informagoes e alternativas possfveis, nao se pode 
afirmar que o resultado favoravel foi otimizado. 7 

Vale destacar que as Arvores de Decisao garantem maior flexibilidade do que os 
mecanismos tradicionais na medida em que as decisoes sao definidas, tal como fora 
mencionado anteriormente, com maior grau de flexibilidade. Nas arvores, uma serie de 
eventos podem ser mapeados ao longo dos diversos ramos, envolvendo varias decisoes 
de sequenciamento. A Figura 2.1 apresenta uma Arvore de Decisao simplificada. 

Constata-se com base na Figura 2.1 que os valores de “p” referem-se as proba- 
bilidades de ocorrencia de um evento. A Arvore de Decisao incorpora conceitos de 
incerteza na analise de uma alternativa de investimento ou mesmo de financiamentos, 
o que a torna uma tecnica superior as tecnicas tradicionais. 



5. Gomes; Gomes; Almeida (2009). 

6. Gomes; Gomes; Almeida (2009). 

7. Corrar; Theophilo (2004). 
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A ESCOLHA DE UM FORNECEDOR COM O AUXILIO 
DA ARVORE DE DECISAO 

Diante das necessidades operacionais de um gerente da importadora brasileira Meta¬ 
lurgica e Distribuidora de Ago, o exemplo a seguir ordena diretamente a utilizagao de 
uma ferramenta que visa auxilia-lo na tomada de decisao. O objetivo e determinar a 
melhor escolha entre tres fornecedores de ago para importagao pautando-se na tecnica 
da Arvore de Decisao. A determinagao do fornecedor torna-se muito importante 
porque existem diversas variaveis que, caso nao sejam cumpridas, vao comprometer 
a imagem da empresa distribuidora perante o mercado. As variaveis nao sao somente 
os custos, mas tambem o prazo de entrega, a qualidade de produto para aceitagao do 
mercado importador e outros, que sao fatores determinantes na escolha do fornecedor. 

A Arvore de Decisao e um metodo dinamico, facil de utilizar e pode ser aproveitado 
para futuras consultas, criando assim um banco de dados das decisoes. O desafio e 
aproximar os dados descritivos a uma analise quantitativa e evitar um julgamento 
intuitivo na tomada de decisao. 

A analise do processo de decisao se desdobra em duas etapas: a primeira e a decisao 
de importar ou nao, e a segunda refere-se a escolha entre tres fornecedores estrangeiros, 
com exploragao previa de atributos de avaliagao. 

A Teoria da Decisao estabelece procedimentos para um processo de tomada de 
decisao. Os passos a seguir analisam o problema proposto a luz dos conceitos desse 
processo, procurando, de forma didatica, determinar a melhor altemativa a ser seguida. 8 

Passo 1: Defini^ao do problema 

A Metalurgica e Distribuidora de Ago e uma empresa brasileira de grande porte que usa 
como insumo tarugos, chapas e bobinas de ago para produgao de tubos, cantoneiras, 
perfis e outros produtos; alem disso, e distribuidora desses produtos para metalurgicas 
de menor porte. A diretoria da empresa observou a tendencia de aquecimento da 
demanda para seus produtos devido a alguns fatores. 

Os produtos sao comprados de fabricantes de ago nacionais que, devido a grande 
demanda para consumo proprio, suprem as necessidades das metalurgicas a partir 
das quotas preestabelecidas. Dessa maneira fica impossfvel aumentar os volumes 
de compra de imediato. Caso o aquecimento da demanda se concretize, tal como 
indicant varias pesquisas das associagoes comerciais do setor metalurgico brasileira, 
a Metalurgica e Distribuidora de Ago precisara aumentar seu estoque, devido a sua 
necessidade, passando de 1 mil toneladas por mes para 4 mil toneladas para o mesmo 
perfodo. O problema mais evidente no primeiro momento e onde obter a quantidade 
de insumos alem da quota com fornecedor nacional. 

Passo 2: Alternativas de a^ao 

Sendo impossfvel obter o volume adicional de 3 mil toneladas de ago no mercado 
nacional, a diretoria tern algumas alternativas, que se alinham em duas etapas de 



8. Corrar; Theophilo (2004). 
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tomada de decisao, sendo a segunda somente necessaria caso opte positivamente pela 
primeira alternativa. 

Primeira etapa de decisao: 

1. importar o produto do exterior, avaliando com otimismo a situagao do mercado em 
expansao; ou 

2. permanecer no mesmo volume de 1 mil toneladas mensais, sem recorrer a compra 
adicional via importagao, devido as incertezas de crescimento do mercado de cons- 
trugao. 

Segunda etapa de decisao: 

1. Caso seja tomada a decisao de importar, deve-se escolher um fomecedor estrangeiro 
dentro de varias ofertas. 

O desafio do profissional que toma a decisao e precisar ou obter dados quanto a 
probabilidade de ocorrencia desse cenario. Essas probabilidades sao estimadas em 
termos percentuais. Entre os fatores que identificam o ambiente de incerteza na tomada 
de decisao do problema estudado estao aumento da demanda no mercado interno, 
mudangas na politica industrial no Brasil ou no pals exportador e volatilidade da taxa 
de cambio de dolar norte-americano, entre outros. 

Baseando-se em dados do mercado, o profissional conseguiu atribuir as proba¬ 
bilidades de ocorrencia de cenarios, concluindo que a decisao seja tomada em um 
ambiente de risco. A condigao do cenario de aquecimento do mercado configura 
um exemplo de opgao real, em que a decisao de importar pode ser tomada apos a 
revelagao da informagao. Mas, no caso analisado, o distribuidor deseja estar preparado 
para o aumento da demanda devido aos prazos extensos de entrega da mercadoria na 
importagao e para estar com o produto a disposigao de compradores. 

Para o tomador da decisao, a situagao ideal seria aquela em que a probabilidade 
de ocorrer um cenario fosse estimada em 100%. Nessa situagao, pode-se considerar 
que o tomador de decisao esta diante de uma informagao perfeita, que transformaria 
um ambiente de risco em um ambiente de certeza. Com base em informagoes sobre 
mercado previamente coletadas, o dirigente admite que as probabilidades de ocorrencia 
dos cenarios de primeira etapa sejam consideradas como as da Tabela 2.1. 

Em se tratando de recusa a importagao de um estoque adicional, nao ha necessidade 
de proxima etapa de decisao, ou seja, de segundo nlvel. Mas, ao optar pela importagao 
do produto, a empresa encontra alguns riscos conjunturais do pals exportador e riscos 


TABELA 2.1 Cenarios e suas probabilidades de ocorrencia 


Cenarios para a primeira etapa 
de decisao: Importar 


Probabilidades 


Cenario! (C,): Aquecimento do mercado 


80 % 


Cenario 2 (C 2 ): Mercado Estavel 


20 % 
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proprios do negocio a ser efetuado. Tambem existe um grau de incerteza nas avalia 9 oes 
feitas em rela 9 ao aos fornecedores, o que significa que um fornecedor apenas bom 
pode ter sido avaliado como otimo e vice-versa. 

Um cenario de risco conjuntural presente no trabalho e a situa 9 ao economica do se- 
tor siderurgico na China. Devido a crescente produ 9 ao e demanda do a 90 naquele pais, 
o governo pode intervir a qualquer momento, limitando a produ 9 ao de 390 nacional 
com o intuito de nao desvalorizar o pre 90 das exporta 9 oes chinesas. A probabilidade 
atribuida pelo profissional a ocorrencia dessa a 9 ao e de 60%. 

Portanto, antes de optar pelo fornecedor chines, deve-se observar a probabilidade de 
que este, diante de novas politicas industrials, podera nao cumprir o contrato de venda, 
e o importador brasileiro ficara sem o produto para a data acordada. Tal ocorrencia 
tambem deve ser analisada do ponto de vista fmanceiro, mesmo nao incorrendo em 
gastos, pois pode acarretar custos como coloca 9 ao de pedido, abertura de carta de 
credito e, ainda, o prejuizo a imagem da empresa por nao dispor do produto conforme 
demanda o mercado, o que representa uma limita 9 ao para esse estudo, ja que nao sera 
levada em conta. 

O risco conjuntural para fornecedores franceses e ucranianos deve ser estimado em 
10%, pois sempre ha uma margem de risco. Assim, baseando-se na probabilidade de 
ocorrencia de interven 9 ao do governo, a Tabela 2.2 mostra a probabilidade de sucesso 
e insucesso na entrega de mercadoria. 


Passo 3: Possfveis resultados decorrentes das alternativas 


Cada alternativa de a 9 ao tera um ganho e uma perda determinados pelo sucesso ou 
insucesso da decisao. A margem de contribuhjao estabelecida, a priori, e de 16,6%. 
Uma estimativa de venda do produto do fornecedor nacional acaba resultando nas 
seguintes opera 9 oes matematicas: 


Quantidade 

Receita Lfquida de Vendas (RLV): 

Custo das Mercadorias Vendidas (CMV): 
Margem de contribuigao (MC): 

Margem de contribuigao (%): 


3.000 toneladas 

US$600,00/t*3.000/t = US$1.800.000,00 
US$500,00/t.* 3.000/t. = US$1.500.000,00 
US$1.800.000,00-US$1.500.000,00 = 
US$300,000,00 

US$300.000,00/US$1.800.000,00 = 16,6% 


TABELA 2.2 Probabilidades de sucesso e insucesso na 

entrega de mercadoria 

A 

Cenarios para a segunda etapa de decisao: escolha do 
fornecedor estrangeiro em fun 9 ao do risco conjuntural 

Probabilidade Probabilidade 

de sucesso de insucesso 

Cenario! (c,): Fornecedor frances 

90% 

10% 


Cenario 2 (c 2 ): Fornecedor chines 

40% 

60% 


Cenario 3 (c 3 ): Fornecedor ucraniano 

90% 

10% 


v._ 



J 
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TABELA 2.3 Possfveis resultados decorrentes das alternativas em fungao 
do mercado 

Decisao 

Lucro/Prejufzo Cp 
Aquecimento do mercado 

Lucro/Prejufzo C 2 : Mercado estavel 

Importar 

Margem de contribui^ao 
16,6% (US$300,000,00) 

Retorno do investimento em quatro meses; 
Perda de aplica^ao no mercado financeiro 
do capital investido em estoque importado 
por tres meses (ate uso completo 
do estoque importado): 1,2% a.m. em 
tres meses = (1 + 0,012) 3 = 3,6% de 
perda de oportunidade de investimento. 

Nao importar 

Taxa de retorno financeiro 
de aplicagao do capital por 
um mes, 1,2% a.m.; Perda de 
oportunidade de investimento 
de ganho de 16,6%. 

Taxa de retorno financeiro de aplicagao 
do capital por um mes, 1,2% a.m. 

Probabilidade 
de ocorrencia 
do cenario 

80% 

20% 

"v- 


_^ 


Por outro lado, caso a empresa opte por nao importar, ela podera aplicar seu capital 
em tftulos publicos e ter uma remuneragao media definida, a priori, de 1,2% a.m. O 
resumo das alternativas esta evidenciado na Tabela 2.3. 

Passo 4: Listar o retorno de cada combina^ao de alternativas 
e consequencias (Tabela Payoff) 

O resultado de cada decisao e um valor, payoff, como e conhecido na Teoria da De¬ 
cisao, que e o valor presente dos lucros dos diversos perfodos futuros envolvidos. A 
tabela que evidencia os payoffs foi desenvolvida pautando-se nas seguintes operagoes 
matematicas: 

Margem de contribuigao na decisao de importar: 

1. Importar em cenario de aquecimento do mercado 

US$1.800.000,00 - US$ 1.500.000,00 = US$300.000,00. 

2. Importar em cenario de mercado estavel 

Margem de contribuigao: US$300,000,00. 

Perda de aplicagao no mercado financeiro: US$1.500.000,00*0+ 0,012) 3 = 

US$54,000,00 

Margem de contribuigao esperada: US$300.000,00-US$54.000,00 = 

US$246,000,00 
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O valor monetario esperado da decisao de importar e obtido por meio da multi- 
plicagao das probabilidades de ocorrencia de cada cenario por valores de margem de 
contribuigao, tal como a seguir: 

US$300.000,00 x 0,80 + US$246.000,00 x 0,20 
= US$289,000,00. 

Margem de contribuigao na decisao de nao importar: 

1. Nao importar em cenario de aquecimento do mercado 

US$ 1.500.000,00 x 1,2% = US$ 18.000,00. 

2. Nao importar em cenario de mercado estavel 

US$ 1.500.000,00 x 1,2% = US$ 18.000,00. 

O valor monetario esperado da decisao de nao importar e obtido a partir da multi- 
plicagao das probabilidades de ocorrencia de cada cenario por valores de margem de 
contribuigao, tal como se demonstra a seguir: 

US$ 18.000,00 x 0,80 + US$ 18.000,00 x 0,20 = US$ 18.000,00. 

Na Tabela 2.4, ve-se um resumo das operacoes matematicas realizadas anteriormente, 
considerando as probabilidades associadas as decisoes a serem tomadas. 

Segundo o metodo de analise dos resultados da tabela de payoff, designado ante¬ 
riormente, deve-se escolher a decisao que apresente maior Valor Monetario Esperado 
(VME), ou seja, deve-se decidir pela importagao, pois o maior VME apresentado e 
de US$289,200,00. 

Passo 5: Selecionar um dos modelos utilizados naTeoria da Decisao, 
aplica-lo e tomar a decisao 

Ao analisar as variaveis e o ambiente do problema, foi escolhido o metodo de Arvore 
de Decisao devido ao fato de trabalhar com fatores de analise em um ambiente de 
risco, com atribuigao de probabilidade de ocorrencia de certas condigoes do mercado. 


TABELA 2.4 Valor monetario esperado de acordo com o mercado e as decisoes 

Decisao 

Lucro/Prejuizo 

Gp aquecimento mercado 

Lucro/Prejuizo 

C 2 : mercado estavel 

Valor Monetario 
Esperado - VME 

Importar 

US$300,000,00 

US$246,000,00 

US$289,200,00 

Nao importar 

US$18,000,00 

US$18,000,00 

US$18,000,00 

Probabilidade 

80% 

20% 

- 

V_ 



_ 
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m _ 20% 

^ 0 







L_J 18.000 

^ 18.000 






falso 

jBJARIO: probabilidade de aquecimento do mercado 



1 1 0 

^ 18.000 








f—| 80% 

^ ° 







L_J 18.000 

^ 18.000 





Escolha de 

h DECISAO: Deve-se Importar? 






Fornecedor 

T 

289.200 







\ 



r— 1 20% 

^ 0,2 





\ 



_J 246.000 

^ 246.000 






im VERDADEIRO 

ACNARIO: probabilidade de aquecimento do mercado 



0 

^ 289.200 








80% 

^ 0,8 







L_l 300.000 

^ 300.000 





FIGURA 2.2 Primeira etapa da construcao da arvore de decisao. 


Utilizando a tabela de payoff, tem-se, na Figura 2.2, o desenvolvimento da Arvore de 
Decisao baseando-se nos dados da primeira decisao. 

Observa-se, na Figura 2.2, que o ramo da decisao que leva o nome “VERDA- 
DEIRO” deve ser considerado como a decisao correta. Para finalizar a construcao da 
Arvore de Decisao, incorporando os cenarios de conjuntura do fornecedor estrangeiro, 
e necessario fazer a analise de atributos e custos envolvidos na avaliaqao de uma 
proposta de fornecimento, para explicitar a margem de contribuifao de cada fornecedor 
e, subsequentemente, atribuir a probabilidade de risco conjuntural. 

Passo 6: Levantamento dos atributos no processo de escolha 
do fornecedor 

Conforme os procedimentos descritos, identificou-se, primeiramente, o problema gerado 
pela situafao presente na empresa importadora. Para atender a demanda crescente do 
mercado, a empresa distribuidora de afo para construcao visa aumentar seus estoques. 

O fato de nao conseguir o produto junto aos grandes fabricantes nacionais, que o 
usam para consumo proprio, ou exportam a maior parte de sua producao, leva a em¬ 
presa a buscar alternativas de fornecimento. Segundo o levantamento de fabricantes 
de aqo mundiais, foram selecionados os seguintes pafses que possuem fabricas que 
possam atender a demanda do mercado brasileiro: Franca, China e Ucrania. Os criterios 
importantes na avaliaqao de um fornecedor sao citados e analisados como atributos 
do produto no Quadro 2.1 . 

A decisao a ser tomada consiste na escolha de um dos tres fornecedores fabricantes 
de aqo para construcao - frances, chines ou ucraniano -, sendo que a industria chinesa e 
conhecida por sua producao crescente e seus precos acessfveis; a industria ucraniana 
e conhecida mundialmente por sua qualidade e seus precos competitivos; e a industria 
francesa, preferida por sua qualidade e tradiqao. 

A simplificacao de nomear fornecedores de acordo com um pais que identifique 
so a sua origem deve-se ao fato de que os precos de varios fornecedores do mesmo 
pais tendem a um valor similar. Entretanto, a media de precos de um pais difere da 
media do outro. 
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QUADRO 2.1 Criterios de avaliagao do fornecedor 


Criterios 

Opgao nao 
desejavel 

Opgao 

intermediaria 

Opgao 

excelente 

Prego comparado com target 
price 9 , baseado no preqo do 
fornecedor nacional 

Acima do target 
price 

Igual ao target 
Price 

Abaixo do target 
price 

Qualidade de acabamento 
final do produto 

Ruim 

Boa 

Excelente 

Qualidade de acordo 
com exigencia do mercado 
brasileiro 

Nao corresponde 

Corresponde 

parcialmente 

Corresponde 

totalmente 

Prazo de entrega 

Proximo semestre 

Entrega no prazo 
de 30 dias 

Entrega imediata 

Forma de pagamento 

Pagamento 

antecipado 

Carta de credito 10 

Pagamento a 
prazo 

Embalagem do produto 

Ausencia da 
embalagem propria 
para transporte 

Embalagem 
normal 
de transporte 

Embalagem 

propria 

para exportagao 


Para o estudo de caso em questao, e mais viavel fazer a analise de tres pafses 
fornecedores do que expor os dados de todos os fornecedores ja analisados por pafs. 
O levantamento dos atributos e feito para avaliagao de cada fornecedor e pesara na 
hora de decisao sobre a melhor opgao. O Quadro 2.2 trata do enquadramento dos 
fornecedores de acordo com os criterios de avaliagao. 

Os tres fornecedores escolhidos apresentaram propostas que se encaixam nas 
opgdes intermediaria e excelente, pre-qualificando os fornecedores para o respectivo 
fornecimento. Os pregos apresentados por todos os fornecedores estrangeiros sao 
descritos de forma detalhada a seguir, estando abaixo daquele pago ao fornecedor 
nacional. Disso decorre que, descontando-se os custos envolvidos, sera possfvel 
obter-se uma margem de contribuigao maior do que com produto nacional, que e 
de aproximadamente 16,6%. A estimativa e analise dos custos e da margem de con¬ 
tribuigao na importagao estao apresentadas a seguir. 

Quantidade a ser importada: 3 mil toneladas 

Valor FOB * 11 da Importagao (FOB): Prego Unitario * Quantidade 

Custo da Carta de Credito (C/C): 3% do valor FOB 


9. Prego desejado pelo comprador. 

10. Documento de garantia bancaria de pagamento para o exportador desde que sejam cumpridas as condi- 
£oes estabelecidas pelo importador; uma forma de pagamento que apresenta seguranija para ambas as partes. 

11. Free On Board , modalidade de venda da mercadoria cujo preijo na fatura comercial do exportador 
incluira os custos de coloca 9 ao dentro do navio (transporte utilizado no estudo). Os custos que ocorrerem 
depois disso (fretes, transbordos, taxas, impostos e outros) sao por conta do importador. 
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QUADRO 2.2 Enquadramento dos fornecedores a partir dos criterios 
de avaliagao 


Criterios 

Opgao nao 
desejavel 

Opgao 

intermediaria 

Opgao excelente 

Prego 



Os tres 

Qualidade de acabamento final 
do produto 



Os tres 

Qualidade de acordo com 
exigencia do mercado brasileiro 



Os tres 

Prazo de entrega 


Chines / Ucraniano 

Frances 

Forma de pagamento 


Os tres 


Embalagem do produto 


Chines 

Frances / Ucraniano 


Custos da Importagao e outros custos (Cl): 47% do valor FOB 

Custo das Mercadorias Vendidas (CMV): (FOB + C/C + Cl) 

Receita Kquida de vendas no mercado nacional (RLV): US$600,00/t * 3.000t = 


US$1,800,000,00 

Margem de contribuigao na importagao (MCI): RLV - CMV 

Margem de contribuigao na importagao (MCI) %: MC/RLV 

Margem de contribuigao a priori (MC): US$300,000,00 (ver Tabela 2. 3) 

Diferenga entre MC a priori e MCI: MC - MCI 

Diferenga entre MC a priori e MCI, %: (MC - MCI)/RLV 


A Tabela 2.5 mostra o resumo dos calculos dos custos da importagao e da margem de 
contribuigao a ser auferida na operagao de importagao com cada um dos fornecedores. 

Para a conclusao da analise, faz-se necessaria a comparagao da margem de con¬ 
tribuigao do produto nacional com a do produto importado. A priori, sabe-se que a 
margem de contribuigao da venda do produto do fornecedor nacional e de 16,6%, o 
que resulta no valor de US$300,000,00 (ver Tabela 2.3). Ja na Tabela 2.6 evidencia-se 
o calculo da margem de contribuigao do produto importado, bem como sua comparagao 
com o produto nacional. 


TABELA 2.5 Calculo da margem de contribuigao na importagao 


Fornecedor Prego Unitario FOB C/C 


Cl 


CMV 


RLV 


MCI 


Frances 


310 


930.000 27.900 437.100 1.395.000 1.800.000 405.000 


Chines 


295 


885.000 26.550 415.950 1.327.500 1.800.000 472.500 


Ucraniano 


300 


900.000 27.000 423.000 1.350.000 1.800.000 450.000 
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TABELA 2.6 Comparativo da margem de contribuigao do produto nacional 
e do importado 

"■N 

Fornecedor RLV 

CMV MCI 

MCI (%) 

MC 

MC- MCI MC 

-MCI (%) 

Frances 1.800.000 

1.395.000 405.000 

22,5% 

300.000 

105.000 

5,9% 


Chines 1.800.000 

1.327.500 472.500 

26,3% 

300.000 

172.500 

9,7% 


Ucraniano 1.800.000 

1.350.000 450.000 

25,0% 

300.000 

150.000 

8,4% 


_ 








Observa-se, na Tabela 2.6, que a rentabilidade do produto importado e maior. Ao 
se trabalhar com o produto dos fornecedores frances, chines e ucraniano, obtem-se 
uma margem de contribuigao aproximadamente 6%, 10% e 8% maior do que a do 
produto nacional, respectivamente. 

A decisao em favor do fornecedor chines fica evidente, ja que a margem de con¬ 
tribuigao apresentada e a maior de todas. Porem, lembrando-se que se trata de uma 
decisao em condigoes de risco, deve-se avaliar a probabilidade de sucesso de cada 
fornecedor. Conforme visto na Tabela 2.2, as probabilidades de sucesso/insucesso para 
cada fornecedor, em fungao do risco conjuntural de intervengao polftica, foram de 90% 
a 10% para o fornecedor frances, de 40% a 60% para o fornecedor chines, e de 
90% a 10% para o fornecedor ucraniano. 

Havendo insucesso no fornecimento da mercadoria, o comprador incorre em 
custo de abertura da carta de credito. A Tabela 2.7 mostra o ultimo elemento dessa 
problematica, o modelo monetario, que transforma a decisao em valores monetarios 
esperados (VME). 

Nota-se que o maior VME de US$402,300,00 indica que se deve optar pelo fornece¬ 
dor ucraniano. A ultima etapa do processo decisorio e preencher a Arvore de Decisao, 
observando-se que o campo que representa o valor monetario de sucesso da decisao de 
importagao ja esta preenchido com a margem de contribuigao estabelecida, a priori , 
como valor igual a US$300,000,00 (ver valores na Tabela 2.4). 



TABELA 2.7 Payoff para decisao de escolha do fornecedor estrangeiro 


Decisao de 

escolha do 

fornecedor 

Probabilidade 

de sucesso 

Lucro/prejuizo 
em cenario de 

sucesso (US$) 

Probabilidade 

de insucesso 

Lucro/prejuizo 
em cenario de 

insucesso (US$) VME (US$) 

Frances 

90% 

405.000 

10% 

-27.900 

361.710 

Chines 

40% 

472.500 

60% 

-26.550 

173.070 

Ucraniano 

90% 

450.000 

10% 

-27.000 

402.300 
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O cuidado a ser tornado e nao incorrer em duplicidade de valores. Por isso, deve-se 
anular o valor de margem de contribuigao no cenario de mercado (US$ 300.000,00 e 
US$246,000,00) e langar o valor da margem de contribuigao no campo correspondente 
de cada fornecedor. Desta maneira, a probabilidade de ocorrencia de sucesso do 
fornecedor com melhor oferta sera a seguinte: 

• Probabilidade de aquecimento do mercado de 0,8 * Probabilidade de sucesso de 
0,9 = 0,72. 

• Probabilidades de sucesso dos outros fornecedores sao calculadas da mesma ma¬ 
neira e o mesmo procedimento vale para o ramo de combinagao decisao - cenario 
“Importar - Mercado estavel”. 

A margem de contribuigao nessa combinagao tera descontado a oportunidade de 
remuneragao do capital perdida (ver Tabela 2.3), conforme se encontra mensurado 
na Tabela 2.8. 

Nota-se que os valores mensurados na Tabela 2.8 vao gerar a melhor alternativa 
em decorrencia de o cenario de mercado permanecer estavel. Preenchendo os dados 
na arvore, tem-se o caminho para a melhor decisao a ser tomada, tal como mostra a 
Figura 2.3. 

A ferramenta apresenta os resultados obtidos na tabela payoff do fornecedor (ver 
a Tabela 2.7) e o ramo “VERDADEIRO” indica a melhor alternativa, com VME de 
US$402,300,00, do fornecedor ucraniano no cenario de aquecimento do mercado. 

A analise da Figura 2.3 apresenta os resultados para ambas as conjunturas, tanto 
para aquecimento do mercado como para permanencia do mercado estavel. Tem-se 
claros os resultados para alternativas de escolha de qualquer opgao analisada. O 
metodo utilizado apresenta resultados monetarios apos revelagao de informagao sobre 
o mercado: 

• Optando por importar e o mercado permanecendo estavel - seguindo o ramo 
“VERDADEIRO”, o maior valor monetario apresentado e de US$358,034. 

• Optando por importar e ocorrendo o aquecimento do mercado - o ramo “VERDA¬ 
DEIRO” que apresenta o maior valor monetario e de US$402,300. 


TABELA 2.8 Margem de contribuigao para decisao importar - 
estavel 

- cenario mercado 

Decisao de 

escolha do 

fornecedor 

Custo das 

mercadorias 

vendidas (US$) 

Oportunidade 

Margem de remuneragao do 
contribuigao capital portres meses, 
(US$) 1,2% a.m. (US$) 

Lucro/prejufzo do 
cenario de mercado 

estavel (US$) 

Frances 

1.395.000 

405.000 

50.825 

354.175 

Chines 

1.325.500 

472.500 

48.292 

424.208 

Ucraniano 

1.350.000 

450.000 

49.185 

400.815 

\_ 




_ 
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A possibilidade de alterar valores e, principalmente, as probabilidades, a medida que 
sao recebidas novas informafoes, proporciona flexibilidade gerencial ao tomador 
de decisao. Desta forma, o modelo proposto torna-se uma ferramenta agil e precisa. 



FIGURA 2.3 Eslruturacao completa da arvore de decisao para escolha de opcocs. 
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Vencida a principal etapa de modelagem do problema, considerada a mais com- 
plexa, o administrador podera recorrer ao modelo ja elaborado quantas vezes seja 
necessario. 

Vale salientar que, devido a complexidade do desenvolvimento do modelo baseado 
no estudo de caso, algumas abordagens foram simplificadas. Dentre as principais 
limitafoes, destacam-se: 

• o numero de dados quantitativos sobre as ofertas de fornecedores estrangeiros, 
reduzido as informacdes essenciais, com o intuito de facilitar o entendimento do 
problema no nfvel macro; 

• a avaliagao do grau de incerteza por meio da atribu^ao de probabilidade, reduzida a 
dois cenarios: de mercado e da conjuntura de polltica industrial no pals de origem. 
Isso foi feito com a ressalva de que arvores muito complexas sao de diffcil analise 
e visualizacao. 

CONSIDERAgOES FINAIS 

Pode-se afirmar que o uso de modelos de decisao no processo decisorio resulta em 
decisoes eficazes, mesmo que o resultado, a primeira vista, nao apresente a alternativa 
mais lucrativa. O modelo traduz a melhor opqao em valores monetarios esperados, 
considerando certas probabilidades de sucesso e insucesso. Portanto, a Teoria da 
Decisao deve ser utilizada nos processos decisorios em coiidicdes de risco para gerar 
respostas assertivas baseadas em metodos e evitar julgamentos puramente intuitivos 
dos tomadores de decisao. 

Um dos impedimentos para a difusao do metodo seria a falta de habilidade 
junto a ferramenta por parte do tomador de decisao, que e adquirida com acompa- 
nhamento do especialista e depois de algum tempo de uso. Os conceitos de no de 
decisao e de no de cenario devem ser bem entendidos e diferenciados na aplicaqao 
do problema. 

Em termos de julgamentos gerenciais no processo de decisao, os decisores nao 
devem deixar de considerar a possibilidade da presen£a de vieses ou pre-conceitos, 
sob pena de afetarem e comprometerem o processo como um todo. Os decisores 
envolvidos no processo de tomada de decisao devem ter em mente que, na busca pela 
melhor decisao, devem trabalhar com todas as informaqoes disponfveis, inclusive 
aquelas consideradas diffceis. 

Nao obstante os problemas de construqao da Arvore de Decisao, seria injusto nao 
mencionar as vantagens desse metodo, que consiste em um processo relativamente 
rapido de aprendizagem, pois a construqao de regras acontece em linguagem de 
usuario, alem de ser posslvel gerar regras em areas em que um profissional tern 
dificuldades de estruturar seus conhecimentos. E conveniente para a construqao 
de modelos nao parametricos e mesmo assim possui uma clara classificaqao e alto 
nfvel de predicao. 

Nesse contexto, a Arvore de Decisao e uma ferramenta de apoio no processo de 
tomada de decisao. Muitos softwares destinados a analise de dados ja dispoem de 
tecnicas de constru 5 ao de Arvores de Decisao. Em areas em que o prefo de um erro e 
muito alto, sao de grande ajuda no trabalho do analista ou do gerente. 
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Cabe destacar que as Arvores de Decisao garantem maior flexibilidade sobre os 
mecanismos tradicionais a medida que as decisoes sao definidas com maior grau de 
flexibilidade, visto que muitos eventos podem ser mapeados ao longo dos diversos 
ramos, envolvendo varias decisoes de sequenciamento. Alem disso, a Arvore de 
Decisao acaba incorporando conceitos de incerteza na analise de alternativas de inves- 
timentos ou financiamentos, o que a torna uma tecnica muito eficaz. 

Uma grande vantagem da aplica§ao da ferramenta quantitativa apresentada 
neste capitulo refere-se ao fato de que uma decisao de prioriza^ao de determinado 
investimento em relafao a outros nao e apresentada de forma intuitiva em uma 
“caixa preta”. Ela se encontra amparada em uma modelagem matematica robusta 
e transparente, que apresenta um encadeamento logico e racional de premissas e 
preferences por parte do decisor acerca de cada um dos atributos considerados 
durante o processo. 

EXERCICIOS 

1. Resolva, sem o auxilio de nenhum software, a arvore a seguir: 



2. Novembro de 2012. Um estudante precisava decidir o que fazer com seu Corsa 
2009. O gerente da oficina mecanica com que trabalha lhe telefonara pela manha 
dizendo que, ao fazerem a revisao de rotina no carro, verificaram que o motor 
nao estava em boas cond^oes e, para evitar um dano permanente, era preciso 
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um servigo imediato que custaria, provavelmente, R$ 1.200. Foi solicitada 
autorizagao para fazer o servigo, mas o estudante mostrou-se indeciso, pois 
sua situagao financeira era crftica, pelo menos ate quando voltasse a trabalhar, 
em novembro de 2013. Devido a este fato, ele resolveu minimizar suas des- 
pesas com o carro durante os 12 meses seguintes. Para poder tomar a melhor 
decisao, o estudante passou a manha investigando as alternativas posslveis. Por 
volta do meio-dia, concluiu que deveria optar por ou trocar o seu Corsa por um 
Palio novo, ou manter o Corsa. Se decidisse pela troca, evitaria o conserto. O 
dispendio de caixa liquido com o novo carro, dando o Corsa como troca, seria 
R$ 3.000, que ele precisaria pagar imediatamente. Se decidisse manter o Corsa, 
ainda nao estava convencido se deveria fazer o conserto de R$ 1.200. Mesmo 
que autorizasse o conserto, continuaria tendo que fazer reparos anuais de cerca 
de R$ 1.000 no carro. Se decidisse nao conserta-lo, ele achava que havia uma 
chance de apenas 40% de o carro quebrar e uma chance de 60% de ele nao que- 
brar. Se ocorresse um enguigo serio, parecia igualmente provavel que o conserto 
custaria R$ 2.000 ou R$ 3.000. Como ambos os carros teriam valores de mercado 
mais ou menos iguais em um ano, o estudante achava que este fato nao precisaria 
ser considerado. No entanto, ele achava que os custos com combustfvel eram 
relevantes. Estimou-os em R$ 2.400 para o Corsa e em R$ 1.100 para o Palio, 
no decorrer dos 12 meses seguintes. Analise o problema de decisao e recomende 
uma estrategia a ser seguida. 

3. Um explorador de petroleo precisa decidir se perfura ou nao determinado local. Ele 
nao esta certo se a perfuragao sera seca ou molhada. Seus lucres sao determinados 
pelas suas agoes e pela natureza do pogo: o custo de perfuragao e de R$ 700.000; 
se um pogo molhado for perfurado (50% de probabilidade), o lucre bruto e de 
R$ 1.700.000. 

a. O pogo deve ou nao ser perfurado? 

b. Por um custo de R$ 100.000, o explorador poderia fazer uma serie de sondagens 
sismicas capazes de determinar a estrutura geologica do local. As sondagens 
mostrarao se o terreno subterraneo tem: (i) nenhuma estrutura (ruim) - 40% de 
probabilidade; (ii) estrutura aberta (mais ou menos) - 35% de probabilidade; 
ou (iii) estrutura fechada (muito bom) - 25% de probabilidade. A probabilidade 
de o pogo estar seco depende do resultado do teste: (i) teste sismico revelando 
“nenhuma estrutura” - 75%; (ii) teste sismico revelando “estrutura aberta” - 
50%; ou (iii) teste sismico revelando “estrutura fechada” - 10%. O teste sismico 
deve ser feito? Qual e a agao otima agora? 

c. Se o teste sismico fosse perfeito, ou seja, indicando sempre “nenhuma estrutura” 
para um pogo seco e “estrutura fechada” para um pogo molhado, ate que valor 
estaria o explorador disposto a pagar por este teste? 

4. Um corretor de imoveis recebe uma proposta para colocar tres salas comerciais 
a venda. Se ele aceitar o negocio, tera que desembolsar R$ 700 para aquisigao 
de um terno e realizagao de um pequeno reparo no seu automovel. Alem disso, o 
custo para tentar vender cada imovel (ligagoes telefonicas, deslocamento, tempo 
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despendido etc.) pode ser estimado como equivalente a 0,40% do prego de venda 
do imovel. O corretor e obrigado a comegar as tentativas de venda pelo imovel A 
e nunca deve tentar vender dois imoveis simultaneamente. Se e quando ele vender 
o imovel A, podera escolher qual imovel tentara vender em seguida. Mas se, apos 
um mes, falhar na tentativa de vender o imovel A (ou qualquer outro), nao podera 
continuar tentando vender os outros imoveis. Mas ele tambem pode desistir de 
continuar tentando quando quiser; nao precisa esperar um fracasso. Sua comissao 
de venda e de 4% sobre o valor do imovel e os demais dados relevantes encon- 
tram-se a seguir. 


Imovel 

A 

B 

C 


Prego 

R$ 25.000,00 
R$ 50.000,00 
R$ 100.000,00 


Comissao 

R$ 1.000,00 
R$ 2.000,00 
R$ 4.000,00 


Custo Pr (sucesso) 

R$ 800,00 70% 

R$ 200,00 60% 

R$ 400,00 50% 


a. Ele deve aceitar o negocio? Se sim, em que ordem deve tentar vender os outros 
imoveis? Em que momento deve parar? 

b. O corretor esta inseguro quanto a probabilidade de vender o imovel A. Para 
quanto este parametro deveria mudar para mudar a decisao do corretor de 
aceitar ou nao o negocio? Quanto passaria a valer o negocio para o corretor 
para diferentes valores deste parametro? 

c. Faga uma analise de sensibilidade completa em relagao a probabilidade de venda 
dos tres imoveis. Considere que a probabilidade de vender A e, no minimo, 30% 
e, no maximo, 90%; a de vender B e, no minimo, 30% e, no maximo, 80%; e a 
de vender C e, no minimo, 25% e, no maximo, 75%. Qual dos tres parametros 
apresenta maior potencial para impactar o valor do negocio? 

d. Agora considere o perfil de aversao a risco do corretor: para ele, ter 50% de 
chance de ganhar R$ 1.600 e 50% de chance de ganhar R$ 5.600 nao equivale 
a ganhar, com certeza, R$ 3.600 (o valor esperado), mas, sim, R$ 2.500. Res- 
ponda novamente as questoes levantadas na letra a. 
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Objetivos de aprendizagem do capftulo 









Saber modelar urn problema real complexo em situagao de risco atraves da Simulagao 
de Monte Carlo. 

Aprender a gerar numeros aleatorios e a transforma-los em variaveis aleatorias 
com suas proprias distribuigoes de probabilidades. 

Entender como replicar a simulagao e coletar e analisar as estatfsticas de safda. 
Saber conduzir uma Analise de Sensibilidade acerca dos parametros mais importantes. 
Entender as vantagens de uma abordagem experimental em comparagao a uma 
abordagem analftica no tratamento de problemas complexos. 


INTRODU^AO 

A logfstica vem mostrando que seu uso adequado dentro das empresas pode trazer 
beneffcios aos processos, com consequente redugao nos custos. O processo mais 
custoso para as empresas e o de transporte. A atividade de transporte representa 
ate 60% dos custos logfsticos. 1 Nao so no Brasil, mas em todo o mundo, busca-se a 
melhoria na infraestrutura e nos processos dos transportes, na tentativa de redugao 
nos custos. 

Ja ouvimos ha muito tempo que o escoamento da produgao no Brasil nao funciona 
adequadamente, seja pelo estado de ma conservagao das estradas, seja por nao existir ou- 
tro modal, com a mesma disponibilidade do rodoviario, que possa escoar esta produgao. 

Como a malha ferroviaria brasileira esta concentrada nas regioes Sul e Sudeste, 
seria necessaria a utilizagao de outros modais que permitissem acessar as demais 


1. Gomes; Ribeiro (2004). 
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TABELA 3.1 Compara^ao das principals caracterfsticas dos modais 2 




Ferroviario 

Rodoviario 

Maritimo 

Dutoviario 

Aereo 


Velocidade 

3° 

2° 

4° 

5° 

1° 


Disponibilidade 

2° 

1° 

4° 

5° 

3° 


Dependencia 

3° 

2° 

4° 

1° 

5° 


Capacidade 

2° 

3° 

1° 

5° 

4° 


Frequencia 

4° 

2° 

5° 

1° 

3° 


Custo fixo 

4° 

1° 

2° 

5° 

3° 


Custo variavel 

3° 

4° 

2° 

1° 

5° 


2 Adaptagao de Bowersox; Closs (2001). 
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regioes. Estes modais poderiam ser o rodoviario e o hidroviario, por exemplo. Na 
verdade, por que o transporte de carga precisa ficar restrito a uma modalidade apenas 
de transporte se cada uma delas apresenta diferentes vantagens, tornando-se mais 
adequada, portanto, a trechos, circunstancias e situacdes especfficas? 

Para ilustrar essas diferentjas de adequa£ao entre os modais, e apresentada, a seguir, 
a Tabela 3.1, que compara as caracterfsticas mais relevantes dos modais de transportes, 
dando, assim, uma visao das diferen£as entre eles, possibilitando uma melhor decisao 
em rela£ao a que modal utilizar e em que momento. 

Partindo desse ponto, surge a pergunta: seria possfvel reduzir os custos com o uso 
da intermodalidade de transporte? 3 Ou, ainda, seria possfvel reduzir os custos logfs- 
ticos de transportes, para uma grande quantidade de carga e em um longo percurso, 
com o uso da intermodalidade de transporte, sem impactar de forma significativamente 
prejudicial a qualidade do nfvel de servico? 

Mas quais seriam as possibilidades adicionais de modais a serem utilizados? Dos 
cinco modais utilizados para o transporte de carga, deu-se enfase ao rodoviario, ao 
ferroviario e ao hidroviario. Os outros dois ficaram de fora, pois o aereo apresenta 
custo demasiadamente alto e o por dutos e por demais limitado a produtos bem es- 
pecfficos. 

Dessa forma, a sugestao consiste em dividir os percursos mais longos entre os 
modais rodoviario, ferroviario e hidroviario (nao necessariamente nessa ordem) e, 
assim, utilizar as vantagens de cada um para compensar as desvantagens apresentadas 
pelo outro. Ou seja, se nos detivermos ao que ha de melhor em cada modal, poderemos 
utilizar os pontos fortes de cada um, sobrepujando os pontos fracos de cada modal 
utilizado. 


3. Integra^ao dos servi^os de mais de um modo de transporte, desde o momento da carga inicial ate a des- 
carga final. 
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Para ajudar a responder as perguntas colocadas anteriormente, propomos a utili- 
za 9 ao da ferramenta de Simula 9 ao de Monte Carlo, uma vez que esta metodologia ja 
e utilizada em estudos logfsticos com bastante frequencia. 

SIMULA^AO DE MONTE CARLO 

A Simula 9 ao de Monte Carlo e um metodo de avalia 9 ao interativa de um modelo 
determinfstico que usa numeros randomizados como entradas. Esse metodo e mais 
utilizado quando o modelo e complexo, ou nao linear, ou quando envolve um numero 
razoavel de parametros de incerteza. Uma simula 9 ao pode envolver mais de 10 mil 
avalia 9 oes do modelo estudado, uma tarefa dificil que no passado so poderia ser 
realizada por supercomputadores. 4 

O metodo de Monte Carlo, um entre varios voltados para analise da propaga 9 ao 
da incerteza, tem como grande vantagem poder determinar como uma varia 9 ao ran- 
domizada, ja conhecida, ou o erro afetam o desempenho ou a viabilidade do sistema 
que esta sendo modelado. 5 

Atualmente, a complexidade dos processos logfsticos faz com que as ferramentas 
estaticas e determinfsticas tenham aplicabilidade pontual. A Simula 9 ao e uma ferra¬ 
menta que leva em considera 9 ao a dinamica e a aleatoriedade presentes nos processos 
logfsticos, avaliando-os com precisao muito maior e permitindo ganhos substanciais 
nas opera 9 oes. 6 

A Simula 9 ao e a tecnica mais utilizada no planejamento de sistemas logfsticos 
pelo fato de poder abordar questoes referentes a escolha dos modais de transporte e 
da localiza 9 ao de armazens, aos custos de transportes e de imobiliza 9 ao, ao nfvel de 
servi 90 aos clientes etc. Torna-se possfvel, assim, avaliar os diversos componentes do 
sistema e optar pelas formas mais adequadas de opera 9 ao. 7 

A Simula 9 ao tem acompanhado a evolu 9 ao da importancia da Logfstica no meio 
empresarial. Necessidades de aloca 9 ao de espa 90 , requisitos de movimenta 9 ao cada 
vez mais variados sao apenas algumas das questoes de distribu^ao e armazenagem que 
podem ser tratados. Uma vez que a Logfstica, em muitos aspectos, esta relacionada a 
movimenta 9 ao de opera 9 oes envolvendo o tempo, as mudan 9 as de locais e aos eventos 
probabilfsticos, a modelagem dinamica se torna fundamental. A principal vantagem da 
aplica 9 ao da Simula 9 ao em sistemas logfsticos consiste no fato de permitirem testar 
virtualmente as alternativas de opera 9 ao. 8 

Podemos entender, entao, a Simula 9 ao como uma das ferramentas mais apreciadas 
para os estudos de processos logfsticos. Ela nos permite representar os problemas de 
forma mais proxima a realidade, em fun 9 ao da quantidade de dados utilizados. E os 
resultados, entao, sao avaliados com maior precisao. de forma mais precisa. 



4. Shamblin; Stevens Jr. (1979). 

5. Moore; Weatherford (2005). 

6. Chwif (2008). 

7. Johnson; Wood (1996). 

8. Johnson; Wood (1996). 






C 44 ) PARTE | 1 Metodos de Apoio a Tomada de Decisao 


ELSEVIER 


A SIMULAgAO DO ESCOAMENTO DA PRODUgAO DE ARROZ 

Comparada a outras ferramentas, a Simulagao demonstra boa representagao das 
situagoes complexas e incertas do mundo real. Por isso, e considerada de grande 
importancia neste estudo para ajudar na verificagao, entre outros impactos, da posslvel 
redugao de custos logisticos a partir da utilizagao da intermodalidade de transporte. 

Assim, foi escolhida para este estudo a tecnica de Simulagao como o mecanismo 
metodologico, em fungao da sua capacidade para proporcionar a visao sistemica sobre 
a realidade operacional. A ideia consiste em estabelecer uma metodologia capaz de 
obter uma confirmagao referente a possibilidade de redugao dos custos logisticos 
atraves do uso da intermodalidade de transporte. Apesar de procurar avaliar tambem 
outros indicadores de desempenho, a enfase esta na analise dos custos de transportes, 
de modo a auxiliar a identificagao da melhor forma de agrupar os modais de transportes 
no escoamento da produgao, podendo assim alcangar a otimizagao na distribuigao dos 
produtos. 

A Simulagao foi aplicada a partir de parametros reais e ambientada em uma planilha 
do Excel, ilustrada parcialmente na Figura 3.1. Estamos considerando o transporte da 
produgao de arroz, partindo do Vale do Jacul (RS) ate a Regiao dos Lagos (RJ). 

Na Figura 3.1 sao apresentadas as variaveis aleatorias utilizadas para aplicagao da 
ferramenta de Simulagao no caso do transporte unimodal 9 rodoviario. 

A carga a ser transportada e de 1,200t, por um trecho de 1,922km, percurso a ser 
percorrido em dois dias, de acordo com a exigencia do cliente. A velocidade media do 
caminhao foi considerada como uma variavel aleatoria, com Distribuigao Triangular, de 
mlnimo de 40km/h, maximo de 70km/h e moda de 60km/h. Alem da velocidade media, 
foram consideradas como variaveis aleatorias (tambem com Distribuigao Triangular): 
as horas uteis por dia, sendo o mlnimo de 12 horas, maximo de 17 horas e moda (valor 
mais provavel) de 15 horas; e o custo de transporte por tonelada, sendo este modelado 
com um mlnimo de R$ 55, maximo de R$ 78 e moda de R$ 60. 

Foram simuladas mil ocorrencias para que o experimento pudesse ser validado es- 
tatisticamente. De cada uma delas sorteou-se a velocidade media e as horas uteis por 
dia, tendo sido calculada a partir delas a duragao do trajeto em horas e o consequente 
nlvel de servigo. Para efeitos deste estudo, foi considerado apenas o aspecto de nlvel 
de servigo referente ao percentual de vezes em que o prazo estipulado pelo cliente - 
dois dias - foi cumprido. Tambem foi sorteado o custo por tonelada e, a partir dele, 
calculado o custo total. 

Em seguida, foi abordado o caso da intermodalidade, em que os parametros tambem 
foram baseados em dados reais. Para este caso, tambem foi utilizada a Distribuigao 
Triangular para modelar o comportamento de todas as variaveis aleatorias. 

O primeiro trecho foi considerado como realizado pelo modal rodoviario, e a 
distancia do percurso que vai de Cachoeira do Sul ate Porto Alegre foi levantada, 
respeitando os outros parametros utilizados pelo modal em questao na primeira 
simulagao. O trecho realizado pelo modal ferroviario foi de Porto Alegre (RS) ate 
o Porto de Sao Francisco do Sul (SC). Por ultimo, para o modal hidroviario foi 


9. Utilizando apenas um modo de transporte. 






Prazo 

2 

dias 













Horas uteis por dia 













Mfnimo 

12 

horas 








Carga 

1.200 

ton 

Maximo 

17 

horas 








Distancia 

1.922 

km 

Moda 

15 

horas 
















Velocidade madia 


Custo/ton 









Mfnimo 

40 


Mfnimo 

R$ 55,00 









Maximo 

70 


Maximo 

R$ 78,00 









Moda 

60 


Moda 

R$ 60,00 

















Rodoviario 

Simulacao 

Aleatdrio 

Velocidade 

Duragao (h) 

Aleatdrio 

Horas dteis por dia 

Duragao (d) 

Chegou? 

Aleatdrio 

Custo/ton 

Custo Total 

988 

0,97 

66,82 

28,77 

0,44 

14,56 

1,98 

1 

0,48 

R$ 

62,71 

R$ 

75.251,31 

989 

0,99 

67,96 

28, 28 

0,31 

14,15 

2,00 

1 

0,74 

R$ 

61,08 

R$ 

73.297,22 

990 

0,29 

53,12 

36,19 

0,95 

16,33 

2,22 

0 

0,24 

R$ 

73,66 

R$ 

88.395,31 

991 

0,75 

61,31 

31,35 

0,18 

13,62 

2,30 

0 

0,21 

R$ 

59,50 

R$ 

71.398,74 

992 

0,99 

68,12 

28,21 

0,61 

15,04 

1,88 

1 

0,09 

R$ 

63,40 

R$ 

76.085,28 

993 

0,05 

45,25 

42,47 

0,15 

13,50 

3,15 

0 

0,21 

R$ 

59,15 

R$ 

70.983,56 

994 

0,27 

52,76 

36,43 

0,01 

12,30 

2,96 

0 

0,36 

R$ 

57,72 

R$ 

69.258,20 

995 

0,74 

61,19 

31,41 

0,35 

14,31 

2,20 

0 

0,47 

R$ 

61,65 

R$ 

73.981,70 

996 

0,88 

63,90 

30,08 

0,47 

14,65 

2,05 

0 

0,20 

R$ 

62,33 

R$ 

74.792,02 

997 

0,33 

53,97 

35,61 

0,44 

14,57 

2,44 

0 

0,04 

R$ 

62,12 

R$ 

74.538,70 

998 

0,66 

59,97 

32,05 

0,68 

15,20 

2,11 

0 

0,80 

R$ 

66,43 

R$ 

79.713,64 

999 

0,36 

54,67 

35,16 

0,43 

14,53 

2,42 

0 

0,47 

R$ 

62,59 

R$ 

75.105,07 

1.000 

0,26 

52,45 

36,64 

0,79 

15,56 

2,35 

0 

0,12 

R$ 

64,55 

R$ 

77.461,73 

Media 

0,50 

56,59 

34,39 

0,48 

14,59 

2,37 

10% 

0,51 

RS 

63,52 

R$ 

76.222,93 


FIGURA 3.1 Ilustra^ao parcial da planilha de simulacao (alternativa unimodal). 
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considerado o percurso do Porto de Sao Francisco do Sul (SC) ate o Porto do Forno 
em Arraial do Cabo (RJ). Ainda foram levados em considera 9 ao os custos de trans- 
bordo, do caminhao para o trem e do trem para o navio, considerando as variaveis 
mencionadas anteriormente. 

Torna-se relevante ressaltar que foi considerado o uso do container para o trans- 
porte da carga, uma vez que a sua utilizafao facilita o transbordo. 

Tambem no caso da intermodalidade, foram simuladas mil ocorrencias do ex- 
perimento. Para o modal rodoviario, utilizou-se o mesmo procedimento do caso 
unimodal. Para cada um dos outros modais, o ferroviario e o hidroviario, sorteou-se 
apenas a velocidade media e o custo por quilometro por cada mil toneladas, uma vez 
que consideramos uteis as 24 horas do dia. Em seguida, foi calculada a durafao da 
viagem. 

No caso do transbordo para ambos, houve o sorteio da duracao unitaria (em minutos 
por tonelada) e, a partir dela, o calculo da duragao total em dias. O custo de cada 
transbordo foi calculado a partir da respectiva dura 9 ao e do custo unitario (em reais 
por minuto). 

Por fim, foram somadas as duracoes (em dias) das viagens empreendidas pelos 
tres modais envolvidos na intermodalidade e dos transbordos. A partir deste total, foi 
gerado o nlvel de servigo intermodal (percentual de vezes em que o prazo estipulado 
pelo cliente - dois dias - foi cumprido). Tambem foram somados todos os custos 
envolvidos na simula^o intermodal, incluindo transportes e transbordos. 

Os principais indicadores de desempenho considerados na analise - duracao media 
da viagem (em dias), nlvel de serv^o e custo medio - foram coletados para o conjunto 
de mil simula^oes; os resultados estao apresentados na Tabela 3.2 a seguir. 

Os resultados das alternativas unimodal e intermodal estao apresentados lado a 
lado (para compara 9 ao), alem da diferen 9 a percentual ocasionada pelo uso da in¬ 
termodalidade em cada indicador de desempenho. 

Esses resultados, na verdade, sao dinamicos e mudam a cada vez que os valores 
da planilha sao atualizados. O conteudo da Tabela 3.2 e apenas uma fotografia, nao 
captando, naturalmente, o aspecto dinamico desses valores. 

Como pode ser observado, os valores obtidos para cada um dos modos de trans¬ 
portes - intermodalidade (rodoviario, ferroviario e hidroviario) e unimodalidade 
(rodoviario) - sao razoavelmente diferentes. No caso da dura 9 ao media em dias, o 
intermodal apresenta 2,12 dias, ao passo que o unimodal apresenta 2,36 dias. Isto 
surpreende, pois o intermodal se mostra com um tempo medio de entrega menor; 
enquanto esperava-se que o unimodal o superasse neste indicador. No caso do nlvel 
de serv^o, como consequencia da dura 9 ao media, a intermodalidade tambem supera 
a unimodalidade: temos 29,20% e 11,60%, respectivamente. Ja em rela 9 ao ao custo 
medio, era esperado um custo bem menor para a intermodalidade em rela 9 ao ao 
unimodal e, o que podemos observar, embora seja um custo menor e significante, a 
diferen 9 a nao e tao grande assim, sendo de R$ 63.910,78 para o intermodal e de R$ 
76.473,26 para o unimodal. 

Apos a simula 9 ao do cenario basico - carga de 1.200t e dois dias de prazo de 
entrega - foi feita uma analise de sensibilidade em rela 9 ao a esses dois parametros - 
variando a carga para 600t e 2.400t e o prazo de entrega para tres dias. 






TABELA 3.2 Indicadores de desempenho para cada possibilidade de transporte (intermodal e unimodal) - utilizando dois dias 
como prazo de entrega e 1.200t para a carga 




Intermodal 


Unimodal 

Diferenga (%) 





Rodoviario 

Transbordo 

Ferroviario 

Transbordo 

Fiidroviario 

Total 

Rodoviario 

Diferenga (%) 

Duragao media (d) 

0,16 

0,04 

0,92 

0,09 

0,91 

2,12 

2,36 

-10,18% 

Nfvel de servigo 






29,20% 

11,60% 

151,72% 

Custo medio 

R$ 10.258,49 

R$ 87,81 

R$ 28.507,94 

R$ 401,90 

R$ 24.654,63 

R$ 63.910,78 

R$ 76.473,26 

-16,43% 



4 ^ 
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TABELA 3.3 Indicadores de desempenho para cada possibilidade de transporte 
intermodal e unimodal - analise de sensibilidade para o prazo de entrega 
e a carga transportada 


Duragao media (d) Nfvel de servigo Custo medio 

1.1.200t-2 dias 

Intermodal 

2,12 

29,20% 

R$ 63.910,78 


Unimodal 

2,36 

11,60% 

R$ 76.473,26 

1.200t - 3 dias 

Intermodal 

2,13 

100,00% 

R$ 63.890,69 


Unimodal 

2,36 

96,40% 

R$ 76.662,44 

600t - 2 dias 

Intermodal 

2,08 

37,20% 

R$ 32.074,10 


Unimodal 

2,36 

11,80% 

R$ 38.270,63 

2.400t-2 dias 

Intermodal 

2,23 

7,30% 

R$ 127.873,30 


Unimodal 

2,36 

12,70% 

R$ 153.659,15 

- 




_ 


A simulagao foi rodada novamente para cada cenario, e os resultados foram co- 
letados, estando apresentados em conjunto com os do cenario basico na Tabela 3.3. 

E valido destacar o ultimo cenario, em que a altemativa intermodal mostrou-se mais 
rapida (2,23 dias) - na media! - do que a alternativa unimodal (2,36), mas encerrou 
um nlvel de servigo inferior (7,30% contra 12,70%). Tal aparente paradoxo deve-se 
a maior variabilidade das duracdes das viagens unimodais, que apresentam valores 
em media um pouco mais altos, mas com ocorrencias extremas, fazendo com que um 
percentual maior das viagens dure muito ou muito pouco, ficando abaixo da barreira 
dos dois dias mais vezes do que no caso intermodal. 

Tal tipo de constatagao so foi possfvel devido ao uso da Simulagao, visto que, nor- 
malmente, atraves de abordagens analiticas so se consegue observar o comportamento 
medio das grandezas relevantes. 

CONSIDERAgOES FINAIS 

Os resultados deste estudo foram obtidos a partir da Simulagao com dados reais e per- 
mitiram verificar a viabilidade da intermodalidade em uma situagao real de escoamento 
da produgao de arroz - de Cachoeira do Sul (RS) ate Arraial do Cabo (RJ). Este era o 
principal objetivo da pesquisa - a analise a respeito das vantagens e das desvantagens 
oriundas do uso da intermodalidade, neste caso especffico, alem de validar a utilizagao 
da Simulagao, instrumento que foi capaz de considerar as incertezas presentes nesta 
operagao logfstica e de subsidiar com informagoes uteis para o processo decisorio, 
como o impacto no custo e em outros indicadores de desempenho. 

Em relagao a incerteza, o exemplo real apresentado neste trabalho considerou es- 
tocastico o comportamento de diversas grandezas capazes de impactar o desempenho 
do sistema: velocidade media de transporte e de transbordo, horas dispomveis por 
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dia para o modal rodoviario e custos unitarios de transporte e de transbordo. Tais 
grandezas foram modeladas a partir da Distribuigao Triangular de probabilidades, 
mas a flexibilidade da Simulagao permite que qualquer distribuigao - principalmente 
nos casos em que nenhuma distribuigao teorica represente bem o comportamento da 
grandeza - seja usada. Tal caracterfstica da ferramenta permitiu que a realidade fosse 
modelada de forma bem acurada, uma vez que houve um esforgo na coleta dos dados 
necessarios para uma representagao realista do problema em questao. 

A partir desses resultados, foi possfvel observar que, no caso da duragao media, a 
alternativa intermodal apresenta uma quantidade em dias inferior a da unimodal. Es- 
perava-se que, com o uso da intermodalidade, o nfvel de servigo cafsse um pouco, sendo 
tal queda irrelevante diante da redugao dos custos. Mas em relagao a esse atributo, como 
consequencia da duragao media, a intermodalidade tambem supera a unimodalidade 
na maioria dos cenarios, o que causou certa surpresa. Tal fato pode ser explicado em 
parte: (i) pela presenga - no trecho em questao - de muitas montanhas no caminho dos 
caminhoes, diminuindo um pouco a velocidade media; (ii) e pelo desenho recortado 
do litoral brasileiro e da sua malha viaria costeira, que impossibilita os caminhoes de 
percorrerem uma linha reta (o que nao acontece com os navios), tornando a distancia 
unimodal maior do que o somatorio das distancias dos trechos da alternativa intermodal. 

Ja em relagao ao custo medio, esperava-se um custo menor para a intermodalidade em 
relagao ao unimodal, o que foi possfvel se observar em todos os cenarios. Cabe aqui uma 
ressalva, no sentido de que tais conclusoes sejam consideradas com precaugao. Esta pes- 
quisa visou analisar apenas questoes objetivas e quantitativas relacionadas a distribuigao, 
como custo e duragao da viagem. Mas e claro que a alternativa intermodal - apesar de ter 
se apresentado mais vantajosa nesses aspectos - envolve alguns inconvenientes potenciais 
- como aspectos humanos relacionados ao comodismo e aversao a complexidade, alem 
de outros tambem menos tangfveis, envolvidos nos transbordos e em outros momentos da 
operagao intermodal - que nao foram considerados e que podem constituir parte da causa 
da restrita utilizagao da intermodalidade no mundo real, ao menos no Brasil. 

EXERCICIOS 

1. Existe interesse no calculo da area sob a curva e x2 , entre os pontos 0 e 2. 
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1,000 
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Como nao existe uma fun£ao predefinida para o calculo da area de uma figura 
irregular como essa, precisarfamos obter a integral da fun 9 ao e calcular o seu valor 
entre os pontos 0 e 2. O problema e que tal funfao nao pode ser integrada; em 
outras palavras, nao e posslvel obter uma fun^ao analftica geral que represente a 
integral da fun 9 ao em questao. O que se sabe (e pode ser demonstrado) a respeito 
dessa fun 9 ao e que a sua integral definida - entre 0 e +oo - e igual a metade da raiz 
quadrada de pi, ou 0,8862, aproximadamente. No entanto, nao queremos saber a 
area desde 0 ate o infinito, mas apenas ate o ponto em que x = 2. Obviamente, este 
valor que buscamos e menor do que 0,8862, que pode ser considerado um limite 
superior para a area procurada. 

a. “Distor^a” a realidade e de a figura uma forma da qual sejamos capazes de 
calcular a area: o triangulo. Para tal, basta “fingirmos” que a curva e uma reta e, 
assim, teremos um triangulo formado por tal segmento de reta e os dois eixos. 
Calcule a area desse triangulo. 

b. Agora descubra a area procurada de maneira experimental. Simule o sorteio 
de 10 mil pontos dentro do retangulo imaginario (0;0) - (2; 1), que engloba a 
area procurada. Verifique que percentual (desses 10 mil pontos) caiu abaixo 
da curva (ou seja, dentro da area visada), para ter uma boa ideia da area que 
procuramos, ja que a area do retangulo e conhecida. 

c. Compare as duas abordagens. 

2. O proprietario de uma empresa que vende material de camping esta empenhado em 
melhorar o atendimento da demanda de barracas, produto que tern importante parti- 
cipa£ao na receita da empresa. A margem de contribu^ao bruta unitaria do produto 
e de R$ 25,00. As perdas de vendas por falta de estoques tem sido constantes, ja 
que, em geral, os clientes que procuram o produto querem leva-lo no mesmo dia e 
dificilmente aceitam esperar a chegada de novos pedidos. A capacidade do galpao 
(alugado) utilizado para estocagem das barracas e suficiente, mas R$ 2,00 sao 
cobrados por dia por barraca armazenada no imcio do perfodo (ou seja, voce deve 
considerar - para efeitos de calculo do custo de estoque - o estoque inicial em vez 
do estoque medio). Isso tem levado a empresa a ter certa cautela nas quantidades 
do produto incluidas nos pedidos realizados. A dificuldade do proprietario esta 
em equacionar os pedidos de barracas, de forma a atingir um nfvel de satisfafao 
que considera aceitavel - pelo menos 95% - e tentar maximizar o seu lucro. As 
demandas diarias do produto, em termos de probabilidades de ocorrencia, baseadas 
em dados historicos dos ultimos meses, sao: 

Demanda de barracas 11 12 13 14 15 16 17 18 19 

Probabilidade 7% 9% 14% 20% 17% 13% 9% 7% 4% 

Para cada pedido, e cobrado um frete de R$ 150. O prazo de entrega dos pedidos 
tambem e variavel e apresenta a seguinte distribuifao de probabilidades: 

Prazo de entrega dos pedidos (em dias) 2 3 4 

Probabilidade 25% 55% 20% 

No momento, o estoque e de 30 unidades e nao existem pedidos a serem recebi- 
dos (pendentes). Voce foi contratado pela empresa para elaborar um modelo de 
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administragao de estoques (ou seja, definir o ponto e a quantidade de ressupri- 
mento), objetivando maximizar o lucro e respeitando o nfvel mmimo de service) 
(atendimento da demanda) exigido. 

a. Simule a experiencia para 100 dias, alterando as variaveis referentes ao modelo 
de administragao de estoques (ponto e quantidade de ressuprimento) e buscando 
atingir os objetivos da empresa. Como se comportam os componentes do lucro 
e o atendimento da demanda em face das mudangas nas variaveis citadas? 

b. O que acontece com a busca da solugao otima, se o valor do frete passa a ser 
R$ 250,00? 

c. Refaga a simulagao mantendo o valor do frete original, mas mudando o valor do 
aluguel diario do espago para armazenar as barracas para R$ 3,00 por unidade. 

d. Elabore um relatorio final com as suas recomendagoes. 

3. Uma empresa comercial necessita dimensionar a capacidade de sua frota de entrega 
visando satisfazer certo nfvel de servigo. As entregas da empresa sao realizadas ape- 
nas em dias uteis (adota-se a premissa de que o prazo normal - que ocorre quando 
existe capacidade de entrega - e de um dia util) e as necessidades sao aleatorias, 
com distribuigao de probabilidades suposta uniforme, porem discreta, dependendo 
do dia da semana. A tabela a seguir descreve as necessidades de entrega minima e 
maxima, de acordo com o dia da semana. 



Segunda 

Terga 

Quarta 

Quinta 

Sexta 

Entrega minima 

21 

10 

11 

11 

11 

Entrega maxima 

45 

22 

23 

23 

23 


O objetivo e estudar a distribuigao do prazo de entrega dos produtos em fungao da 
capacidade diaria de entrega, que podera ser de 20, 22, 24, 26, 28 ou 30 volumes. 
Para isso, deseja-se estimar a media e o desvio-padrao desse prazo para cada uma 
das capacidades de entrega consideradas. Depois faga uma comparagao entre as 
possfveis capacidades de entrega. 
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Objetivos de aprendizagem do capi'tulo 

• Consolidar o aprendizado sobre modelagem e resolugao de problemas segundo as 
abordagens vistas nos tres primeiros capitulos do livro. 

• Aprender a abordar problemas em condigoes de incerteza de forma nao determinfs- 
tica, mas sem considerar (e precisar de) as probabilidades de ocorrencia dos eventos. 

• Saber atacar um problema de decisao por meio de diferentes abordagens: Pro- 
gramagao Linear, Criterios nao probabilfsticos, Arvore de Decisao e Simulagao de 
Monte Carlo. 

• Verificar que tipo de conclusoes cada uma dessas abordagens e capaz de fornecer. 

• Entender as diferengas, vantagens, desvantagens e aplicabilidade de cada uma dessas 
abordagens. 

• Visualizar um quadro comparativo entre essas abordagens de modo a saber identificar 
qual a abordagem mais adequada em vista das caracterfsticas do problema a ser 
atacado. 


INTRODUgAO 

O ambiente de tomada de decisao e dependente de quem decidira e das informagoes 
que se tem a respeito do problema. Podemos identificar tres ambientes no processo 
decisorio: o de certeza, em que o tomador de decisao (administradores de empresas 
ou gerentes) identifica com certeza os efeitos de todas as alternativas e decisoes; o 
de incerteza, em que o tomador de decisao nao e capaz de atribuir as probabilidades de 
ocorrencias aos varios resultados das alternativas identificadas; e, por fim, o ambiente 
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de risco, quando o tomador de decisao e capaz de atribuir probabilidades de ocorrencia 
aos cenarios. 1 

Apesar da enorme contribuigao das tecnicas quantitativas aplicadas na solugao de 
problemas gerenciais, em algumas situagoes elas tem potencial limitado. Uma decisao 
precisa ser tomada sempre que estamos diante de um problema que possui mais de 
uma alternativa para sua solugao. Num modelo de processo decisorio, e necessario 
analisar e identificar a situagao e o problema, desenvolver e comparar as alternativas, 
classificar os riscos de cada uma delas, escolher a melhor e executar (implantar) e 
avaliar (validar ou nao o processo utilizado), e para isso necessitamos de um modelo. 2 

Algumas situagoes se apresentam para a tomada de decisao, como: decisao em 
condigoes de certeza, que ocorre quando a decisao e assumida com pleno conhecimento 
de todos os estados da natureza; decisao em condigoes de risco, que ocorre quando sao 
conhecidas as probabilidades de cada estado, podendo variar entre 0% e 100% para 
cada estado; decisao em condigoes de incerteza, que ocorre quando nao se obteve total 
ou parcialmente os estados conhecidos e a probabilidade associada a elas; decisao em 
condigoes de competigao ou conflito, quando ha dois ou mais decisores envolvidos e 
o resultado depende da escolha de cada um dos decisores. Nestes tres ultimos casos, 
a decisao e considerada como nao deterministica. 3 

O aspecto nao determimstico pode resultar de uma fonte de imprecisao causada 
pela dificuldade de avaliar as agoes sob a influencia de um ou mais criterios, pela in- 
determinagao dos metodos de avaliagao dos resultados, pela classificagao ou opinioes 
amblguas ou ainda por eventos probabillsticos. O problema e que, quando existe a 
incerteza, frequente no mundo empresarial, algumas tecnicas nao sao adequadas. A 
solugao deterministica (sem considerar as incertezas) para um problema pode ser bas- 
tante desfavoravel para a situagao gerencial com a qual se depara a empresa. 

Existem diversas ferramentas quantitativas ja consagradas para adotar uma abor- 
dagem deterministica para resolugao de problemas, dentre elas a Programagao Linear 
(ou o seu caso particular, a Programagao Binaria). 4 Porem, em um contexto real das 
organizagoes, as incertezas existem, como, por exemplo, quanto a rentabilidade de 
um fundo de agoes, e neste caso nao deveria ser utilizado este tipo de abordagem. 

Tomemos a situagao de uma empresa que tem atualmente cinco propostas para 
desenvolvimento de projetos, cada uma com um percentual de expectativa de se 
tornar um contrato, baseado na experiencia dos profissionais do setor comercial da 
empresa. 

A empresa trabalha hoje com 33 analistas e 87 programadores, sendo que 21 
analistas e 61 programadores estao alocados em projetos atuais e 12 analistas e 26 pro¬ 
gramadores estao terminando seus projetos. Cada proposta apresentada pela empresa 
para se tornar um contrato tem as caracterlsticas mostradas na Tabela 4.1 . 

Quais propostas a empresa deveria desenvolver para obter maior lucro total, visto 
que a disponibilidade de recursos (analistas e programadores) e razoavelmente escassa? 


1. Gomes; Gomes; Almeida (2009). 

2. Gomes; Gomes; Almeida (2009). 

3. Gomes; Gomes; Almeida (2009). 

4. Hillier (2006); Medeiros (2008). 
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TABELA 4.1 Caracterfsticas de cada proposta 


A, 

Proposta Perspectiva 

Recursos necessarios 

Lucro (em R$) Qtd Horas 


de se tornar 

contrato 

Analistas 

Programadores 



PI 

50% 

5 

8 

780.000 

17.350 

P2 

100% 

4 

6 

450.000 

10.080 

P3 

10% 

5 

10 

900.500 

20.010 

P4 

30% 

6 

12 

1.633.000 

36.288 

P5 

60% 

3 

5 

303.500 

6.720 

V 





V 


Se a empresa tivesse certeza de que todas as propostas langadas seriam aceitas, 
transformando-se em projetos, a Programagao Binaria 5 poderia ser utilizada e a solugao se 
tornaria trivial. Neste tipo de situagao - ausencia de incerteza - o uso de aplicagoes como 
a Programagao Binaria - que utiliza variaveis de decisao do tipo 0 ou 1 - pode desempe- 
nhar um papel particularmente importante na busca da solugao. Inicialmente, entao, o 
problema sera tratado como se fosse deterministico, atraves da Programagao Binaria. 

Para obtermos a solugao otima para o modelo de escolha da melhor proposta (ainda 
sem considerarmos a incerteza quanto a efetivagao dos contratos), utilizamos a opgao 
de restrigao binaria do Solver do Excel. No caso especifico do problema, definimos 
os tres elementos basicos para o modelo: 

• Variaveis de decisao: fazer ou nao cada proposta. 

• Variavel-objetivo: lucro total obtido, que deve ser maximizado. 

• Restrigoes: recursos humanos exigidos/dispont'veis (analistas e programadores). 

Utilizando a Programagao Binaria, ainda partindo da premissa de que todas as pro¬ 
postas vao se transformar em contratos, a solugao do problema e obtida a partir da 
aplicagao do Solver. Neste caso, nao ha nenhum risco envolvendo a transformagao 
das propostas em contrato. Pode-se observar que, quando nao ha nenhuma incerteza, 
podemos escolher os projetos 3 e 4 que, alem de tudo, sao os de maior lucro esperado. 

A inadequagao dessa abordagem nada tem a ver com o calculo para descobrirmos, 
dentre as cinco propostas, quais deverao ser trabalhadas. O que falta nesse estudo sao 
as incertezas. Em um contexto organizacional real, nao e possfvel saber se as propos¬ 
tas vao se transformar em contrato: apenas se pressupoe a probabilidade de cada um 
deles se concretizar. 

A solugao para esse problema deterministico (sem a presenga de riscos) pode ser 
bastante desfavoravel para a situagao gerencial com a qual se depara a empresa. Ex¬ 
plicate: na verdade, e pouco provavel (10%) que a proposta 3 se transforme em projeto. 
Se a empresa elaborar apenas as propostas 3 e 4 (que e o que sugere a solugao obtida 


5. Moore; Weatherford (2005); Hillier (2006); Leao et al. (2006). 
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atraves da abordagem primaria), e bem possfvel que o projeto 3 nao se concretize, 
ficando a empresa com recursos ociosos, que poderiam ter sido utilizados, por exemplo, 
no projeto 1, apenas um pouco menos rentavel, mas com probabilidade de efetivagao 
muito maior. Por outro lado, nao podem ser elaboradas propostas demais, porque entao 
passa a existir o risco de nao haver pessoal suficiente para atender a todos os projetos. 

Evoluindo, entao, para uma abordagem mais realista - ja que a metodologia sim¬ 
ples e trivial da Programagao B inaria (que utiliza um modelo determimstico) nao se 
aplica adequadamente nesse caso em que temos risco e incertezas - precisaremos 
utilizar um modelo que os considere, necessitando ser tratado atraves de uma tecnica 
nao determinlstica como: os criterios Maximax, Maximin e Minimax; 6 a Arvore de 
Decisao; 7 ou a Simulagao. 8 

Dessa forma, a proposta principal deste capftulo consiste em analisar os resultados 
obtidos pela aplicagao de diferentes abordagens que considerem incertezas e riscos ao 
problema de escolha das propostas, de forma a mostrar que tipo de conclusoes essas 
abordagens sao capazes de obter em casos normalmente tratados de forma determinfs- 
tica, a fim de se constituir um quadro comparativo. 

A seguir e apresentada uma breve revisao dos conceitos acerca de Programagao 
Binaria, Criterios nao probabilfsticos, Arvore de Decisao e Simulagao. 

PROGRAMAgAO BINARIA 

Uma area de aplicagao muito importante da Programagao Linear e aquela que envolve 
problemas com uma serie de decisoes sim ou nao, fazer ou nao fazer. Tais variaveis sao 
chamadas variaveis binarias e os problemas de Programagao Inteira contendo apenas 
variaveis binarias sao denominados Programagao Binaria (ou PI 0-1 ). 9 

As variaveis binarias, 0 ou 1, desempenham um papel particularmente importante, 
ja que possibilitam incorporar instrugoes condicionais a decisoes do tipo sim ou nao, de 
forma semelhante ao Excel. Estas decisoes sao, por vezes, chamadas de dicotomicas, 
em um modelo de otimizagao. 10 

Os gerentes frequentemente precisam enfrentar decisoes do tipo sim ou nao. 
Portanto, a Programagao Inteira Binaria (PIB) e amplamente usada para ajudar nessas 
decisoes. Elas possibilitam lidar com um problema com cuja formulagao natural e 
dificil de se trabalhar e reformula-lo como um problema de PI pura ou mista. * 11 

CRITERIOS NAO PROBABILISTICOS 

O criterio Maximax se baseia em uma visao otimista do problema. Escolhido determina- 
do modelo, supoe-se que ocorrera o melhor evento possfvel. Neste criterio, o problema 
decisorio pode levar-nos a solugao otima, em que uma alternativa sera escolhida pela 


6. Gomes; Gomes; Almeida (2009); Corrar; Theophilo (2004); Moore; Weatherford (2005). 

7. Corrar; Theophilo (2004); Moore; Weatherford (2005). 

8. Hertz (1980); Corrar; Theophilo (2004); Moore; Weatherford (2005); Mendes eta/. (2006). 

9. Hillier (2006). 

10. Moore; Weatherford (2005). 

11. Hillier (2006). 
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maximiza§ao da fun§ao-objetivo. Este metodo e indicado quando buscamos a altemativa 
com o melhor desempenho de maximizagao. O criterio Maximax carrega consigo 
uma visao de mundo muito otimista, ou seja, deve-se escolher o melhor resultado de 
cada alternativa, e dentre eles o melhor dos melhores resultados. 12 O metodo pode ser 
definido como criterio de decisao otimista, que tem por objetivo encontrar o melhor 
resultado possfvel, a alternativa de decisao que maximizara o maximo resultado. 13 

O criterio Maximax avalia cada decisao pelo que melhor pode acontecer se o 
gerente tomar essa decisao. Neste caso eles avaliam cada decisao pelo maximo retorno 
possfvel associado aquela decisao, porem nao se deve confundir a decisao com o es- 
tado de natureza que produz o resultado otimo. 14 

Ja o criterio Maximin se baseia em uma visao pessimista do problema. Supoe-se 
que, escolhido determinado modelo, ocorrera o pior evento possfvel. A alternativa sera 
escolhida como aquela que tem a melhor entre as piores op 9 oes de todas as alternativas. 
Deve-se determinar o lucro mfnimo para cada alternativa e, em seguida, escolher a 
alternativa com o maior lucro mfnimo. 

O criterio Maximin, tido como extremamente conservador ou pessimista ao tomar 
decisoes, tem por objetivo encontrar a decisao que maximizara o mfnimo resultado 
entre varios possfveis, ou seja, ele indicara a alternativa que produz o maior resultado 
entre os mfnimos encontrados. 15 

Ao usar esse criterio, o tomador de decisoes deve ficar alerta para o que significa 
maximo ou mfnimo, dependendo de como sao expressos os resultados da matriz de 
decisao. 16 

Finalmente, a estrutura do criterio Minimax e baseada inicialmente desenvolven- 
do-se uma tabela de perdas de oportunidade, a seguir encontrando-se a oportunidade 
maxima de perdas em cada alternativa e selecionando-se entre elas a alternativa que 
apresenta a menor perda. 17 

Moore e Weatherford (2005, p. 385) destacam que “arrependimento e sinonimo 
de ‘custo de oportunidade’ de nao tomar a melhor decisao para um dado estado de 
natureza”. Em uma tomada de decisao, o gerente gostaria de tomar uma decisao que 
minimizasse o arrependimento, mas nao sabe que estado de natureza acontecera. Se 
o gerente nao conhece a probabilidade, o que se sugere e usar o criterio Minimax 
conservador, selecionando a decisao que faz o melhor do pior caso, ou seja, a decisao 
que tem o menor arrependimento maximo. 

ARVORE DE DECISAO 

A Teoria da Decisao e um conjunto de tecnicas quantitativas que tem por objetivo 
ajudar o gerente a sistematizar o problema e soluciona-lo. Nao ha solugao de um 
problema sem criterio, e e neste sentido que a Teoria da Decisao baseia-se em 


12. Moreira (2007). 

13. Corral" Theophilo (2004). 

14. Moore; Weatherford (2005). 

15. Corrar; Theophilo (2004); Moore; Weatherford (2005). 

16. Gomes; Gomes; Almeida (2009); Moreira (2007). 

17. Corrar; Theophilo (2004). 
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criterios preestabelecidos, havendo sempre espa 90 para novos criterios e novas 
contributes. 18 

Para construir uma Arvore de Decisao e necessario desenhar um grafico de modo que 
as informacoes adicionais obtidas pelo tomador de decisao, no transcorrer do processo, 
possam ser inclufdas como novas variaveis. A constn^ao da Arvore de Decisao respeita 
a sequencia utilizada na leitura da Tabela de Decisao, em que parte das alternativas leva 
em conta os cenarios e culmina com o indicador monetario que definira a escolha do 
tomador de decisao, que e o resultado que retornara em decorrencia da decisao. 19 

Existem dois criterios usualmente utilizados no ambiente de risco: o Valor Mo¬ 
netario Esperado (VME), que maximiza os resultados esperados e que, para cada 
alternativa de decisao, sua apura 9 ao e possfvel pelo calculo da media dos valores de 
seus diversos resultados, ponderada pelas probabilidades de ocorrencia dos eventos; 
e a Perda de Oportunidades Esperada (POE), que minimiza as perdas e refere-se a 
diferen 9 a entre o retorno otimo e o retorno recebido. 20 

Para determinar o valor esperado para cada alternativa e escolher a alternativa que 
tiver o melhor valor esperado, deve-se atribuir uma probabilidade a cada acontecimento 
(mutuamente exclusivos), calcular os valores esperados de cada a 9 ao, multiplicando 
cada valor consequente pela correspondente probabilidade, somar esses produtos e 
escolher a a 9 ao cujo valor esperado seja o maximo. 21 

Na analise de uma decisao cujas posslveis consequencias sao valores monetarios 
relativamente proximos, o VME pode ser utilizado como criterio de decisao quando 
da escolha de uma dentre varias alternativas. No entanto, nos casos em que os valores 
monetarios das varias alternativas sao substancialmente diferentes, como numa decisao 
estrategica, o VME nao mais reflete necessariamente como um tomador de decisao avalia 
uma estrategia de a 9 ao, pois nao leva em conta a sua possfvel aversao a grandes perdas. 22 

Uma alternativa ao VME e a Utilidade, que consiste na atratividade dos possfveis 
resultados de decisao, ja que a maioria das pessoas e avessa ao risco. A analise decisoria 
lida com esse comportamento de aversao ao risco. A fun 9 ao Utilidade e definida como 
uma medida de satisfa 9 ao, diferente para cada tomador de decisao, e que substitui o 
VME no momento de considerar o valor equivalente de um no de decisao. 23 

Uma Arvore de Decisao deve ser complexa o suficiente para incorporar todos os 
elementos-chave e valores que sao importantes para os pacientes e, ao mesmo tempo, 
simples o suficiente para ser compreensfvel e operacional. 24 

SIMULAgAO 

A origem da Simula 9 ao e muito antiga. Os Chineses ja a utilizavam em seus jogos 
de guerra ha 5 mil anos e os povos prussios utilizaram esses jogos no final do 
seculo XVIII para auxilia-los no treinamento militar de suas tropas. Durante a 


18. Moreira (2007). 

19. Corrar; Theophilo (2004). 

20. Corrar; Theophilo (2004). 

21. Ferreira (2004). 

22. Vatter etal. (1978). 

23. Moore; Weatherford (2005). 

24. Campolina; Ciconelli (2006). 
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Segunda Guerra Mundial, o matematico John Von Neumann criou o conceito 
denominado Simulagao de Monte Carlo, que consistia na simulagao direta de 
problemas probabilisticos. Atualmente com o desenvolvimento dos recursos 
computacionais, esse metodo e usado em diversas areas, desde a simulagao 
de fenomenos ffsicos complexos ate menos complexos, como a simulagao do 
resultado de loterias. 25 

Na Simulagao, procura-se obter informagoes sobre os relacionamentos entre as 
variaveis que compoem o sistema no sentido de predizer seus futuros desempenhos 
sob novas condigoes. O objetivo e criar um ambiente no qual a informagao sobre agoes 
alternativas possfveis e conseguida atraves da experimentagao. 26 

A DETERMINAgAO DO MIX IDEAL DE PROJETOS ATRAVES 
DEABORDAGENS NAO DETERMINISTICAS 

Alem das variaveis utilizadas na abordagem do modelo binario (decisao de fazer ou 
nao cada proposta, lucro esperado e restrigoes de pessoal), inclufmos, nas analises 
probabilisticas (ver Arvore da Decisao e Simulagao, a seguir), a variavel expectativa 
de transformar as propostas em contratos, estimada com a ajuda da experiencia dos 
gerentes, ja apresentadas na Tabela 4.1. 

Criterios Nao Probabilisticos 

O problema de escolha das propostas foi resolvido segundo os tres criterios: Ma¬ 
ximax, Maximin e Minimax. Para resolver esse problema, foram elencadas as 31 
solugoes possfveis criadas pela combinagao das cinco propostas, e os possfveis 
resultados para cada uma delas em fungao das propostas se transformarem ou nao 
em contrato. 

De acordo com o criterio Maximax, foi escolhido o melhor resultado possfvel; 
com o criterio Maximin, o resultado que produz o maior resultado entre os mfnimos 
encontrados; e com o criterio Minimax, a alternativa que minimizara a maior perda 
de oportunidade. 

Vale ressaltar que, segundo esses tres criterios nao probabilisticos, conforme o 
nome sugere, as probabilidades nao estao contempladas; levou-se em conta apenas 
que o projeto pode ou nao ser feito. Alem disso, foi considerada uma multa para as 
situagoes em que a empresa nao dispuser de recursos humanos suficientes para lidar 
com todas as propostas aprovadas. 

Criterio Maximax 

As variaveis de decisao dependem exclusivamente do tomador de decisao. Ele pode 
decidir fazer somente a proposta 1, ou so a 2 ou ainda optar por qualquer uma das 
31 possibilidades de solugao. Supondo que tenha escolhido fazer a 3 e a 4, o que vai 
acontecer com o lucro se somente a 1 viesse a ser aprovada? Nesse caso, a proposta 



25. Corrar; Theophilo (2004). 

26. Moore; Weatherford (2005). 
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1 nao foi feita e nao pode ser aprovada, e as outras duas foram feitas, mas nao apro- 
vadas, gerando um lucro nulo. Se a proposta 3 for aprovada, a empresa vai ganhar o 
lucro referente a ela. A melhor possibilidade e fazer as cinco propostas (1, 2, 3, 4 e 5, 
ou seja, C12345 x 12345), o que demanda 23 analistas e 41 programadores, e portanto 
nao ha recursos suficientes para fazer todas as propostas, pois o total de recursos dis- 
ponfveis e de 23 programadores e 12 analistas. 

Quanto seria perdido pelo desgaste comercial de assumir fazer as cinco propostas 
e nao ter cond^ao de executa-las? A empresa vive um momento instavel, com muitas 
ocorrencias, e por este motivo, para resolver este problema, foi criada uma variavel que 
definimos como “% residual de lucro no caso de falta de recursos”, que arbitramos como 
50% (ou seja, e como se, em cima do lucro, fosse aplicada uma “Multa” de 100% - 
50% = 50%) para o caso de a proposta tornar-se contrato e nao haver recursos suficientes 
para executar o projeto. Foi criada uma formula aplicada em todas as celulas para definir 
de quanto sera o lucro. Assim, tomemos como exemplo a rela£ao “Fazer as propostas” 
1 e 3 x “Propostas que vao virar contrato” 2345. Se a proposta 1 foi feita e nao foi 
aprovada, na celula 2345 e somado apenas o lucro da proposta 3. Aplicando o valor de 
Multa de 50%, o lucro maximo dos maximos foi de R$ 2.533.500,00, para as alternativas 
C34, Cl34, C234, C345, C1234, C1345, C2345 ou C12345. Para efeitos comparativos, 
processamos, respectivamente, os resultados para o caso de a Multa ser 70% (cenario 
pessimista) e 30% (cenario otimista). 

Os resultados encontrados com o valor da Multa de 70% (cenario pessimista) sao 
os mesmos encontrados para valor da Multa de 50%, porem, os resultados apresentados 
para o valor da Multa de 30% (cenario otimista) tiveram uma mudan 9 a significativa, 
mostrando que apenas a alternativa C12345 tera o valor maximo dos maximos, ou 
seja, R$ 2.846.900,00. 

Criterio Maximin 

Os parametros sao os mesmos considerados no criterio Maximax. Novamente 
utilizamos a Multa de 50%, e neste caso o valor maximo dos mfnimos apontou 
para fazer, indistintamente, qualquer conjunto de propostas, todos com lucro 
de R$ 0,00. Em razao da natureza pessimista desta abordagem, nao faz sentido 
elaborar a analise de sensibilidade no que se refere a variavel Multa, ja que, se 
tudo der errado, nenhuma proposta virara contrato e o valor da Multa nao estara 
influenciando a decisao. 

Criterio Minimax 

Os parametros e os resultados intermediaries sao os mesmos encontrados para o criterio 
Maximax. Novamente utilizamos a Multa de 50%, e neste caso o valor maximo dos 
mmirnos apontou para fazer as propostas 1, 3 e 4 com lucro de R$ 876.750,00. 

Para efeito de compara 5 ao, tambem foram feitas analises para a Multa com 70% e 
30%, respectivamente, e neles os resultados diferem do resultado apresentado para a 
Multa normal de 50%, sendo, respectivamente, R$ 1.533.500,00 para fazer as propostas 
3 e 4, e R$ 445.050,00 para as propostas PI, P2, P3 e P4. 

Para facilitar o entendimento, o Quadro 4.1 mostra o resumo das melhores situates 
encontradas por cada criterio/cenario. 
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QUADRO 4.1 Quadro comparative Maximax x Maximin x Minimax 
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Arvore de Decisao 

Foi feita uma Arvore de Decisao contemplando as 31 possfveis solucoes (as mesmas 
da sefao anterior), porem, nesse caso, foram utilizadas as probabilidades (expectativas) 
para considerar a chance de uma proposta virar ou nao contrato. Em cada uma das 
solu 9 oes foram consideradas as probabilidades de as propostas se transformarem em 
contrato para avaliar o resultado financeiro. A mesma multa mencionada na se 9 ao 
anterior foi considerada para esta metodologia. A Arvore de Decisao foi desenvolvida 
inicialmente segundo o criterio do Valor Esperado. Fazer as propostas e uma decisao 
do tomador de decisao. A arvore foi desenvolvida de modo a comportar a decisao de 
fazer uma proposta, duas propostas, tres propostas, quatro propostas ou ainda de fazer 
as cinco propostas. Para cada posslvel decisao, verificamos quais as probabilidades 
de as propostas virarem contrato. Uma das possibilidades e nenhuma virar contrato. 
Isto nao significa que nao foi feita nenhuma proposta, e sim que, nesse caso, nenhuma 
vai virar contrato. 

Se todas as propostas fossem feitas e todas virassem contrato, a empresa ganharia 
R$ 4.067.000,00, porem, nao haveria recursos humanos suficientes para o desenvolvi- 
mento de todos os projetos, causando atraso e desgaste comercial, que foi considerado 
a partir da aplica 9 ao de uma “Multa” que reduz (em um determinado percentual) o 
lucro obtido. Este percentual foi fixado inicialmente em 50% e a Multa foi aplicada 
somente nos casos especlficos em que nao houve recursos. 

Assim, a Arvore de Decisao, utilizando o software TreePlan, escolhe a solu 9 ao 
com melhor Valor Esperado. O no de decisao (fazer as cinco propostas) foi escolhido 
como o melhor resultado, com Valor Esperado de R$ 1.283.349,00. 

Segundo a execu 9 ao do software TreePlan, dado que haja 50% de chance de a em¬ 
presa ganhar R$ 780.000,00, e 50% de chance de ela ganhar R$ 0,00, o Valor Esperado 
calculado e de R$ 390.000,00, mostrado na Figura 4.1. Para o calculo da arvore atraves 
da Fmujao de Utilidade, foi coletado com o tomador de decisao o seu Nivel de Aversao 
ao Risco, tendo sido indagado qual seria a equivalente-certeza para esta situa 9 ao de 
risco descrita anteriormente. O valor de R$ 250.000,00 apareceu como resposta e foi 
aplicado na fun 9 ao Atingir Meta de forma a se obter o valor de Tolerancia ao Risco 
de R$ 494.827,00, conforme ilustra a Figura 4.2. 

A arvore tambem nos mostrou algumas informa 9 oes importantes. Se forem fei¬ 
tas quatro propostas, o melhor seria fazer as propostas 1, 2, 4 e 5, com lucro de 



FIGURA 4.1 No de evento - Valor Esperado. 
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F1GURA 4.2 No de evento - Tolerancia ao Risco. 



R$ 930.898,00. Se for escolhido fazer tres propostas, a melhor opgao seriam a 1, a 2 e a 
5, com lucro de R$ 891.905,00. Caso se decidisse por duas propostas, a melhor opgao 
consistiria nas propostas 1 e 2, com lucro de R$ 700.000,00. Por ultimo, se a opgao e 
fazer apenas uma proposta, a melhor seria a mimero 2, com lucro de R$ 450.000,00. 

Com a Arvore de Decisao construlda, procedeu-se a analise de sensibilidade 
com as variaveis mais importantes. Primeiro foi trabalhada a Multa, variando de 0% 
a 100%, de 1% em 1%. O output e a solugao otima, em que conclulmos que de 0% a 
73% o valor da Multa nao importa para efeitos de determinagao da melhor solugao e 
podem ser feitas as cinco propostas. Porem ao aplicarmos o percentual de Multa de 
74% o melhor seria fazer tres propostas, sendo o conjunto 1, 2 e 5 o melhor deles, com 
lucro de R$ 891.509,00, valor mostrado na Figura 4.3, quando aplicamos a analise de 
sensibilidade para o valor final da arvore. 

Percebe-se na Figura 4.3 que, a medida que a Multa aumenta, o valor da solugao 
com cinco propostas vai diminuindo; mas o valor da Multa passa a ser irrelevante 
quando a solugao otima consiste de apenas tres propostas, pois nunca faltarao recursos 
neste caso. 












r^J PARTE | 1 Metodos de Apoio a Tomada de Decisao 


ELSEVIER 


A Arvore de Decisao tambem foi analisada com cenarios pessimista e otimista. Em 
um cenario mais pessimista, utilizamos as perspectivas de transforma 9 ao das propostas 
em contrato como sendo: 40%; 80%; 5%; 20%; e 50%. A partir desta situa 9 ao, para 
uma Multa de 50%, encontramos um lucro de R$ 750.164,00, sendo a melhor solu 9 ao 
fazer as cinco propostas. 

A Analise de Sensibilidade do impacto da Multa na decisao otima mostra que de 
0% a 90% o valor da Multa nao importa e podem ser feitas as cinco propostas, porem, 
ao aplicarmos o percentual de Multa de 91%, o melhor seria fazer tres propostas, sendo 
o conjunto 1, 2 e 5 o melhor deles, com lucro de R$ 639.198,00, quando aplicamos a 
Analise de Sensibilidade para o valor final da arvore. 

A Analise de Sensibilidade do impacto da Multa foi tratada tambem com o cenario 
otimista (probabilidades de efetiva 9 ao dos contratos iguais a 60%, 100%, 15%, 40% e 
70%). O Quadro 4.2 mostra o resumo das melhores solu 9 oes encontradas pelo metodo 
da Arvore de Decisao (abordagens do Valor Esperado e da Fun 9 ao de Utilidade) em 
diferentes cenarios. 

Simula^ao 

Nesta abordagem, tambem foi incluida a variavel expectativa de transformar as propos¬ 
tas em contratos. Neste caso, o lucro total dependeu efetivamente de uma combina 9 ao 


QUADRO 4.2 Quadro comparativo - Arvore de Decisao 
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5 
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R$ 1.283.349,00 
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Realista 


Cenario 1 

0 - 73% 

5 

1,2, 3, 4 e 5 

R$ 1.103.547,00 
a R$ 898.056,00 


Cenario 2 

74% ou mais 

3 

1,2, 5 

R$ 891.509,00 

Pessimista 


Cenario 1 

0 - 90% 

5 

1,2, 3, 4 e 5 

R$ 766.253,00 
a R$ 639.912,00 


Cenario 2 

91% ou mais 

3 

1,2, 5 

R$ 639.198,00 

Otimista 


Cenario 1 

0-71% 

5 

1,2, 3, 4 e 5 

R$ 1.268.477,00 
a R$ 964.963,00 


Cenario 2 

72% ou mais 

3 

1,2, 5 

R$ 961.258,00 
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especffica de um grande numero de variaveis diferentes. A simulagao foi necessaria 
para verificar se cada proposta elaborada iria se transformar efetivamente em projeto. 

Tal verificagao ocorreu a partir da geragao de numeros aleatorios que, em conjunto 
com as probabilidades de sucesso na efetivagao de cada proposta, foi capaz de simular 
se cada contrato se efetivou, para cada experimento realizado. 

O trabalho procurou simular o que aconteceria com o lucro total e com adequagao 
do pessoal exigido ao pessoal disponfvel em cada uma de todas as possfveis solugoes 
para o problema de decidir quais propostas deverao ser feitas. Sao 31 solugoes pos¬ 
sfveis: fazer todas as propostas (PI, P2, P3, P4, P5), fazer somente uma (PI; P2; P3; 

P4 ou P5), fazer duas, tres ou quatro propostas. Foram realizados 5 mil experimentos 
para cada uma das solugoes. 

Para cada proposta foi sorteado um numero aleatorio, em cada simulagao. Em 
fungao deste numero aleatorio e da probabilidade de efetivagao de cada proposta, foi 
verificado em cada uma das 5 mil simulagoes se a proposta iria virar projeto ou nao, 
o que esta representado pelos valores 1 ou 0 na planilha, respectivamente. 

Tratando apenas os projetos que vao ser efetivamente feitos, foram calculadas as 
quantidades de analistas e programadores necessarios para a sua execugao em fungao 
das necessidades basicas de cada projeto. Estes valores foram comparados aos recursos 
disponfveis, para verificar se cada conjunto de recursos necessarios e viavel. Se ambos 
os conjuntos foram viaveis, a viabilidade geral (celula “OK geral”) recebeu o valor 
1. A celula “% de OK” contabilizou em quantas das 5 mil simulagoes os recursos 
disponfveis puderam atender a demanda por eles. Assim, uma medida de risco pode 
ser construfda (risco = 100% - % OK) para contemplar qual a probabilidade de tal 
conjunto de propostas feitas gerar um desgaste comercial por conta de nao haver 
recursos humanos necessarios para atender a demanda de todos os contratos efetivados. 

O outro indicador de qualidade da solugao foi o lucro total, fungao exclusivamente 
dos projetos que foram efetivados e do lucro esperado de cada um. Ao final da corrida 
de 5 mil simulagoes, foi possfvel ter uma ideia do lucro medio (ou esperado) da solugao 
analisada. 

Aqui e muito importante deixar claro uma premissa que foi utilizada na modelagem: 
mesmo que os recursos disponfveis nao sejam suficientes para atender a demanda de 
todos os projetos, ainda assim o lucro esperado de cada um deles sera considerado no 
calculo do lucro total. Tal consideragao e razoavel, porque a expectativa e que todos os 
projetos aprovados venham a ser feitos em algum momento, a partir da liberagao gradual 
dos recursos oriundos dos projetos que forem chegando ao fim. A realizagao nao imediata 
do projeto, no entanto, pode gerar um desgaste comercial junto ao cliente, e portanto 
esse risco precisou ser monitorado atraves da medida de risco introduzida anteriormente. 

A planilha apresentada na Tabela 4.2 ilustra a simulagao de parte dos 5 mil ex¬ 
perimentos realizados para a solugao C 2345. A partir da simulagao de todas as 
solugoes exequfveis, foi possfvel fazer um quadro sumarizando o lucro e o nfvel de 
risco previstos para cada uma. 

Vale ressaltar que esses resultados sao, na realidade, dinamicos. Cada vez que a 
planilha e atualizada, os valores aleatorios sao recalculados e novos resultados sao 
obtidos. As ilustragoes apresentadas consistem apenas em “fotografias” de uma dessas 
atualizagoes. Felizmente, com 5 mil simulagoes, os resultados sao robustos, ou seja, 
nao apresentam grandes variagoes entre as atualizagoes. 
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A Figura 4.4 mostra os mesmos resultados, porem, abordando a relagao Lucro x 
Risco e indicando um ponto de corte para um nfvel maximo de risco aceitavel (tendo 
a area comercial da empresa sinalizado ser interessante nao haver mais do que 5% de 
probabilidade de haver desgaste comercial). 

A linha pontilhada na Figura 4.4 indica o que seria um patamar desejavel de risco 
maximo toleravel (2,5%). A partir da analise do quadro e das figuras, verificamos 
que o conjunto de propostas PI + P2 + P5 nao traz risco algum e fornece uma 
perspectiva de lucro de R$ 1.000.000, aproximadamente, sendo esta, a princlpio, 
uma boa opgao. 

Se apenas a analise minimizadora de risco for utilizada, os executivos poderao 
ir adiante, sem a preocupagao de estarem correndo risco. Neste caso, o conjunto 
PI + P2 + P5 ou o P2 + P3 + P5 poderiam ser escolhidos. Verificamos que as propostas 
P2 e P5 aparecem, nas duas composigoes, com 0% de risco. 

Para executar todas as propostas, terfamos um nfvel de risco muito grande, de 
aproximadamente 29%. Comportando um pequeno risco (aproximadamente 2,5%), 
duas solugoes aparecem: PI + P3 + P4 e P3 + P4 + P5 (as propostas P3 e P4 aparecem 


TABELA 4.2 Resultado (parcial) da Simula^ao para o caso C 2345 

~"' s S 


PI 

P2 

P3 

P4 

P5 

Totais 

Fazer ou nao 

0 

1 

1 

1 

1 

para 

Expectativa 

50% 

100% 

10% 

30% 

60% 

100% 

Lucro R$ 

780.000 

450.000 

900.500 

1.633.000 

303.500 
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TABELA 4.2 Resultado (parcial) da Simula^ao para o caso C 2345 (cont.) 


ALEATORIOS PROJETOS ANALISTAS PROGRAM. Lucro 


Simulagao 

1 

2 

3 

4 

5 

PI 

P2 

P3 

P4 

P5 

Utilizados 

Ok? 

Utilizados 

Ok? 

OK geral 

Total 

1 

0,74 

0,75 

0,75 

0,56 

0,11 

0 

1 

0 

0 

1 

7 

1 

11 

1 

1 

R$ 753.500,00 

2 

0,73 

0,96 

0,16 

0,09 

0,25 

0 

1 

0 

1 

1 

13 

0 

23 

1 

0 

R$ 2.386.500,00 

3 

0,11 

0,15 

0,01 

0,70 

0,67 

0 

1 

1 

0 

0 

9 

1 

16 

1 

1 

R$ 1.350.500,00 

4 

0,68 

0,65 

0,74 

0,52 

0,13 

0 

1 

0 

0 

1 

7 

1 

11 

1 

1 

R$ 753.500,00 

5 

0,44 

0,48 

0,73 

0,10 

0,77 

0 

1 

0 

1 

0 

10 

1 

18 

1 

1 

R$ 2.083.000,00 

6 

0,72 

0,53 

0,67 

0,50 

0,67 

0 

1 

0 

0 

0 

4 

1 

6 

1 

1 

R$ 450.000,00 

7 

0,29 

0,51 

0,14 

0,45 

0,50 

0 

1 

0 

0 

1 

7 

1 

11 

1 

1 

R$ 753.500,00 

8 

0,34 

0,75 

0,62 

0,80 

0,51 

0 

1 

0 

0 

1 

7 

1 

11 

1 

1 

R$ 753.500,00 

9 

0,44 

0,51 

0,70 

0,18 

0,48 

0 

1 

0 

1 

1 

13 

0 

23 

1 

0 

R$ 2.386.500,00 

10 

0,82 

0,78 

0,04 

0,31 

0,07 

0 

1 

1 

0 

1 

12 

1 

21 

1 

1 

R$ 1.654.000,00 
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FIGURA 4.4 Lucro e risco de cada solucao, com ponto de corte. 
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nas duas combinagoes). A primeira delas apresenta um lucro esperado de aproxima- 
damente R$ 1.000.000. 

Ainda dentro do risco delimitado pela area comercial (5%), qualquer combi nagao 
que leve em conta a execugao de tres propostas pode ser recomendada, exceto os 
conjuntos PI + P2 + P4 e P2 + P4 + P5. As solugoes PI + P2 + P3 e PI + P4 + P5 
estao fora da zona aceitavel, porem, em algumas simulagoes, elas sao consideradas. 

Dentre as solugoes aceitaveis com tres propostas, as combinagoes PI + P2 + P5 
e P2 + P3 + P5 sao as mais lucrativas (em torno de R$ 1.000.000). Porem, a solugao 
PI + P2 + P3 + P5 tambem aparece na zona aceitavel de risco (5%, aproximadamente), 
mas com um lucro esperado de cerca de R$ 1.100.000. 

Verificamos que, alem de ser de interesse da empresa apresentar o maior numero 
possivel de propostas (por questoes estrategicas), devemos fazer as propostas PI, P2, 

P3 e P5, dentro do mvel de risco aceitavel e obtendo a maior expectativa de lucro. 

Como as probabilidades de as propostas se transformarem em contrato foram 
baseadas na experiencia dos gerentes e sao, portanto, subjetivas, convem realizar uma 
analise de sensibilidade a respeito do seu impacto. 

Tratamos tambem essa abordagem em um cenario mais pessimista, em que as pers- 
pectivas de transformagao das propostas em contrato poderiam ser, respectivamente: 
40%; 80%; 5%; 20%; e 50%, e em um cenario mais otimista, no qual as perspectivas 
de transformagao das propostas em contrato poderiam ser, respectivamente: 60%; 
100%; 15%; 40%; e 70%. 

O Quadro 4.3 mostra uma comparagao das melhores solugoes encontradas pelo 
metodo de Simulagao para os cenarios realista, otimista e pessimista. 

Analise Comparativa das Abordagens 

Nao Probabilfsticas x Probabilfsticas 

Como vantagem das abordagens nao probabilfsticas (Maximax, Maximin, Minimax), 
podemos considerar que nelas nao ha dependencia da subjetividade (e de uma eventual 
imprecisao) das probabilidades utilizadas na Arvore de Decisao e na Simulagao. Se 
estas probabilidades estiverem erradas (elas sao subjetivas), isto pode comprometer 
um pouco estas duas abordagens. Para as nao probabilfsticas, essa consideragao nao 
entra na composigao. 

De qualquer forma, trata-se de uma pequena vantagem; porem, em contrapartida, as 
abordagens nao probabilfsticas apresentam uma desvantagem porque estaremos tratando 
todas as propostas da mesma maneira, ou seja, todas as propostas podem ou nao se 
transformar em contrato, quando, na realidade, algumas tern mais chance do que outras. 

Comparagao entre os tres criterios nao probabilfsticos 
Como desvantagem, a abordagem Maximax apresentou multiplas solugoes otimas 
em todos os cenarios. A abordagem Maximin simplesmente nao funcionou para este 
problema. 

A abordagem Minimax foi a melhor: apresentou apenas uma solugao otima para 
cada cenario e revelou certa sensibilidade em relagao a Multa, ou seja, para cada valor 
da Multa, ela disponibilizou uma solugao otima diferente. 






VI 
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QUADRO 4.3 Quadro comparative - Simulagao 


Realista 

Solugao 

C 123 

C 125 


C 145 


C 235 

C 245 


C 1235 


Lucro R$ 

719.007,00 

839.369,50 


783.392,00 


564.151,00 

837.167,50 


875.276,50 


Risco 

1,40% 

0,00% 


2,10% 


0,00% 

7,80% 


2,10% 

Otimista 

Solugao 


C 125 

C 134 



C 235 


C 345 



Lucro R$ 


1.125.421,00 

1.282.459,00 



785.221,50 


984.605,50 



Risco 


0,00% 

4,40% 



0,00% 


3,80% 


Pessimista 

Solugao 




C 145 

C 234 



C 1234 



Lucro R$ 




783.392,00 

780.942,44 



1.011.647,50 



Risco 




2,10% 

1,40% 



3,00% 
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Arvore de Decisao x Simulagao 

Podemos considerar como uma pequena desvantagem da Simulagao o fato de a planilha 
ficar muito “pesada”. Alem disso, ela nao apresentou diretamente a solugao otima, 
diferentemente da Arvore de Decisao; ela apenas mostra o desempenho de cada uma 
das possfveis solugoes. Mas como a quantidade de solugoes possfveis nao e muito 
grande (apenas 31), foi possfvel verificar o desempenho de todas elas na Simulagifo, 
o que nao e uma desvantagem tao marcante. 

A confecgao da Arvore de Decisao e muito trabalhosa, e isso pode ser considerado 
uma pequena desvantagem; porem, como vantagem, ela nao ficou muito “pesada”. Na 
Arvore de Decisao, foi possfvel contemplar o perfil de risco do tomador de decisao. 

Na Simulagao isso nao foi possfvel. 

Podemos colocar como desvantagem da Arvore de Decisao o fato de que assumi- 
mos uma Multa para penalizar o lucro nos casos em que a quantidade de recursos dis- 
ponfveis nao foi suficiente para atender todas as demandas de projeto. A desvantagem 
deve-se ao fato de esta multa ter sido arbitrada. 

Na Simulagao, a analise do problema baseou-se na verificagao da probabilidade 
da demanda por recursos superar a oferta. Por este motivo, o tratamento foi menos 
arbitrario quando comparado com o da Arvore de Decisao: nesta ultima, arbitramos a 
Multa; na Simulagao, mensuramos qual a probabilidade de isto acontecer. A Simulagao 
cuidou melhor esta situagao, porem na Arvore de Decisao foi possfvel fazer uma 
analise de sensibilidade em relagao a esse parametro arbitrado, o que diminui bastante 
a desvantagem. 

Outra vantagem para a Simulagao e o fato de a Arvore de Decisao precisar conside¬ 
rar uma unica variavel como sendo o objetivo (no caso, o lucro a ser maximizado). Na 
Simulagao, e possfvel verificar o desempenho de cada uma das 31 possfveis solugoes 
nos criterios que forem considerados necessarios - no caso, alem do lucro, o risco 
envolvido na solugao e a propria quantidade de propostas elaboradas. Estes ultimos 
foram usados como criterios de desempate, possibilitando uma analise mais rica e, 
assim, uma decisao mais embasada. 

Vantagens e Desvantagens das Metodologias 

No Quadro 4.4 apresentamos um resumo das vantagens e desvantagens detalhadas 
anteriormente. 

CONSIDERAgOES FINAIS 

Este capftulo apresentou um problema real de determinagao de mix de projetos de 
consultoria, em que nao se sabia, a priori, quais das propostas elaboradas se trans- 
formariam efetivamente em projetos a serem feitos. Tal presenga de incerteza e 
suficiente para inviabilizar a abordagem mais tradicional para este tipo de problema: 
a Programagao Binaria. 

Tratou-se tal problema de forma nao probabilfstica pelos criterios Maximax, 
Maximin e Minimax, e de forma probabilfstica a partir de uma Arvore de Decisoes 
e da analise de risco viabilizada pela Simulagao, que se revelou adequada para tratar 
tal tipo de problema, sobretudo pelo fato de permitir que riscos sejam incorporados 
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QUADRO 4.4 Quadro comparative das metodologias 


Nao probabilfsticas 


Vantagem 

Desvantagem 


Nao ha dependencia 
de subjetividade presente 
nas probabilidades 

Todas as propostas sao tratadas 
da mesma maneira 

Criterios 

Maximax 


Apresentou multiplas solugoes 
em todos os cenarios 

Maximin 


Nao funcionou para este 
problema 

Minimax 

Apresentou apenas uma 
solugao diaria para cada 
cenario 


Sensibilidade em relagao 
a multa 

Probabilfsticas 


Vantagem 

Desvantagem 


Cada proposta e tratada 
de forma diferente 

Ha dependencia 
de subjetividade 

Abordagens 

Arvore da Decisao 

Apresentou uma solugao 
otima 

Confecgao rnuito trabalhosa 


Contemplou o perfil do risco 
do tomador de decisao 

Foi arbitrada uma multa 
para penalizar o lucro 

Foi escolhida apenas 
uma variavel objetivo 

Simulagao 

Foi verificado o desempenho 
de cada uma das 31 solugoes 
por varios criterios 

Planilha muito "pesada" 


Foi verificada a probabilidade 
da demanda por recursos 
superar a oferta 

Nao apresentou diretamente 
uma solugao otima 

Nao contemplou o perfil de 
risco do tomador de decisao 
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de forma satisfatoria, garantindo maior representatividade para o modelo e mais 
confiabilidade para os resultados. 

Por causa da sua natureza, a Simulagao apresenta uma desvantagem: a solugao 
otima nao e revelada. Neste caso em particular, isso nao representou um grande pro¬ 
blema, pois a quantidade de solugoes possfveis era finita e razoavelmente pequena, e 
foi possfvel verificar o desempenho de cada uma delas. Isso revelou outra vantagem da 
ferramenta: a possibilidade de considerar, para a escolha da solugao, nao apenas uma 
variavel-objetivo (o lucro esperado, no caso), mas tambem outros aspectos importantes, 
como o risco de desgaste comercial e a quantidade de propostas elaboradas, referentes 
a cada uma das possfveis solugoes. 


EXERCICIOS 


1. Diariamente, um jornaleiro compra Q exemplares de um jornal, ao custo unitario 
de R$ 1,00. Cada jornal e vendido por R$ 3,00, sendo as sobras ao final do dia 
perdidas (valor residual nulo). A demanda diaria de jornais, que e independente e 
identicamente distribufda para todos os dias, e descrita pela seguinte distribuigao 
de probabilidades: 

Demanda diaria (D) Probabilidade (%) 


16 

19 

24 

22 

19 


60 

70 

80 

90 

100 


O objetivo do problema e determinar o valor otimo de Q (a quantidade diaria de 

jornais comprada pelo jornaleiro). 

a. Trate o problema de forma determinfstica. Para tal, considere a demanda diaria 
como um parametro fixo, arredondando o valor esperado da sua distribuigao. 

b. Agora trate o problema de forma nao determinfstica; mas ainda nao considere 
as probabilidades, apenas a incerteza de a demanda poder assumir qualquer 
dos valores apresentados na tabela. Resolva utilizando cada um dos criterios 
nao probabilfsticos. 

c. Inclua agora as probabilidades e resolva o problema atraves de uma Arvore de 
Decisao. 

d. Ainda abordando de forma probabilfstica, simule, para os valores razoaveis de 
Q (60, 70, 80, 90 e 100), a atividade do jornaleiro por um perfodo de 100 dias. 
Para cada dia simulado, liste a demanda observada, a quantidade vendida, a 
receita, o custo e o lucro diario. Ao final da corrida, calcule a media para cada 
uma destas variaveis e o desvio-padrao do lucro. Faga um grafico da media e 
do desvio-padrao do lucro diario em fungao de Q. 

e. Reflita sobre as vantagens e desvantagens de cada abordagem, assim como 
sobre a sua adequagao a este e a outros problemas. 
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Parte II 


Metodos de previsao 


Metodos de previsao podem ajudar as empresas a, entre outros objetivos: 

• planejar e estimar melhor os valores dos investimentos a serem realizados, 
os estoques a serem criados, a capacidade de atendimento necessaria para 
a prestagao de um servigo, o tamanho da produgao etc.; 

• verificar o efeito da entrada de concorrentes; 

• descobrir onde e necessario o esforgo de venda e se criar cronogramas 
diferenciados de promogoes e descontos. 

Como resultado do processo de previsao, estimativas superestimadas (previsao 
acima do real) acarretam o chamado custo de excesso, que pode se manifestar 
na forma de, por exemplo, custo fixo maior, custo desnecessario de estoque, 
obsolescencia, perecibilidade ou gasto maior de propaganda. 

Ja estimativas subestimadas (previsao abaixo do real) ocasionam o chamado 
custo de falta, que pode assumir a forma de, por exemplo, perda da margem de 
contribuigao de um produto ou servigo, ordens em espera na fila de produgao 
ou consequencias negativas de uma demanda nao atendida. 

Por isso, o objetivo consiste em realizar previsoes que sejam o mais proximo 
possfvel da realidade, ou seja, criar niodelos que apresentem menores erros de 
previsao possfveis. 

Existem metodos qualitativos e quantitativos de previsao. Os metodos quan- 
titativos, que sao o objeto desta Parte II, podem ser subdivididos em dois tipos: 
os causais e os de serie temporal. 

Os metodos causais sao recomendados quando a variavel a ser predita esta 
fortemente correlacionada com outras variaveis que estao sujeitas a grandes 
variagdes. Eles costumam apresentar resultados bastante satisfatorios, mas 
carregam, intrinsecamente, dois inconvenientes: 

1. e necessario o acesso ao historico nao so da variavel a ser prevista, como tambem 
das variaveis capazes de impactar a grandeza que esta se querendo estimar; 

2. e preciso saber prever o valor futuro destas variaveis potencialmente impac- 
tantes, chamadas variaveis explicativas, para que a grandeza em questao 
possa ser bem estimada. 


75 





f 76 ) PARTE | II Metodos de previsao 


ELSEVIER 


O representante classico dos metodos causais e a Regressao Multipla, abordada 
nos Capitulos 5 e 6. Mediante a Regressao Multipla, tenta-se prever, em um call 
center que recebe chamadas de clientes com reclamagoes ou duvidas acerca de 
sua conta telefonica, respectivamente, o volume e a duragao media de ligagoes 
recebidas diariamente, em fungao de diversos fatores, como o dia da semana, a 
ocorrencia de feriados e a data da chegada e do vencimento da conta. 

Ja os metodos de serie temporal 1 pressupoem que os fatores temporals que 
influenciaram padroes da atividade no passado continuarao a faze-lo, de forma 
similar, no futuro. Em vista disso, eles envolvem a projegao de valores futuros 
de uma variavel, com base, inteiramente, em observagoes do historico dessa 
variavel; em outras palavras, outras variaveis nao entram na analise e, natural- 
mente, nao precisam ter seus valores coletados nem previstos. 

Portanto, esses metodos acabam sendo mais simples e, o que e mais importan- 
te, podem ser aplicados em uma gama maior de situagoes, 2 3 quando comparados 
aos modelos causais. No entanto, algumas vezes, seus resultados nao costumam 
ser tao satisfatorios. 

Existem diversos metodos de series temporais, 1 sendo um dos mais completos 
o da Decomposigao Classica, que tern, como consequencia natural, a sua com- 
plexidade, que pode ser considerada relativamente alta, sobretudo no que diz 
respeito aos procedimentos operacionais. Em fungao disso, este livro oferece, 
como ultimo material desta Parte II, um tutorial multimidia (texto-planilha) sobre 
o metodo da Decomposigao Classica, que esta apresentado no Apendice. 


1. Uma serie temporal consiste em um conjunto de dados numericos obtidos durante perfodos 
regulares ao longo do tempo. 

2. Notadamente quando nao estao disponfveis informagoes acerca de outras variaveis. 

3. Media Movel, Amortecimento Exponencial Simples, Analise de Tendencia, Metodo de Holt, 
Modelo de Winter etc. 






-Capitulo 5 

Previsao da demanda de liga^oes 
em um call center por meio 
de um modelo de Regressao 
Multipla 


Marco Aurelio Carino Bouzada 

MADE/UNESA 


Objetivos de aprendizagem do capitulo 

• Saber construir um modelo causal, capaz de explicar o comportamento de uma 
grandeza real com base no comportamento de outras variaveis. 

• Entender as vantagens de um modelo causal em relagao a um modelo de serie 
temporal. 

• Aprender a quantificar variaveis qualitativas e incluf-las no modelo por meio de 
variaveis dummy. 

• Conseguir aferir a capacidade de explicagao do modelo. 

• Conseguir aferir a capacidade de explicagao de cada variavel e entender como cada 
uma impacta a variavel dependente. 

• Saber refinar o modelo, tornando-o mais parcimonioso com a retirada de variaveis 
pouco explicativas. 

• Entender como obter e interpretar a equagao capaz de prever a variavel dependente. 

• Aprender a mensurar o erro de previsao decorrente da utilizagao de modelo desen- 
volvido. 


INTRODUgAO 

Call centers consistent em centros operacionais, instalados para o uso da telecomuni- 
cagao e das tecnologias computacionais, para automatizar uma quantidade volumosa 
de diferentes atividades e servigos telefonicos, tanto recebidos quanto originados pelo 
centra. Os centros do tipo inbound, onde as chamadas sao originadas pelos clientes, 
sao caracterizados como um sistema constitufdo de varias pessoas (atendentes), que 
recebem ligagoes telefonicas de outras pessoas, normalmente clientes - ainda que 
potenciais, que desejam obter alguma informagao, contratar algum servigo, comprar 
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algum produto, responder a alguma pesquisa, atualizar dados, registrar ocorrencias 
ou fazer reclamafoes, entre outras demandas. 1 

Os gerentes e planejadores dos call centers exercem t'uncoes que, hoje em dia, 
exigem muito mais empenho e capacidade do que no passado. Com mais produtos e 
serviqos sendo criados, disponibilizados no mercado, vendidos e assistidos, em relafao 
a momentos anteriores, eles precisam entender o que esta havendo nos call centers 
para saber como as suas caracterfsticas influenciam os indicadores de desempenho. 
Certamente a utilizaqao de simples palpites, experiencias de tentativa e erro, intuigao 
ou softwares de funcionamento invisfvel para o usuario nao representam a melhor 
alternativa para os gerentes das centrais de empresas que querem ser bem-sucedidas. 2 

Os supervisores dos call centers, responsaveis pelo planejamento operacional, 
precisam dimensionar eficientemente as suas equipes de modo a proporcionar nfveis de 
serv^o satisfatorios a custos razoaveis. A administraqao apropriada das centrais requer 
a estima 5 ao de algumas caracterfsticas operacionais que levam ao dimensionamento re- 
comendado da equipe de atendimento; dentre essas caracterfsticas, destaca-se a demanda 
de liga 9 oes de clientes, questao a ser permanentemente abordada neste tipo de negocio. 3 

Considerando-se que o custo de pessoal representa aproximadamente 70% do 
total da industria, a materia-prima para o processo de gerenciamento eficiente (nfveis 
satisfatorios de servi 90 a custos razoaveis) de um call center consiste na previsao da 
demanda de liga 9 oes e da consequente carga de trabalho, que precisam, entao, ser 
estimadas da forma mais acurada possfvel. Essa necessidade confere a previsao de 
demanda de liga 9 oes o status de atividade crftica durante a busca pela eficiencia e 
pela efetividade do negocio. Uma previsao acurada dos parametros do sistema e um 
pre-requisito para um mvel de servi 90 consistente, representando um dos passos mais 
importantes do processo como um todo. 4 

A previsao do volume de chamadas e um importante aspecto das opera 9 oes em call 
centers, sendo necessaria para a previsao de carga de trabalho, aloca 9 ao de mao de obra 
e planejamento de capacidade. Em geral, existem algumas regularidades nos dados de 
chegada das chamadas, como a dependencia entre os dias, ou seja, a correla 9 ao entre 
os volumes de dias consecutivos e possfveis ciclos semanais, mensais e ate anuais nos 
volumes de liga 9 oes; e a dependencia intradiaria, ou seja, a correla 9 ao entre chegadas 
em diferentes perfodos de tempo - manha, tarde ou noite - dentro do mesmo dia. 5 

O processo de previsao consiste tanto em arte quanto em ciencia. Em arte, porque 
o futuro, afinal de contas, esta sendo previsto e a acuracia da previsao vai depender 
em parte do julgamento e da experiencia de quern estiver conduzindo o processo. Mas 
tambem em ciencia: um processo matematico, passo a passo, que utiliza a historia 
passada para prever eventos futuros. Assim, um entendimento de tecnicas estatlsticas 
especializadas e necessario para o processo ser benfeito. Os administradores que 


1. Grossman et al. (2001); Hawkins et al. (2001). 

2. Mehrotra; Profozich; Bapat (1997). 

3. Weinberg; Brown; Stroud (2006); Klungle; Maluchnik (1997). 

4. Hall; Anton (1998); Alam (2002); Gans; Koole; Mandelbaum (2003); Reynolds (2005); Bouzada (2006); 
Weinberg; Brown; Stroud (2006). 

5. Shen; Huang (2005). 
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possuem softwares de gerenciamento de forga de trabalho, que automatizam o processo 
de previsao, nao podem pensar que isso apenas e suficiente; entender estes calculos e 
tao critico quanto possuir o software: nao so para verificar a acuracia dos resultados 
como tambem, e talvez de forma mais importante, para explicar os numeros para a 
gerencia. Ou seja, mesmo tendo as ferramentas, vale a pena entender os fundamentos 
do processo de previsao. 6 

Para gerar previsoes acuradas dos volumes futuros de chamadas, os gerentes dos 
call centers precisam juntar informacoes de varias fontes, desde previsoes do mercado 
ate, eventualmente, a duragao dos intervalos comerciais veiculados na mfdia. No entan- 
to, parte importante das informagoes provem das estatlsticas do historico de chamadas 
que, usualmente, podem ser ajustadas para as tendencias atuais ou esperadas. 7 

Os modelos de series temporais tern sido a abordagem mais utilizada para o trata- 
mento da previsao de chamadas telefonicas. No entanto, outros modelos de previsao 
foram aplicados com sucesso em alguns ambientes de call center, como, por exemplo, 
centrais de vendas cuja previsao de chamadas se baseia nas remessas de produtos 
recebidas e nas propagandas enviadas aos potenciais clientes. 8 

A atividade de previsao e inerentemente diffcil no caso dos call centers em razao 
do tamanho relativamente pequeno dos intervalos de tempo utilizados na pratica. 
Tecnicas de previsao estabelecidas, como o Amortecimento Exponencial de Winters 
e a Regressao Linear Multipla, sao uteis para este tipo de negocio. Geralmente, usar 
uma abordagem de previsao apropriada ira reduzir pela metade os erros de previsao. 
No entanto, muitos call centers encontram dificuldades nesse processo por conta do 
conhecimento tecnico necessario para adequar estas tecnicas aos complexos padroes 
de chamadas e ao aspecto pouco organizado dos dados. 9 

Em geral, se a variavel a ser predita esta fortemente correlacionada com outras 
variaveis que estao sujeitas a grandes variagoes, urn modelo multivariado e necessario. 
O metodo da Regressao Multipla, em particular, e relativamente robusto em relagao a 
afastamentos dos pressupostos da normalidade e da homocedasticidade. Alem disso, a 
eventual ocorrencia de multicolinearidade entre as variaveis explicativas nao prejudica 
a capacidade preditiva do modelo, causando apenas problemas quanto a sua parcimonia 
e de compreensao do real efeito das variaveis independentes sobre o comportamento 
da variavel dependente. 10 

O PROCESSO DE PREVISAO DE DEMANDA NO CALL CENTER 
ESTUDADO SEM A REGRESSAO MULTIPLA 

A principal ferramenta que a empresa estudada utiliza para auxiliar as decisoes geren- 
ciais de planejamento operacional e um software totalmente integrado a infraestrutura 
de telefonia. O produto consiste em uma ferramenta de forga de trabalho utilizada 


6. Reynolds (2005); Bouzada (2006). 

7. Matan; Nourbakhsh (1998). 

8. Mehrotra (1997). 

9. Grossman et al. (2001). 

10. Silver (2000); Levine; Berenson; Stephan (2000); Corrar; Paulo; Dias Filho (2007). 
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para dimensionar e controlar o atendimento, que auxilia os gerentes - entre outras 
atividades - na previsao de demanda (por horario e dia da semana). 

De acordo com a gerencia da empresa, a principal preocupa9ao consiste em con- 
seguir melhorar as previsoes de demanda de Iigardes. A acuracia das previsoes nao e 
boa para alguns produtos. Para o produto que abrange clientes que buscam servhjos 
relacionados a telefonia fixa, como mudan9a de endereco de cobranga e duvidas ou 
problemas com a conta telefonica, a empresa dispoe de todas as informa9oes que 
poderiam fornecer uma melhor previsao, estando apenas carente de um metodo mais 
adequado. Essas informa9oes sao basicamente as quantidades de contas enviadas para 
os clientes em cada data, normalmente cinco a seis dias antes do vencimento, e as suas 
proprias datas de vencimento. A impressao da gerencia e que a demanda de liga9oes 
para o produto e influenciada pelos eventos de chegada da conta e pelo seu proprio 
vencimento, datas em que haveria maior acesso ao servi90 por parte dos clientes. O 
problema e que o sistema utilizado para previsao de demanda nao leva em conta essas 
informa9oes. 

A magnitude do erro absoluto percentual medio, segundo a otica diaria, incorrido 
durante o processo de previsao de demanda para o produto em questao era da ordem 
de 5 %. Na opiniao dos responsaveis pela previsao de demanda, esse indicador de 
erro encontra-se em um patamar baixo, sendo considerado mais do que aceitavel. No 
entanto, estes responsaveis entendem que existe um potencial de melhoria no processo, 
de forma que esses erros podem se tornar ainda menores. 

O primeiro passo do processo de previsao de demanda para o produto em questao 
consiste na gera9ao de uma previsao basica, realizada pelo software de gerenciamento 
de for9a de trabalho. O segundo passo consiste em uma crftica, por parte dos analistas 
da coordena9ao, a esses numeros gerados pelo software. O terceiro passo consiste em 
incorporar aos valores previstos o possfvel impacto causado por alguns eventos es- 
pecificos, como a segunda-feira de Carnaval ou uma ocorrencia diferente em algum 
sabado, por exemplo. A maneira de incorporar tais impactos depende de o evento ser 
novo ou desconhecido e do previo conhecimento das suas consequencias. 

Alguns dos eventos especiais podem ter seu impacto previsto de uma forma melhor, 
quando a equipe de analistas e subsidiada por algumas premissas recebidas de outros 
setores capazes de impactar a demanda de chamadas, como: previsao de crescimento 
da base de clientes; implanta9ao de novos pianos e serv^os na telefonia fixa, como o 
piano de minutos, por exemplo. 

Essas premissas tambem envolvem o comportamento dos ciclos das contas tele- 
fonicas dos clientes, uma informa9ao oriunda do setor de faturamento da empresa 
telefonica e que ajuda a previsao de demanda do produto em questao. A informa9ao 
consiste, basicamente, nas quantidades de contas enviadas para os clientes em cada 
data, normalmente cinco a seis dias antes do vencimento. Ela e util porque a demanda 
de liga9oes para o produto e influenciada, segundo a equipe de previsao, pelos eventos 
de chegada da conta no destino e seu vencimento, datas em que haveria maior acesso 
ao servi90 por parte dos clientes, para fins de esclarecimento e/ou reclama9ao. 

Na verdade, durante todo o perfodo de quatro a cinco dias em que a conta esta 
com o cliente (desde a chegada ate o vencimento), ocorre um acrescimo no volume 
de liga9oes, mas o maior impacto acontece no momento da chegada da conta na 
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residencia do cliente, na opiniao da coordenagao responsavel pelo produto. Segundo 
os analistas, o impacto pos-vencimento existe quando o cliente que esqueceu de pagar 
a conta liga para saber como proceder, mas e pequeno, nao sendo significativo em 
termos de aumento de demanda. 

O problema e que o software nao leva em conta isso; a sazonalidade do dia da 
semana e considerada, mas nao a distancia do dia em relagao a chegada da conta ou 
ao seu vencimento. Por exemplo, a conta de determinado cliente vence no dia 10 de 
cada mes, e nao na segunda quarta-feira; entao, no mes em que o dia 10 cafsse em uma 
terga-feira, esse dia da semana deveria ser mais impactado do que os outros, enquanto 
no mes em que o dia 10 caisse em uma quarta-feira, o impacto maior deveria ser 
inclufdo neste dia. Mas o software considera, em termos de sazonalidade, todas as 
tergas-feiras da mesma forma, seguindo o mesmo raciocfnio para os outros dias da 
semana. Alem disso, ha meses em que as contas nao sao enviadas exatamente cinco 
dias antes do vencimento, e a chegada da conta ao destino tambem deve ser considerada 
como um evento especial. 

Essas informagoes sao levadas em conta no processo de previsao, mas nao de 
forma automatica. Toda a analise e feita pela equipe de previsao, e nao pelo software. 
Mas durante esse processo, os analistas nao levam em conta um possfvel impacto da 
vespera e da antevespera do vencimento da conta e nem do perfodo pos-vencimento; 
apenas o efeito dos eventos chegada e vencimento. 

O PROCESSO DE PREVISAO DE DEMANDA COM A REGRESSAO 
MULTIPLA 

A ideia consiste em fugir dessa pos-analise nao metodica - que depende muito da 
intuigao do analista - e poder contar com uma ferramenta que use as informagoes do 
faturamento no processo de elaboragao das previsoes, para que se possa evitar o traba- 
lho manual, que aumenta a possibilidade de erro. A ferramenta nao tem a pretensao 
de substituir a intuigao e a experiencia da equipe de previsao, mas servir de apoio a 
tomada de decisao. Ou seja, a opiniao dos especialistas nao deve ser descartada; ela 
deve, sim, ser utilizada para enriquecer o modelo que busca, basicamente, objetivizar 
as informagoes subjetivas que sao relevantes para o processo de previsao. 

O metodo escolhido foi o da Regressao Multipla, pois e capaz de captar o impacto 
especffico de cada evento e informagao na demanda de cada dia. Por exemplo, qual o 
impacto, na demanda de ligagoes, do fato de o dia em questao cair em uma quinta-feira? 
Qual o impacto de o dia consistir em um feriado? E de chegarem 200 mil contas para os 
clientes no dia em questao? E de o dia ser a vespera do vencimento de 500 mil contas? 

A ideia e verificar a distancia do dia especffico em relagao aos eventos especiais, a 
chegada e o vencimento das contas, e calcular para quantas contas esse e o dia C (dia de 
chegada), o V (dia do vencimento), o V - 1 (vespera do vencimento), o V - 2 (anteves¬ 
pera do vencimento),..., V + 1 (dia seguinte ao vencimento), V + 2,... entre outros. 

De posse do cronograma de faturamento fornecido pela equipe de previsao, para o 
perfodo de julho de 2005 a julho de 2006, foi possfvel consolidar as informagoes em 
termos da data de chegada e da data de vencimento das contas, conforme apresentado 
parcialmente nas Tabelas 5.1 e 5.2, respectivamente. 
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TABELA 5.1 Contas agrupadas por data de chegada prevista (setembro a outu- 
bro de 2005) 

Data de Chegada Prevista 

Total 

06/09/05 

3.360.569 

13/09/05 

801.143 

14/09/05 

432.670 

15/09/05 

709.437 

26/09/05 

1.315.342 

27/09/05 

820.855 

28/09/05 

1.555.097 

29/09/05 

2.506.038 

04/10/05 

1.808.735 

05/10/05 

1.167.447 

07/10/05 

800.105 

11/10/05 

1.890.364 

18/10/05 

709.437 

26/10/05 

1.731.518 

28/10/05 

587.118 

Total 

20.195.875 

- 

_^ 


Para poder ser desenvolvido o modelo de regressao multipla para prever a quan- 
tidade de ligagoes em cada data, foi preciso ter acesso ao historico desta grandeza. 
Por motivo de sigilo estrategico, esses valores nao puderam ser revelados, em sua 
plena magnitude; no entanto, a eles foi aplicada uma constante de proporcionalidade, 
e os resultados, de julho de 2005 ate margo de 2006, fornecidos pela empresa, estao 
apresentados graficamente na Figura 5.1. 

Como se observa, a serie nao apresenta tendencia de crescimento ou de decres- 
cimento, mas apenas um comportamento erratico, dispensando, assim, a necessidade 
de inclusao de uma variavel associada a passagem do tempo. 

De posse de todas essas informagoes foi posslvel montar uma planilha com os 
dados necessarios para a construgao do modelo de regressao multipla dentro do qual 
o comportamento da variavel dependente quantidade de ligagoes em um dia especlfico 
poderia estar sendo explicado - inicialmente - por 10 variaveis independentes: 

• dia da semana; 

• feriado - se o dia especifico consiste em um feriado ou nao; 
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• chegada (C) - quantas contas estao previstas para chegarem ao cliente no dia es¬ 
pecffico; 

• vencimento (V) - quantas contas vencem no dia especffico; 

• V - 3 - para quantas contas aquele dia especffico corresponde a tres dias antes 
do vencimento, ou, em outras palavras, quantas contas irao veneer tres dias 
depois da data especffica; 

• V - 2 - para quantas contas aquele dia especffico corresponde a antevespera 
do vencimento; 

• V - 1 - para quantas contas aquele dia especffico corresponde a vespera do 
vencimento; 

• V + 1 - para quantas contas aquele dia especffico corresponde a um dia depois 
do vencimento, ou, em outras palavras, quantas contas terao vencido na vespera 
daquela data especffica; 
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Quantidade de Iiga9des por dia (produto 103) 



Dia 


FIGURA 5.1 Quantidade de liga^oes por dia (julho de 2005 a mar^o de 2006). 

• V + 2 - para quantas contas aquele dia especifico corresponde a dois dias depois 
do vencimento; 

• V + 3 - para quantas contas aquele dia especifico corresponde a tres dias depois 
do vencimento. 

A Tabela 5.3 apresenta, em parte, a planilha de dados em questao. 

Cabe salientar que as variaveis independentes dia da semana e feriado foram traba- 
lhadas como sendo do tipo dummy (0 = Nao; 1 = Sim). Foi criada uma variavel dummy 
para a ocorrencia de feriado e seis variaveis dummy para o dia da semana: domingo, 
segunda-feira, ter 9 a-feira, quarta-feira, quinta-feira e sexta-feira. 

Dessa forma, para o dia 28/02/2006, por exemplo, a variavel dummy feriado as- 
sumiu o valor 1, por se tratar da terca-feira de Carnaval; as variaveis dummy domingo, 
segunda-feira, quarta-feira, quinta-feira e sexta-feira assumiram o valor 0, enquanto 
a variavel dummy ter 9 a-feira assumiu o valor 1. Seguindo o mesmo raciocinio, para o 
dia 04/02/2006, por exemplo, a variavel dummy feriado assumiu o valor 0; e as variaveis 
dummy domingo, segunda-feira, ter 9 a-feira, quarta-feira, quinta-feira e sexta-feira 
assumiram o valor 0, caracterizando, portanto e por exclusao, o dia em questao como 
sendo um sabado. 

Apos a regressao ser executada, a primeira analise aponta para um modelo com 
um potencial de explica^o muito bom, ja que o R-quadrado ajustado indica que 84% 
da varia 5 ao dos dados pode ser explicada. Neste caso, pode-se inferir a real influencia 
das variaveis explicativas na quantidade de I igardes e aceitar o modelo com elevada 
confianca, ja que o nlvel de significance da estatistica F foi igual a 0,0000. 

A contribuigao de cada variavel independente e expressa pelos coeficientes es- 
timados, juntamente com os respectivos p-values (ou valores p), apresentados na 
Tabela 5.4. 
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TABELA 5.3 Planilha de dados do modelo de 

regressao 

multipla para prever a 

quantidade de ligagoes (fevereiro 

de 2006) 


Dia 

Dia da 

semana 

Feriado? 

Chegada (C) 

V - 3 

V - 2 

V-1 

Vencimento 

(V) 

V + 1 

V + 2 

V + 3 

Liga^oes 

01/02/200 

4 a -feira 

Nao 

2.506.038 

2.142.215 

- 

- 

2.136.197 

- 

- 

- 

199.908 

02/02/200 

5 a -feira 

Nao 

- 

- 

2142.215 

- 

- 

2.136.197 

- 

- 

195.497 

03/02/200 

6 a -feira 

Nao 

1.122.892 

- 

- 

2.142.215 

- 

- 

2.136.197 

- 

185.542 

04/02/200 

Sabado 

Nao 

- 

1.918.920 

- 

- 

2.142.215 

- 

- 

2.136.197 

116.657 

05/02/200 

Domingo 

Nao 

- 

- 

1.918.920 

- 

- 

2.142.215 

- 

- 

65.558 

06/02/200 

2 a -feira 

Nao 

587.354 

- 

- 

1.918.920 

- 

- 

2.142.215 

- 

216.104 

07/02/200 

3 a -feira 

Nao 

1.484.473 

- 

- 

- 

1918.920 

- 

- 

2.142.215 

221.831 

08/02/200 

4 a -feira 

Nao 

1.238.119 

2.177.552 

- 

- 

- 

1.918.920 

- 

- 

213.584 

09/02/200 

5 a -feira 

Nao 

801.143 

- 

2.177.552 

- 

- 

- 

1.918.920 

- 

200.359 

10/02/200 

6 a -feira 

Nao 

- 

798.630 

- 

2.177.552 

- 

- 

- 

1.918.920 

185.882 

11/02/200 

Sabado 

Nao 

- 

- 

798.630 

- 

2.177.552 

- 

- 

- 

122.764 

12/02/200 

Domingo 

Nao 

- 

1.456.656 

- 

798.630 

- 

2.177.552 

- 

- 

72.033 

13/02/200 

2 a -feira 

Nao 

- 

- 

1.456.656 

- 

798.630 

- 

2.177.552 

- 

197.231 

14/02/200 

3 a -feira 

Nao 

432.670 

801.143 

- 

1.456.656 

- 

798.630 

- 

2.177.552 

181.092 

15/02/200 

4 a -feira 

Nao 

709.437 

- 

801.143 

- 

1.456.656 

- 

798.630 

- 

188.990 

16/02/200 

5 a -feira 

Nao 

- 

- 

- 

801.143 

- 

1.456.656 

- 

798.630 

174.253 

^17/02/200 

6 a -feira 

Nao 

- 

432.670 

- 

- 

801.143 

- 

1.456.656 

- 

156.317/ 

(Continua) 
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TABELA 5.3 Planilha de dados do modelo de regressao multipla para prever a quantidade de ligagoes (fevereiro de 2006) (cont.) 


Dia 

Dia da 

semana 

Feriado? 

Chegada(C) 

V - 3 

CN 

1 

> 

V-1 

Vencimento 

(V) 

V+ 1 

V + 2 

V + 3 

Liga^oes 

1 8/02/200 

Sabado 


Nao. 

- 

432.670 

- 

- 

801.143 

- 

1.456.656 

105.306 

19/02/200 

Domingo 

Nao 

- 

- 

- 

432.670 

- 

- 

801.143 

- 

63.827 

20/02/200 

2 a -feira 

Nao 

- 

709.437 

- 

- 

432.670 


- 

801.143 

171.921 

21/02/200 

3 a -feira 

Nao 

- 

- 

709.437 

- 

- 

432.670 

- 

- 

155.035 

22/02/200 

4 a -feira 

Nao 

2.136.197 

- 

- 

709.437 

- 

- 

432.670 

- 

152.537 

23/02/200 

5 a -feira 

Nao 

1.555.097 

- 

- 

- 

709.437 

- 

- 

432.670 

140.621 

24/02/200 

6 a -feira 

Nao 

2.506.038 

- 

- 

- 

- 

709.437 

- 

- 

139.056 

25/02/200 

Sabado 

Nao 

- 

- 

- 

- 

- 

- 

709.437 

- 

86.746 

26/02/200 

Domingo 

Nao 

- 

2.136.197 

- 

- 


- 

- 

709.437 

49.520 

27/02/200 

2 a -feira 

Sim 

- 

- 

2.136.197 

- 

- 

- 

- 

- 

85.790 

28/02/200 

3 a -feira 

Sim 

- 

- 

- 

2.136.197 

- 

- 

- 

- 

75.801 
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TABELA 5.4 Coeficientes estimados das variaveis independentes do modelo 

de regressao multipla 

para prever a quantidade de ligagoes na 

data em questao 


Coeficientes 

valor-P 

Intersegao 

92.545 

0,00 

Domingo 

(41.228) 

0,00% 

2 a -feira 

79.882 

0,00% 

3 a -feira 

75.694 

0,00% 

4 a -feira 

73.812 

0,00% 

5 a -feira 

73.068 

0,00% 

6 a -feira 

61.763 

0,00% 

Feriado 

(71.088) 

0,00% 

Chegada (C) 

0,0035 

13,57% 

V - 3 

0,0052 

2,21% 

V - 2 

0,0042 

7,81% 

V- 1 

0,0068 

0,47% 

Vencimento (V) 

0,0074 

0,19% 

V+ 1 

0,0071 

0,30% 

V + 2 

0,0064 

0,75% 

V + 3 

0,0009 

68,36% 

s -- 


_ 


Na Tabela 5.4, observa-se um elevado valor p para a variavel V + 3, indicando 
que a quantidade de contas para as quais aquele dia especffico corresponde a tres dias 
depois do vencimento nao pode ser vista como um bom previsor para a quantidade de 
ligagoes na data em questao. 

Assim sendo, o efeito pos-vencimento parece, entao, acabar ja no segundo dia apos 
o vencimento da conta, ja que os valores pdeV+leV + 2 sao baixos; em outras 
palavras, os clientes ainda ligam um ou dois dias apos o vencimento da conta, mas esse 
impacto desaparece a partir do terceiro dia, deixando de ser significativo. 

A variavel V + 3 pode estar contaminando o modelo e convem retira-la da analise. 
Sem ela, um novo modelo de regressao foi gerado, continuando com um potencial 
muito bom e ate ligeiramente melhor, ja que 85% da variagao dos dados pode ser 
explicada por ele. 

Os coeficientes de regressao estimados, juntamente com os respectivos valores p, 
estao apresentados na Tabela 5.5. 

Como pode ser observado, os valores p das outras variaveis mudaram muito 
pouco e, para este modelo modificado, apenas duas variaveis apresentam coeficientes 
angulares com valores p superiores a 5%: Chegada (C) e V - 2. 
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TABELA 5.5 Coeficientes estimados das variaveis independentes do modelo 
(sem V + 3) de regressao multipla para prever a quantidade de ligagoes na data 
em questao 


Coeficientes 

valor-P 

Intersegao 

92.693 

0,00% 

Domingo 

(41.103) 

0,00% 

2 a -feira 

79.874 

0,00% 

3 a -feira 

75.900 

0,00% 

4 a -feira 

73.908 

0,00% 

5 a -feira 

73.082 

0,00% 

6 a -feira 

61.898 

0,00% 

Feriado 

(71.345) 

0,00% 

Chegada (C) 

0,0034 

14,30% 

V - 3 

0,0052 

2,11% 

V - 2 

0,0042 

7,86% 

V- 1 

0,0072 

0,13% 

Vencimento (V) 

0,0077 

0,06% 

V+ 1 

0,0071 

0,30% 

V + 2 

0,0060 

0,71% 



_^ 


Em relagao a esta ultima variavel, ela sera mantida no modelo porque o valor p 
do seu coeficiente esta abaixo de um valor tambem bastante utilizado como nfvel de 
significancia nos testes de hipoteses (10%) e pelo fato de as variaveis V - 3 e V - 1 
apresentarem-se tambem como bons previsores. Ora, nao faria muito sentido man ter no 
modelo V - 3, retirar V - 2 e manter V - 1, pois a continuidade estaria sendo quebrada; 
seria o mesmo que afirmar que os clientes ligam para o produto em questao tres dias 
antes do vencimento e na sua vespera, mas nao na antevespera! 

Ja a variavel Chegada (C) tambem sera mantida no modelo, apesar de apresentar 
um valor p maior do que 10%, pelo fato de ter sido identificada, na opiniao subjetiva e 
calcada na experiencia da equipe de planejamento, como crftica e muito impactante na 
demanda de ligagoes, contrariando o resultado objetivo que indica um coeficiente nao 
significativo para essa variavel. Tambem e valido ressaltar que seu valor p (14,30%) 
nao ficou tao acima de um nfvel de significancia considerado bom (5%). 

Essa intervengao subjetiva em um modelo tao objetivo traz a tona a ja mencionada 
importancia da utilizagao da opiniao de especialistas na construgao de um processo de 
previsao, mesmo nos casos em que esse e essencialmente sistematico. 
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Dessa forma, a retirada da variavel V + 3 foi benefica ao modelo, tornando-o mais 
parcimonioso e permitindo uma interpretagao dos coeficientes obtidos. 

O primeiro coeficiente corresponde ao valor da variavel dependente quantidade 
de ligagoes quando todas as variaveis independentes assumirem o valor zero; ou seja, 
quando domingo, segunda-feira, terga-feira, quarta-feira, quinta-feira, sexta-feira, 
feriado, C, V - 3, V - 2, V - 1, V, V+leV + 2 forem zero - ou, ainda, quando o dia 
for sabado, nao feriado e nao consistir da data de chegada, nem de nenhuma data dentro 
do intervalo de seis dias em torno do vencimento de alguma conta, a quantidade de 
ligagoes demandadas prevista e de 92.693. 

Esse numero pode ser encarado como uma quantidade basica, a qual devem ser 
adicionados os efeitos da ocorrencia das variaveis independentes, quantificados pelos 
seus respectivos coeficientes angulares. 

Assim, o fato de o dia ser domingo contribui com uma diminuigao de 41.103 
ligagoes previstas em relagao ao numero basico (92.693); a segunda-feira espera, 
segundo o modelo, 79.874 ligagoes a mais do que o dia basico (sabado); a terga-feira, 
75.900 ligagoes a mais do que o sabado; a quarta-feira, 73.908 ligagoes a mais; a 
quinta-feira, 73.082 a mais; e a sexta-feira, 61.898 ligagoes a mais do que o dia basico. 

Com essas informagoes, e possfvel elaborar um grafico (Figura 5.2) que permite 
uma melhor visualizagao da sazonalidade do dia da semana. 

A evolugao da demanda de ligagoes ao longo da semana, considerando dias sent 
a ocorrencia de feriados e nao proximos da chegada ou do vencimento de contas, 
comega com um valor bem baixo no domingo, apresentando um grande aumento na 
segunda-feira; a demanda vai caindo suavemente ao longo da semana, tendo uma 
queda um pouco mais brusca na sexta-feira e bem mais intensa no sabado, antes de 
voltar ao patamar de domingo. 


Sazonalidade do dia semana - produto 103 



FIGURA 5.2 Comportamento sazonal do dia da semana - modelo sem V + 3. 
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De forma analoga aos dias da semana, o fato de um dia ser feriado reduz em 
71.345 a quantidade prevista de Iigardes; e cada conta telefonica prevista para chegar 
a casa do cliente naquele dia aumenta em 0,0034 a previsao da quantidade de ligagoes 
demandadas para o dia em questao. Assim, se, por exemplo, 100 mil contas estiverem 
chegando em um dia especifico, e de se esperar que isso implique um aumento de 
100.000 X 0,0034 = 340 ligagoes no contingente diario. Se esse numero for examinado 
de forma percentual, e possfvel concluir que, em media, 0,34% (= 0,0034) dos clientes 
ligam para a empresa telefonica no dia em que recebem a conta. 

Em relagao aos dias proximos a data do vencimento, e analogo o raciocinio que 
explica as consequencias: cada conta vencendo tres dias depois, ou seja, para as quais a 
data em questao representa V - 3, aumenta em 0,0052 a quantidade de ligagoes; cada 
conta vencendo dois dias depois aumenta em 0,0042 a quantidade de ligagoes; cada conta 
vencendo no dia seguinte aumenta em 0,0072 essa quantidade; cada conta vencendo 
no dia em questao aumenta em 0,0077 a mesma quantidade; cada conta vencendo 
na vespera aumenta em 0,0071 a quantidade; e cada conta vencendo na antevespera 
aumenta em 0,0060 a quantidade prevista de ligagoes. 

Naturalmente, os coeficientes nao tem apenas significado isolado, mas tambem 
podem e devem ser analisados em conjunto: se as consequencias dos efeitos dia da 
semana, ocorrencia de feriado e data crftica em relagao a chegada ou ao vencimento da 
conta forem incorporadas ao valor previsto para um dia basico, sera possfvel estabelecer 
uma equagao para prever a quantidade demandada de ligagoes em determinado dia em 
fungao das caracterfsticas da data em questao. Desta maneira, a quantidade de ligagoes 
pode ser prevista conforme a Equagao 5.1. 

EQUAGAO 5.1 Equagao do modelo para prever a quantidade de ligagoes 

Quantidade prevista de ligagoes = 92.693 - 41.103 x domingo+ 

79.874 x segunda-feira + 75.900 x terga-feira + 

73.908 x quarta-feira + 73.082 x quinta-fcira + 61.898 x sexta-fcira 
-71.345 x feriado + 0,0034 xC+0,0052 x V - 3 + 

0,0042 x V - 2 + 0,0072 x V -1 + 

0,0077 x V + 0,0071 x V +1 + 0,0060x V + 2 

• as sete primeiras variaveis sao do tipo dummy, assumindo o valor 0 (nao) ou 1 (sim); 

• as sete ultimas representam, respectivamente: a quantidade de contas chegando no 
dia; vencendo tres dias depois; dois dias depois; um dia depois; no proprio dia; um 
dia antes; e dois dias antes. 

De posse dessa equagao, e possfvel calcular quanto o modelo de regressao teria 
previsto para a quantidade de ligagoes em cada um dos dias acerca dos quais o his- 
torico disponibiliza o valor real desta grandeza. Desta forma, pode-se comparar os 
valores que teriam sido previstos com os que efetivamente ocorreram, de maneira 
que um erro de previsao possa ser mensurado. Esta comparagao pode ser visualizada 
na Figura 5.3. 
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Ouantidade de ligacoes (reais x previstas) por dia (produto 103) 



17/07/2005 16/08/2005 15/09/2005 15/10/2005 14/11/2005 14/12/2005 13/01/2006 12/02/2006 14/03/2006 


Dia 

FIGURA 5.3 Valores reais multiplicados por uma constante e previstos pelo modelo para a quan- 
tidade de ligagoes por dia (julho de 2005 a mar^o de 2006). 

A diferenga entre os valores reais e os gerados pelo modelo tambem e conhecida 
como resfduo. Para ser usado um modelo de regressao miiltipla, e necessario supor 
que a correlagao entre os residuos, ao longo do espectro das variaveis independentes, 
e zero; ou seja, que os residuos sao independentes entre si, nao apresentando, conse- 
quentemente, autocorrelagao serial. 11 

Uma forma de verificar a validade dessa suposigao e por meio do teste de Dur¬ 
bin-Watson. No modelo que esta sendo analisado, foi encontrado o valor 1,838 para 
a Estatistica DW. Esta estatistica obtida deve ser comparada com valores criticos que 
variam em fungao do nivel de significancia do teste, do tamanho da amostra e do 
numero de variaveis independentes. 12 

Com 5% de significancia, 234 dados na amostra e 14 variaveis independentes, o 
valor critico inferior (d L ) obtido foi 1,661 e o valor critico superior (d ([ ) encontrado 
foi 1,914. 13 

Como o valor obtido para a Estatistica DW (1,838) se situa entre os dois valores 
criticos encontrados, o teste para verificar a existencia de autocorrelagao dos residuos e 
nao conclusivo: nao e possivel garantir que os residuos sao independentes, mas tambem 
e incorreto afirmar que existe uma relagao de dependencia entre eles. 14 

Sendo assim, parece valido examinar o grafico do residuo contra cada variavel 
independente a fim de procurar uma suposta autocorrelagao e a identificagao de uma 


11. Corral'; Paulo; Dias Filho (2007). 

12. Corral'; Paulo; Dias Filho (2007). 

13. Valores encontrados em tabela disponivel no enderego <http://www.stanford.edu/~clint/bench/dwcrit. 
htm>. Acesso em 24/01/2009. 

14. Corrar; Paulo; Dias Filho (2007). 
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eventual lei de formagao nos valores dos reslduos a medida que o valor da variavel 
independente avanga. A Figura 5.4 ilustra o caso da variavel dummy quarta-feira. 

Como pode ser observado, nao parece existir uma lei de formagao; ou seja, a magni¬ 
tude dos reslduos nao parece depender do valor da variavel independente (quarta-feira). 

Os reslduos apresentaram comportamento semelhante em relagao as outras va¬ 
riaveis dummy (feriado e as referentes aos outros dias da semana). E visto que o 
modelo trabalha com 14 variaveis independentes, seria um tanto inviavel apresentar 
aqui tal quantidade de graficos. Dessa forma, serao apresentadas graficamente apenas 
as relagoes entre os reslduos e as outras variaveis independentes mais marcantes: 
Chegada e Vencimento. Tais redoes podem ser visualizadas nas Figuras 5.5 e 5.6, 
respectivamente. 

O comportamento dos reslduos em relagao as outras variaveis (V - 3, V - 2, V - 1, 
V + 1 e V + 2) e bastante similar ao verificado nas Figuras 5.5 e 5.6: parecem estar 
distribuldos aleatoriamente em torno do eixo das abscissas, nao apresentando nenhuma 
tendencia. 

Tudo isso gera uma seguranga maior para os resultados do modelo, ja que o teste 
para a verificagao de autocorrelagao serial dos reslduos foi nao conclusivo. E claro que 
esta ultima analise e meramente visual, carecendo de um cientificismo maior. Mas e 
importante proceder dessa maneira principalmente em situagoes em que os criterios 
objetivos apontam para uma indefinigao, como foi o caso. 

Ja em relagao a qualidade preditiva do modelo, foi posslvel perceber, na Figura 5.3, 
a menos de alguns pontos especlficos, uma boa aderencia visual dos valores previstos 
em relagao aos dados reais. Esta aderencia foi aqui formalizada pelo MAPE (Mean 
Absolute Percentual Error), ou erro medio percentual absoluto. 


Plotagem de reslduos 
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FIGURA 5.4 Reslduos do modelo X variavel dummy quarta-feira. 
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Plotagem de residuos 



FIGURA 5.5 Residuos do modelo X variavel Chegada. 


Plotagem de residuos 



FIGURA 5.6 Residuos do modelo X variavel Vencimento. 


A tftulo de ilustragao, a Tabela 5.6 apresenta dados relativos ao mes de dezembro 
de 2005, nos quais constam a quantidade real de ligagoes, os valores que o modelo de 
regressao teria previsto para esta grandeza, assim como os erros percentuais absolutos 
oriundos desta previsao. Pode-se calcular que o modelo teria errado em 15,5%, em 
media, na previsao dos volumes diarios de ligagoes durante o mes em questao. 
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TABELA 5.6 Valores reais, previstos e erros do modelo de regressao multipla 
para prever a quantidade de liga^oes (dezembro de 2005) 


Dia 

Dia da semana 

Feriado? 

Liga^oes 

Liga^oes (previstas Erro de previsao 
pelo modelo) do modelo 

01/12/2005 

5 a -feira 

Nao 

189.179 

193.288 

2,2% 

02/12/2005 

6 a -feira 

Nao 

179.607 

178.716 

0,5% 

03/12/2005 

Sabado 

Nao 

119.738 

121.003 

11% 

04/12/2005 

Domingo 

Nao 

63.118 

77.992 

23,6% 

05/12/2005 

2 a -feira 

Nao 

209.398 

203.208 

3,0% 

06/12/2005 

3 a -feira 

Nao 

203.525 

197.677 

2,9% 

07/12/2005 

4 a -feira 

Nao 

201.497 

186.681 

7,4% 

08/12/2005 

5 a -feira 

Nao 

163.731 

190.634 

164% 

09/12/2005 

6 a -feira 

Nao 

165.708 

175.358 

5,8% 

10/12/2005 

Sabado 

Nao 

108.033 

112.476 

41% 

11/12/2005 

Domingo 

Nao 

65.321 

71.680 

97% 

12/12/2005 

2 a -feira 

Nao 

204.689 

201.258 

1,7% 

13/12/2005 

3 a -feira 

Nao 

202.289 

196.754 

2,7% 

14/12/2005 

4 a -feira 

Nao 

187.186 

188.341 

0,6% 

15/12/2005 

5 a -feira 

Nao 

168.579 

185.651 

10,1% 

16/12/2005 

6 a -feira 

Nao 

155.912 

170.652 

95% 

17/12/2005 

Sabado 

Nao 

102.895 

109.900 

68% 

18/12/2005 

Domingo 

Nao 

58.292 

59.078 

1,3% 

19/12/2005 

2 a -feira 

Nao 

170.636 

180.520 

5,8% 

20/12/2005 

3 a -feira 

Nao 

155.824 

172.655 

10,8% 

21/12/2005 

4 a -feira 

Nao 

145.547 

170.260 

1 70% 

22/12/2005 

5 a -feira 

Nao 

127.727 

169.413 

32,6% 

23/12/2005 

6 a -feira 

Nao 

116.546 

155.684 

33,6% 

24/12/2005 

Sabado 

Nao 

69.839 

93.699 

342% 

25/12/2005 

Domingo 

Sim 

39.038 

118.895) 

1 484% 

26/12/2005 

2 a -feira 

Nao 

152.909 

172.567 

12,9% 

27/12/2005 

3 a -feira 

Nao 

158.507 

168.593 

64% 

28/12/2005 

4 a -feira 

Nao 

164.109 

166.601 

1,5% 

29/12/2005 

5 a -feira 

Nao 

154.690 

165.775 

7,2% 

30/12/2005 

6 a -feira 

Nao 

124.406 

154.591 

24,3% 

31/12/2005 

Sabado 

Nao 

67.541 

92.693 

37,2% 
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Mas dezembro e um mes atipico, com um comportamento geralmente mais diffcil 
de prever. Tal afirmacao pode ser verificada neste exemplo particular a partir da 
Figura 5.7, no qual e possrvel perceber que todos os outros meses relacionados, com 
excegao de janeiro, outro periodo atipico, apresentam erros medios de previsao pelo 
menos quatro pontos percentuais menores do que o erro medio de dezembro de 2005. 

Em media, o modelo esta incorrendo num erro de previsao da ordem de 10,98%. 
Embora nao possa ser considerado um valor elevado, ele e substancialmente maior que 
os 5% em que a equipe de planejamento costuma incorrer na previsao de demanda do 
produto em questao. No entanto, cabem aqui algumas observagoes. 

E valido lembrar, por exemplo, que, durante o processo padrao de previsao, 
os analistas recebem os valores previstos pelo software e acrescentam a eles os 
impactos esperados por eventos e premissas, como o comportamento do ciclo de 
contas e a ocorrencia de feriados, assim como outros eventos especiais relacionados 
ao produto, como piano de minutos, campanhas na mldia, crescimento de base e ou¬ 
tros potenciais aspectos impactantes. Os 5% de erro de previsao so sao conseguidos 
apos a consideragao de todas essas subjetividades e dentro de um processo diario de 
revisao da previsao. 

Ora, o modelo de regressao proposto visa substituir apenas os valores sugeridos 
pelo software, acrescidos da influencia dos feriados e do comportamento dos ciclos 
de contas; nunca foi a ideia substituir o papel do analista encarregado da previsao - 
dotado de toda a sua valiosa experiencia e subjetividade mas, sim, de fornecer-lhe um 
subsldio complementar para a tomada de decisao. Assim, e visto que o analista estara 
apto a acrescentar aos valores sugeridos pelo modelo de regressao os impactos que 
julgar relevantes, e mais justo comparar os erros de previsao do modelo aqui proposto 


Erro medio de previsao do modelo (mes a mes) 



Mes/ano 

FIGURA 5.7 Erro medio mensal do modelo de regressao multipla para prever a quantidade de 
ligagoes (julho de 2005 a mar^o de 2006). 















































































(jhT) PARTE | II Metodos de Previsao 


ELSEVIER 


com os erros oriundos do modulo de previsao do software, com o intuito de verificar 
a ocorrencia ou nao de melhorias no processo de previsao. 

Assim sendo, o MAPE do software para o mesmo perfodo analisado foi de 13,54%, 
superior, portanto, aos 10,98% do modelo de regressao. 

Alem disso, uma analise mais detalhada dos erros de previsao do modelo proposto 
revela que alguns dos dias mais diffceis de prever e, portanto, com erros de previsao 
potencialmente maiores, consistem em feriados ou datas que se comportam quase 
como um feriado. Por exemplo, os dias 24/12/2005 e 31/12/2005 nao foram clas- 
sificados como feriados, mas certamente tambem nao sao dias comuns e devem gerar 
uma demanda de ligagoes inferior ao de um dia padrao. O problema e que o modelo de 
regressao nao levou em conta isso, ja que era posslvel classificar um dia apenas como 
sendo feriado ou nao; e visto que esses dias nao eram efetivamente feriados, foram 
tratados pelo modelo como dias comuns. Por este motivo, muito provavelmente, suas 
demandas foram superestimadas e seus erros de previsao foram bem acima da media: 
34,2% e 37,2%, respectivamente, conforme ja apresentado na Tabela 5.6. 

Problemas semelhantes, mas talvez em menor escala, certamente ocorreram em 
alguns dias do ano que consistiram em feriados locais, mas nao nacionais. Como o 
modelo propoe uma previsao consolidada da demanda, so foi possfvel considerar os 
feriados nacionais; por exemplo, o dia 20/01/2006, um feriado local na Cidade do 
Rio de Janeiro, teve uma redugao na demanda nacional provocada pela redugao na 
demanda da cidade mas, ao ser tratado como um dia normal pelo modelo, teve sua 
demanda superestimada. Este problema tambem se observa nos feriados prolongados 
e nos dias que se situam entre um feriado e o fim de semana, em que a demanda foge 
do padrao. 

Um problema ainda mais grave ocorreu com o dia 25/12/2005. Seu erro de previsao 
foi de quase 150%, ja que a demanda prevista assumiu um valor negativo, -18.895 
ligagoes, conforme a Tabela 5.6. Isto aconteceu porque o dia em questao apresentava 
dois efeitos negativos e redutores de demanda, que foram somados pelo modelo de 
natureza aditiva: ser domingo e ser feriado. Mas, na pratica, o efeito de um dia ser 
feriado certamente e reduzido se ele cair em um domingo ou sabado; e isto nao foi 
levado em consideragao, ate porque so e possfvel considerar um dia como sendo feriado 
ou nao, e nao como uma especie de feriado moderado. 

Na verdade, esses problemas aumentaram os erros de previsao, mas isso nao 
quer dizer que o processo futuro de previsao tera que incorrer nesses erros de grande 
magnitude quando se tratar de dias com comportamentos diferentes. Na pratica, quando 
os analistas de previsao estiverem utilizando o modelo aqui proposto, eles poderao 
manipular os valores da variavel dummy feriado para os dias diferentes do padrao. 

Assim, por exemplo, o dia 24/12 pode receber manualmente o valor 0,7 para a sua 
variavel Feriado, em vez de ser tratado como um dia normal; o dia 20/01 pode receber 
o valor 0,4 da mesma forma; um feriado que caia num domingo pode ter o seu valor 
da variavel feriado reduzido de 1 para 0,6. 

Depois dessa manipulagao, basta considerar as informacbes do setor de faturamento 
da empresa telefonica sobre os ciclos de contas e usar a Equagao 1 apresentada ante- 
riormente para prever a quantidade de ligagoes para o produto em questao, a cada dia. 
Uma vez gerados os outputs da equagao, os analistas podem incorporar os impactos 
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esperados pelos eventos especiais, como campanhas de rmdias, novos servigos etc. e, 
desta forma, reduzir substancialmente os erros de previsao. 

EXERCICIOS 

1 . A planilha Exercicio 1 apresenta algumas caracteristicas de alguns imoveis que 
foram vendidos por uma corretora: a quantidade de banheiros, a area (em m 2 ) e o 
prego de venda (em US$). 

a. Desenvolva um modelo de Regressao Multipla capaz de explicar o prego de 
venda do imovel em fungao da quantidade de banheiros e da area. 

b. Qual a capacidade explicativa deste modelo? 

c. Como se comportam as variaveis independentes como explicadoras do prego 
de venda? Trabalhe com 10% de significancia. 

d. Aplique a equagao gerada pelo modelo para “prever” o prego de venda de cada 
imovel. 

e. Calcule o erro percentual absoluto medio dessas previsoes. Compare com o 
resultado encontrado na letra b. 

f. Segundo o modelo, qual seria o prego de venda de um imovel com dois ba¬ 
nheiros e 200 m 2 ? 

g. Segundo o modelo, qual seria o prego de venda de um imovel com tres banheiros 
e 200 m 2 ? Qual e o valor de um banheiro adicional? 

h. Segundo o modelo, qual seria o prego de venda de um imovel com tres banheiros 
e 300 m 2 ? Qual e o valor de 100 m 2 adicionais? 

2. Uma cadeia de lojas que vende CDs de musica esta tentando determinar a eficiencia 
de sua propaganda via radio e televisao. Foram coletados dados de algumas lojas 
da cadeia sobre vendas e gastos mensais em propaganda tanto por radio quanto 
por televisao. Os dados estao na planilha Exercicio 2. 

a. Existe evidencia, a 5% de significancia, de alguma relagao entre os gastos em 
cada tipo de propaganda e as vendas? 

b. Que tipo de propaganda parece ser mais eficiente em relagao ao custo? 

c. Qual o nivel de explicagao do modelo como um todo? 

3. Uma empresa que presta servigo de jardinagem esta interessada em prever a 
quantidade mensal media que sera gasta com servigos de jardinagem em uma certa 
vizinhanga. Uma amostra de vizinhangas e selecionada, e a quantidade media gasta 
e registrada, juntamente com o tamanho medio dos terrenos (em lotes), a renda 
familiar media e o valor medio das casas. Os dados estao na planilha Exercicio 
3. Desenvolva e avalie um modelo de Regressao Multipla para estimar o que 
impacta, a 5% de significancia, os gastos medios por residencia em servigos de 
jardinagem. 
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Capftulo 6 


Prevendo o tempo medio 
de atendimento em um call 
center por meio de um modelo 
de Regressao Multipla 


Marco Aurelio Carino Bouzada 

MADE/UNESA 


Objetivos de aprendizagem do capftulo 







Saber construir um modelo causal, capaz de explicar o comportamento de uma 
grandeza real com base no comportamento de outras variaveis. 

Entender as vantagens de um modelo causal sobre um modelo de serie temporal. 
Aprender a quantificar variaveis qualitativas e a incluf-las no modelo por meio de 
variaveis dummy. 

Conseguir aferir a capacidade de explicagao do modelo. 

Conseguir aferir a capacidade de explicagao de cada variavel e entender como cada 
uma impacta a variavel dependente. 

Saber refinar o modelo, tornando-o mais parcimonioso com a retirada de variaveis 
pouco explicativas. 

Entender como obter e interpretar a equagao capaz de prever a variavel dependente. 
Aprender a mensurar o erro de previsao decorrente da utilizagao de modelo desen- 
volvido. 


INTRODU^AO 

Tendo sido conceituados os call centers no Capftulo 5, assim como informados os seus 
desafios operacionais, destacou-se que o custo de pessoal nessa industria representa 
aproximadamente 70% do total. A materia-prima para o processo de gerenciamento 
eficiente (nfveis satisfatorios de service) a custos razoaveis) de um call center consiste 
na previsao da demanda - tanto em termos da quantidade de ligagoes quanto do tempo 
de atendimento - e da consequente carga de trabalho. 

Dessa forma, outra caracterfstica operacional que ganha importancia e precisa ser 
estimada com acuracia para o correto dimensionamento da equipe de atendimento 
consiste no tempo medio de atendimento, tambem conhecida como TMA. 
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Uma previsao pode ser necessaria para dimensionar as opera 9 oes cotidianas 
ou para lidar com situa 9 oes especiais, como, por exemplo, a implementa 9 ao de 
uma nova tecnologia capaz de afetar nao so o volume, mas tambem o padrao das 
chamadas. 1 

Como o proposito real do processo de previsao e predizer quantos agentes sao 
necessarios para lidar com o volume de chamadas, uma previsao precisa incluir - para 
ser completa - nao apenas o volume de liga 9 oes, mas tambem a grandeza tempo de 
atendimento. Para prever o tamanho da equipe e os pianos de horarios posteriormente, 
e necessaria uma fotografia da carga de trabalho, que consiste no numero de chamadas 
multiplicado pelo TMA. 2 

O TMA e um dos alvos atacados pelo metodo proposto por Shen e Brown (2002), 
que e aplicado - entre outros propositos - para modelar o padrao (dependente da hora 
do dia) do TMA das chamadas de clientes recebidas no call center de um banco is- 
raelense. O entendimento de tal comportamento variavel desta grandeza e, na opiniao 
dos autores, essencial para compreender o ambiente operacional de um sistema, e 
tambem para prever dinamicamente a sua carga de trabalho futura. 

O PROCESSO DE PREVISAO DO TMA NO CALL CENTER 
ESTUDADO SEM A REGRESSAO MULTIPLA 

Conforme dito no Capltulo 5, a Gerencia de Planejamento de Trafego da empresa 
preocupa-se sobretudo em conseguir melhorar as previsoes de demanda, mas nao so 
em termos da quantidade de liga 9 oes, como tambem em termos do TMA. 

E a impressao da gerencia quanto ao TMA e a mesma no que se refere ao volume 
de liga 9 oes: o tempo de atendimento para o produto 103 e influenciado pelos eventos 
de chegada da conta ao destino e pelo dia do seu vencimento, datas em que haveria 
um acesso diferenciado ao servi 90 por parte dos clientes. 

O processo de previsao para o TMA segue os mesmos tres passos, descritos no 
Capltulo 5, do processo de previsao do volume de chamadas. Mas os eventos es¬ 
peciais capazes de impactar a grandeza a ser prevista podem ser um pouco diferentes: 
mudan 9 as no procedimento (roteiro de atendimento), que podem levar o operador a 
falar mais ou menos durante o atendimento; mudan 9 as na Unidade de Resposta Audlvel 
(URA); entre outros. Mas, assim como no outro processo, aqui o comportamento dos 
ciclos das contas telefonicas dos clientes tambem consiste em uma informa 9 ao que 
ajuda a previsao. 

Sob os olhos de uma analise mais minuciosa, a magnitude do TMA tem variado 
muito ao longo do dia, assumindo, na parte da tarde, valores 35 a 40 segundos maiores 
que os tempos matutinos, segundo a coordena 9 ao da area. A impressao geral e que 
isso ocorre, principalmente, pela maior concentra 9 ao de operadores novos na parte 
da manha. Durante a noite, o TMA e ainda mais alto, principalmente - na opiniao da 
equipe de previsao - por causa do menor nlvel de servi 90 oferecido neste perlodo, 
o que pode fazer com que o cliente aguarde mais tempo na flla e reclame a respeito 


1. Reynolds (2005). 

2. Reynolds (2005); Wise (2006). 
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com o atendente antes de comegar a tratar do assunto principal, propriamente dito, 
aumentando o tempo de atendimento. 

O PROCESSO DE PREVISAO DOTMA COM A REGRESSAO 
MULTIPLA 

A ideia novamente consiste em poder contar com uma ferramenta que use as infor¬ 
magoes sobre as contas dos clientes no processo de elaboragao das previsoes. Assim 
como anteriormente, o metodo escolhido foi o da Regressao Multipla. 

As informagoes em termos da data de chegada ao destino e data de vencimento 
das contas foram consolidadas e parcialmente ja apresentadas nas Tabelas 5.1e 5.2, 
respectivamente. 

Para desenvolver o modelo de regressao multipla para prever o TMA em cada data, 
foi preciso ter acesso ao historico desta grandeza. Por motivo de sigilo estrategico, 
esses valores nao puderam ser revelados em sua plena magnitude; no entanto, a eles foi 
aplicada uma constante de proporcionalidade; os resultados, de julho de 2005 ate margo 
de 2006, fornecidos pela empresa, estao apresentados graficamente na Figura 6.1 . 

Como se observa, a serie nao apresenta tendencia de crescimento ou decres- 
cimento, mas apenas um comportamento erratico, dispensando, assim, a inclusao de 
uma variavel associada a passagem do tempo. 

De posse de todas essas informagoes foi possfvel montar uma planilha com os 
dados necessarios para a construgao do modelo de Regressao Multipla dentro do qual 
o comportamento da variavel dependente TMA em um dia especffico poderia es- 
tar sendo explicado - inicialmente - por nove variaveis independentes: (i) dia da 


TMA por dia (produto 103) 



Dia 


FIGURA 6.1 TMA por dia (julho de 2005 a mar^o de 2006). 
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semana; (ii) feriado - se o dia especffico consiste em um feriado ou nao; (iii) chegada 
(C) - quantas contas estao previstas para chegarem ao cliente no dia especffico; (iv) 
vencimento (V) - quantas contas vencem no dia especffico; (v) V - 3 - para quantas 
contas aquele dia especffico corresponde a tres dias antes do vencimento, ou, em outras 
palavras, quantas contas irao veneer tres dias depois da data especffica; (vi) V - 2 - para 
quantas contas aquele dia especffico corresponde a antevespera do vencimento; (vii) 
V - 1 - para quantas contas aquele dia especffico corresponde a vespera do vencimento; 
(viii) V + 1 - para quantas contas aquele dia especffico corresponde a um dia depois do 
vencimento, ou, em outras palavras, quantas contas terao vencido na vespera daquela 
data especffica; e (ix) V + 2 - para quantas contas aquele dia especffico corresponde 
a dois dias depois do vencimento. 

A Tabela 6.1 apresenta, em parte, a planilha de dados em questao. 

Cabe salientar que as variaveis independentes dia da semana e feriado foram traba- 
lhadas como sendo do tipo dummy (0 = Nao; 1 = Sim). Foi criada uma variavel dummy 
para a ocorrencia de feriado e seis variaveis dummy para o dia da semana: domingo, 
segunda-feira, terga-feira, quarta-feira, quinta-feira e sexta-feira. 

Dessa forma, para o dia 06/08/2005, por exemplo, a variavel dummy feriado as- 
sumiu o valor 0; e as variaveis dummy domingo, segunda-feira, terga-feira, quarta-feira, 
quinta-feira e sexta-feira assumiram o valor 0, caracterizando, portanto e por exclusao, 
o dia em questao como sendo um sabado comum (nao feriado). 

Apos a regressao ser executada, a primeira analise aponta para um modelo com 
um potencial de explicagao bom, ja que o R-quadrado ajustado indica que 79% 
da variagao dos dados podem ser explicados. Neste caso, pode-se inferir a real 
influencia das variaveis explicativas no tempo medio de atendimento e aceitar o 
modelo com elevada confianga, visto que o nfvel de significancia da estatfstica F 
foi igual a 0,0000. 

A contribuigao de cada variavel independente e expressa pelos coeficientes es- 
timados, juntamente com os respectivos p-values (ou valores p), apresentados na 
Tabela 6.2. Observa-se um altfssimo valor p para a variavel Chegada (C), indicando 
que a quantidade de contas chegando ao cliente naquele dia especffico nao pode ser 
vista como um bom previsor para o TMA na data em questao. 

Apesar de ser uma variavel crftica, ela deve ser retirada da analise, pois seu valor 
p esta muito acima do razoavel (10%) e porque o sinal do seu coeficiente angular 
(negativo) estaria indicando que o TMA em determinado dia diminuiria para cada 
conta que chegasse ao cliente no dia em questao, constituindo um efeito contrario ao 
esperado. 

As variaveis V-3eaV-2 tambem apresentam valores p altos, indicando que a 
quantidade de contas para as quais aquele dia especffico corresponde a dois ou tres 
dias antes do vencimento nao pode ser vista como um bom previsor para o TMA. 
O efeito pre-vencimento para o TMA parece, entao, comegar de forma significativa 
apenas na vespera (V - 1). 

Assim como a variavel Chegada (C), as variaveis V - 3 e V - 2 podem estar 
contaminando o modelo e convem retira-las da analise. Sem ela, um novo modelo de 
regressao foi gerado, continuando com um bom potencial, ja que 79% da variagao dos 
dados podem ser explicados por ele. 




TABELA 6.1 Planilha de dados do modelo de regressao multipla para prever o TMA (agosto de 2005) 


Dia 

Dia da semana 

Feriado? 

Chegada (C) 

V - 3 

V - 2 

V-1 

Vencimento (V) 

V + 1 

V + 2 

TMA (seg) 

01/08/2005 

2 a -feira 

Nao 

- 

2142.215 

- 

- 

2.136.197 

- 

- 

145,50 

02/08/2005 

3 a -feira 

Nao 

1.918.920 

- 

2.142215 

- 

- 

2.136.197 

- 

144,00 

03/08/2005 

4 a -feira 

Nao 

2.177.552 

- 

- 

2.142.215 

- 

- 

2.136.197 

144,00 

04/08/2005 

5 a -feira 

Nao 

- 

1.918.920 

- 

- 

2.142.215 

- 

- 

147,00 

05/08/2005 

6 a -feira 

Nao 

218.438 

- 

1.918.920 

- 

- 

2.142.215 

- 

144,75 

06/08/2005 

Sabado 

Nao 

- 

- 

- 

1.918.920 

- 

- 

2.142.215 

129,00 

07/08/2005 

Domingo 

Nao 

- 

- 

- 

- 

1.918.920 

- 

- 

113,25 

08/08/2005 

2 a -feira 

Nao 

- 

2.177.552 

- 

- 

- 

1.918920 

- 

146,25 

09/08/2005 

3 a -feira 

Nao 

2.036.848 

- 

2.177.552 

- 

- 

- 

1.918.920 

148,50 

10/08/2005 

4 a -feira 

Nao 

801.143 

798.630 

- 

2.177.552 

- 

- 

- 

149,25 

11/08/2005 

5 a -feira 

Nao 

- 

- 

798.630 

- 

2.177.552 

- 

- 

147,00 

12/08/2005 

6 a -feira 

Nao 

- 

1.456.656 

- 

798.630 

- 

2.177.552 

- 

141,75 

13/08/2005 

Sabado 

Nao 

- 

- 

1.456.656 

- 

798.630 

- 

2.177.552 

125,25 

14/08/2005 

Domingo 

Nao 

- 

801.143 

- 

1456.656 

- 

798.630 

- 

110,25 

15/08/2005 

2 a -feira 

Nao 

- 

- 

801.143 

- 

1.456.656 

- 

798.630 

138,75 

16/08/2005 

3 a -feira 

Nao 

296.642 

- 

- 

801.143 

- 

1.456.656 

- 

142,50 

1 7/08/2005 

4 a -feira 

Nao 

246.901 

296.642 

- 

- 

801.143 

- 

1.456.656 

143,25 

18/08/2005 

5 a -feira 

Nao 

- 

- 

296.642 

- 

- 

801.143 

- 

142,50 ) 

( Continua) 
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TABELA 6.1 Planilha de dados do modelo de regressao multipla para prever o TMA (agosto de 2005) (cont.) 


Dia 

Dia da semana 

Feriado? 

Chegada(C) 

V - 3 

V - 2 

V-1 

Vencimento (V) 

V + 1 

V + 2 

TMA (seg) 

19/08/2005 

6 a -feira 

Nao 

- 

- 

- 

296.642 

- 

- 

801.143 

142,50 

20/08/2005 

Sabado 

Nao 

- 

246.901 

- 

- 

296.642 

- 

- 

123,75 

21/08/2005 

Domingo 

Nao 

- 

- 

246.901 

- 

- 

296.642 

- 

109,50 

22/08/2005 

2 a -feira 

Nao 

- 

- 

- 

246.901 

- 

- 

296.642 

139,50 

23/08/2005 

3 a -feira 

Nao 

- 

- 

- 

- 

246.901 

- 

- 

139,50 

24/08/2005 

4 a -feira 

Nao 

2.136.197 

- 

- 

- 

- 

246.901 

- 

139,50 

25/08/2005 

5 a -feira 

Nao 

- 

- 

- 

- 

- 

- 

246.901 

139,50 

26/08/2005 

6 a -feira 

Nao 

- 

- 

- 

- 

- 

- 

- 

133,50 

27/08/2005 

Sabado 

Nao 

- 

- 

- 

- 

- 

- 

- 

118,50 

28/08/2005 

Domingo 

Nao 

- 

- 

- 

- 

- 

- 

- 

100,50 

29/08/2005 

2 a -feira 

Nao 

- 

2.136.197 

- 

- 

- 

- 

- 

139,50 

30/08/2005 

3 a -feira 

Nao 

2.181.734 

- 

2.136.197 

- 

- 

- 

- 

139,50 

31/08/2005 

4 a -feira 

Nao 

1.879401 

- 

- 

2.136.197 

- 

- 

- 

145,50 

_ 
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TABELA 6.2 Coeficientes estimados das variaveis independentes do modelo 

de regressao multipla 

para prever o TMA na data 

em questao 


Coeficientes 

valor-P 

Intersegao 

122,27 

0,00% 

Domingo 

(15,17) 

0,00% 

2 a -feira 

15,37 

0,00% 

3 a -feira 

15,20 

0,00% 

4 a -feira 

17,70 

0,00% 

5 a -feira 

17,86 

0,00% 

6 a -feira 

15,11 

0,00% 

Feriado 

(29,20) 

0,00% 

Chegada (C) 

(0,0000003) 

68,35% 

V - 3 

0,0000008 

29,97% 

V - 2 

0,0000011 

1 6,89% 

V- 1 

0,0000016 

3,44% 

Vencimento (V) 

0,0000029 

0,01% 

V+ 1 

0,0000026 

0,12% 

V + 2 

0,0000022 

0,36% 

V- 


_ 


Os coeficientes de regressao estimados, juntamente com os respectivos valores p, 
estao apresentados na Tabela 6.3. Como pode ser observado, os valores p das outras 
variaveis nao mudaram muito e, para este modelo modificado, apenas uma variavel 
apresenta coeficiente angular com valor p superior (mas nao muito) a 5%: V - 1. No 
entanto, sera mantida no modelo, porque o valor p do seu coeficiente angular esta 
abaixo de 10%. 

Dessa forma, a retirada das variaveis C, V - 3 e V - 2 foi benefica ao modelo, 
tornando-o mais parcimonioso e permitindo uma interpreta£ao dos coeficientes obtidos. 

O primeiro coeficiente corresponde ao valor da variavel dependente TMA quando 
todas as variaveis independentes assumirem o valor zero; ou seja, quando domingo, 
segunda-feira, ter 5 a-feira, quarta-feira, quinta-feira, sexta-feira, feriado, V - 1, V, V 
+ 1 e V + 2 forem zero - ; quando o dia for sabado, nao feriado e nao consistir de 
nenhuma data dentro do intervalo de quatro dias em torno do vencimento de alguma 
conta, o TMA previsto e de 122,64 segundos. 

Esse numero pode ser considerado uma quantidade basica, a qual devem ser 
adicionados os efeitos da ocorrencia das variaveis independentes, quantificados atraves 
dos seus respectivos coeficientes angulares. 
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TABELA 6.3 Coeficientes estimados das variaveis independentes do modelo 
(sem C; V - 3; V - 2) de regressao multipla para prever o TMA na data em questao 


Coeficientes 

valor-P 

Intersegao 

122,64 

0,00% 

Domingo 

(15,15) 

0,00% 

2 a -feira 

15,36 

0,00% 

3 a -feira 

14,96 

0,00% 

4 a -feira 

17,33 

0,00% 

5 a -feira 

17,84 

0,00% 

6 a -feira 

15,03 

0,00% 

Feriado 

(28,89) 

0,00% 

V- 1 

0,0000012 

6,96% 

Vencimento (V) 

0,0000031 

0,00% 

V+ 1 

0,0000032 

0,00% 

V + 2 

0,0000024 

0,04% 

- 


_ 


Assim, o fato de o dia ser domingo contribui com uma diminuiijao de 15,15 segun- 
dos no TMA previsto em rela 9 ao ao numero basico; o fato de o dia ser segunda-feira, 
contribui com um aumento esperado de 15,36 segundos no TMA; se for ter^a-feira, 
com um aumento de 14,96 segundos; se for quarta-feira, com um aumento de 17,33 
segundos; se for quinta-feira, com 17,84 segundos a mais; e se for sexta-feira, contribui 
com um aumento de 15,03 segundos na previsao do tempo medio de atendimento, em 
relagao ao dia basico. 

Com essas informa 9 oes, e possfvel elaborar um grafico (Figura 6.2), que permite 
uma melhor visualiza 9 ao da sazonalidade do dia da semana. A evolu 9 §o da demanda 
de liga 9 oes ao longo da semana, considerando dias sem a ocorrencia de feriados e 
nao proximos do vencimento de contas, come 9 a com um valor baixo no domingo, 
apresentando um certo aumento na segunda-feira; a demanda vai se mantendo cons- 
tante ao longo da semana, tendo uma queda razoavel no sabado, antes de voltar ao 
patamar de domingo, ainda mais baixo. 

De forma analoga aos dias da semana, o fato de um dia ser feriado reduz em 28,89 
segundos o TMA previsto. Esta grande redu 9 ao no tempo medio de atendimento tanto 
para feriados como para finais de semana quando comparados aos dias uteis e, de certa 
forma, estranha. Em princlpio, nao foi possfvel especular os motivos que levariam a 
tal varia 9 ao. A equipe de previsao foi indagada a respeito e tambem nao foi capaz 
de justificar o fato. A investiga 9 ao acerca de tais motivos pode constituir-se em um 
interessante objeto de pesquisa futura. 
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Sazonalidade do dia semana - produto 103 



FIGURA 6.2 Comportamento sazonal do dia da semana - modelo sem C, V - 3 e V - 2. 


Em relacao aos dias que circundam a data do vencimento, os seguintes efeitos 
podem ser explicados, segundo o TMA: cada conta vencendo no dia seguinte (ou 
seja, para as quais a data em questao representa V - 1) aumenta em 0,0000012; cada 
conta vencendo no dia em questao aumenta em 0,0000031; cada conta vencendo na 
vespera aumenta em 0,0000032; e cada conta vencendo na antevespera aumenta em 
0,0000024. 

Naturalmente, os coeficientes nao tern apenas significado isolado, mas tambem 
podem e devem ser analisados em conjunto: se as consequencias dos efeitos dia da 
semana, ocorrencia de feriado e data crftica em relacao ao vencimento da conta forem 
incorporadas ao valor previsto para um dia basico, sera posslvel estabelecer uma 
cquacao para prever o TMA em determinado dia em fungao das caracterfsticas da 
data em questao. Dessa maneira, o TMA pode ser previsto conforme a Equa 9 ao 6.1. 

As sete primeiras variaveis da equagao sao do tipo dummy, assumindo o valor 0 
(nao) ou 1 (sim). As quatro ultimas representam, respectivamente: a quantidade de 
contas vencendo um dia depois; no proprio dia; um dia antes; e dois dias antes. 

EQUA^AO 6.1 Equagao do modelo para prever o TMA 

TMA(emsegundos) = 122,64 -15,15 x domingo + 

15,36 x segunda-feira + 14,96 x terga-feira + 

17,33 x quarta-feira + 17,84 x quinta-feira + 

15,03 x sexta-feira - 28,89 x feriado + 0,0000012 x V - 1 + 

0,0000031 x V + 0,0000032 x V + 1 + 0,0000024 x V + 2 
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De posse dessa equa 9 ao, entao, e posslvel calcular quanto o modelo de regressao 
teria previsto para o TMA em cada um dos dias acerca dos quais o historico dis- 
ponibiliza o valor real desta grandeza. Desta forma, pode-se comparar os valores que 
teriam sido previstos com os que efetivamente ocorreram, de maneira que um erro de 
previsao possa ser mensurado. Esta compara 9 ao pode ser visualizada na Figura 6.3. 

A diferen 9 a entre os valores reais e os gerados pelo modelo tambem e conhecida 
como reslduo. Para ser usado um modelo de Regressao Multipla, e necessario supor 
que a correla 9 ao entre os reslduos, ao longo do espectro das variaveis independentes, 
e zero; ou seja, que os reslduos sao independentes entre si, nao apresentando, conse- 
quentemente, autocorrela 9 ao serial. 3 

Uma forma de verificar a validade dessa suposi 9 ao e atraves do teste de Dur¬ 
bin-Watson. No modelo que esta sendo analisado, foi encontrado o valor 1,756 para 
a Estatlstica DW. Esta estatlstica obtida deve ser comparada com valores crlticos que 
variam em fmnjao do nlvel de significancia do teste, do tamanho da amostra e do 
numero de variaveis independentes. 4 

Com 5% de significancia, 234 dados na amostra ell variaveis independentes, o 
valor crltico inferior (diJ obtido foi 1.688 e o valor crltico superior (d (J ) encontrado 
foi 1,885. 5 


TMA (real x previsto) por dia (produto 103) 

180-,- 



-Real - - - Previsto 


30 


0-1-1-,-,-i-1-,-i-■- 

17/07/2005 16/08/2005 15/09/2005 15/10/2005 14/11/2005 14/12/2005 13/01/2006 12/02/2006 14/03/2006 

Dla 

FIGURA 6.3 Valores reais multiplicados por uma constante e previstos pelo modelo para o TMA 
por dia (julho de 2005 a mar^o de 2006). 


3. Corrar; Paulo; Dias Filho (2007). 

4. Corrar; Paulo; Dias Filho (2007). 

5. Valores encontrados em tabela dispomvel no endere 90 <http://www.stanford.edu/~clint/ bench/dwcrit. 
htm>. Acesso em 18/08/2009. 
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Plotagem de residuos 
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Domingo 

Residuos do modelo X variavel dummy Domingo. 


Como o valor obtido para a Estatistica DW (1,756) se situa entre os dois valores 
crfticos encontrados, o teste para verificar a existencia de autocorrelagao dos residuos e 
nao conclusivo: nao e possivel garantir que os residuos sao independentes, mas tambem 
e incorreto afirmar que existe uma relagao de dependencia entre eles. 6 

Sendo assim, parece valido examinar o grafico do residuo contra cada variavel 
independente a fim de procurar uma suposta autocorrelagao e a identificagao de uma 
eventual lei de formagao nos valores dos residuos a medida que o valor da variavel 
independente avanga. A Figura 6.4 ilustra o caso da variavel dummy domingo. Como 
pode ser observado, nao parece existir uma lei de formagao; ou seja, a magnitude dos 
residuos nao parece depender do valor da variavel independente (domingo). 

Os residuos apresentaram comportamento bem parecido em relagao as outras 
variaveis dummy (feriado e as referentes aos outros dias da semana). E visto que o 
modelo trabalha com 11 variaveis independentes, seria, de certa maneira, inviavel 
apresentar aqui tal quantidade de graficos. Desta forma, sera apresentada graficamente 
apenas a relagao entre os residuos e a outra variavel independente mais marcante: 
Vencimento (Figura 6.5). 

O comportamento dos residuos em relagao as outras variaveis (V-l,V+leV = 2) 
e bastante similar ao verificado na Figura 6.5: parecem estar distribuidos aleatoriamente 
em torno do eixo das abscissas, nao apresentando nenhuma tendencia. 

Tudo isso gera uma seguranga maior para os resultados do modelo, ja que o teste 
para a verificagao de autocorrelagao serial dos residuos foi nao conclusivo. E claro 
que esta ultima analise e meramente visual, carecendo de um cientificismo maior. Mas 
e importante proceder dessa forma principalmente em situagoes em que os criterios 


6. Corrar; Paulo; Dias Filho (2007). 
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Plotagem de residuos 



objetivos apontam para uma indefinicao, como foi o caso. Em relafao a qualidade 
preditiva do modelo, foi possfvel perceber, na Figura 6.3, excluindo alguns pontos 
especfficos, uma boa aderencia visual dos valores previstos em rela£ao aos dados reais. 
Esta aderencia foi aqui formalizada atraves do MAPE (Mean Absolute Percentual 
Error), ou erro medio percentual absoluto. 

A tftulo de ilustra£ao, a Tabela 6.4 apresenta dados relativos ao mes de janeiro de 
2006, em que constam o TMA real, os valores que o modelo de regressao teria previsto 
para esta grandeza, assim como os erros percentuais absolutos oriundos desta previsao. 
Pode-se calcular que o modelo teria errado em 4,2%, em media, na previsao do TMA 
durante o mes em questao. 

Mas janeiro e um mes atfpico, com comportamento geralmente mais diffcil de 
prever, sobretudo nos primeiros dias. Tal afirmacao pode ser verificada neste exemplo 
particular a partir da Figura 6.6, na qual e possfvel perceber que todos os outros meses 
relacionados, com exce 9 ao de dezembro, outro perfodo atfpico, apresentam erros 
medios de previsao menores do que o erro medio de janeiro de 2006. 

Em media, o modelo esta incorrendo em erro de previsao da ordem de 3,61%, valor 
satisfatoriamente baixo. Antes de compara-lo com o erro obtido pelo processo-padrao 
de previsao, e valido lembrar, por exemplo, que, durante este processo, os analistas 
recebem os valores previstos pelo software e acrescentam a eles os impactos esperados 
por eventos e premissas, como o comportamento do ciclo de contas e a ocorrencia de 
feriados, assim como outros eventos especiais relacionados ao produto, como o piano 
de minutos, campanhas na mfdia e outros potenciais aspectos impactantes. 

Ora, o modelo de regressao proposto visa substituir apenas os valores sugeridos 
pelo software, acrescidos da influencia dos feriados e do comportamento dos ciclos de 
contas; nao se teve por pretensao substituir o papel do analista encarregado da previsao, 
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TABELA 6.4 Valores reais, previstos e erros do modelo de regressao multipla 
para prever o TMA (janeiro de 2006) 



Dia da 



TMA (seg) Previsto 

Erro de 

Dia 

semana 

Feriado? 

TMA (seg) 

pelo modelo) 

previsao 

01/01/2006 

Domingo 

Sim 

93,00 

78,60 

15,5% 

02/01/2006 

2 a -feira 

Nao 

138,75 

137,99 

0,5% 

03/01/2006 

3 a -feira 

Nao 

141,75 

137,60 

2,9% 

04/01/2006 

4 a -feira 

Nao 

146,25 

139,97 

4,3% 

05/01/2006 

5 a -feira 

Nao 

145,50 

140,48 

3,5% 

06/01/2006 

6 a -feira 

Nao 

145,50 

137,67 

5,4% 

07/01/2006 

Sabado 

Nao 

134,25 

122,64 

8,6% 

08/01/2006 

Domingo 

Nao 

116,25 

107,49 

7,5% 

09/01/2006 

2 a -feira 

Nao 

150,00 

137,99 

8,0% 

10/01/2006 

3 a -feira 

Nao 

149,25 

137,60 

7,8% 

11/01/2006 

4 a -feira 

Nao 

149,25 

139,97 

6,2% 

12/01/2006 

5 a -feira 

Nao 

144,00 

140,48 

2,4% 

13/01/2006 

6 a -feira 

Nao 

141,75 

137,67 

2,9% 

14/01/2006 

Sabado 

Nao 

126,75 

122,64 

3,2% 

15/01/2006 

Domingo 

Nao 

108,00 

107,49 

0,5% 

16/01/2006 

2 a -feira 

Nao 

140,53 

137,99 

1,8% 

17/01/2006 

3 a -feira 

Nao 

138,33 

137,60 

0,5% 

18/01/2006 

4 a -feira 

Nao 

137,94 

139,97 

1,5% 

19/01/2006 

5 a -feira 

Nao 

137,13 

140,48 

2,4% 

20/01/2006 

6 a -feira 

Nao 

129,92 

137,67 

6,0% 

21/01/2006 

Sabado 

Nao 

119,52 

122,64 

2,6% 

22/01/2006 

Domingo 

Nao 

103,92 

108,36 

4,3% 

23/01/2006 

2 a -feira 

Nao 

135,62 

140,21 

3,4% 

24/01/2006 

3 a -feira 

Nao 

134,42 

139,84 

4,0% 

25/01/2006 

4 a -feira 

Nao 

134,69 

141,70 

5,2% 

26/01/2006 

5 a -feira 

Nao 

137,07 

140,48 

2,5% 

27/01/2006 

6 a -feira 

Nao 

135,58 

137,67 

1,5% 

28/01/2006 

Sabado 

Nao 

122,99 

122,64 

0,3% 

29/01/2006 

Domingo 

Nao 

107,68 

107,49 

0,2% 

30/01/2006 

2 a -feira 

Nao 

133,53 

137,99 

3,3% 

31/01/2006 

3 a -feira 

Nao 

127,10 

140,22 

10,3% 
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Erro medio de previsao do modelo (mes a mes) 



Mds/ano 

FIGURA 6.6 Erro medio mensal do modelo de regressao multipla para prever o TMA (julho de 
2005 a fevereiro de 2006). 

dotado de toda a sua valiosa experiencia e subjetividade, mas, sim, fornecer-lhe um 
subsfdio complementar para a tomada de decisao. Assim, e ja que o analista estara 
apto a acrescentar aos valores sugeridos pelo modelo de regressao os impactos que 
julgar relevantes, e mais justo comparar os erros de previsao do modelo aqui proposto 
com os erros oriundos do modulo de previsao do software, com o intuito de verificar 
a ocorrencia ou nao de melhorias no processo. 

Assim sendo, o MAPE do software para o mesmo perfodo analisado foi de 5,58%, 
razoavelmente superior, portanto, aos 3,61% do modelo de regressao. 

Alem disso, uma analise mais detalhada dos erros de previsao do modelo proposto 
revela que alguns dos dias mais dificeis de prever e, portanto, com erros de previsao 
potencialmente maiores, consistem em feriados ou datas que se comportam quase como 
um feriado. Por exemplo, os dias 24/12/2005 e 31/12/2005 nao foram classificados 
como feriados, mas, certamente tambem nao sao dias comuns e devem gerar um TMA 
diferente daquele de um dia-padrao. O problema e que o modelo de regressao nao levou 
em conta isso, ja que era possfvel classificar um dia apenas como sendo feriado ou 
nao; e ja que esses dias nao eram efetivamente feriados, foram tratados pelo modelo 
como dias comuns. 

Problemas semelhantes, mas talvez em menor escala, certamente ocorreram em 
alguns dias do ano que consistiram em feriados locais, mas nao nacionais. Como 
o modelo propoe uma previsao consolidada do TMA, so foi possfvel considerar 
os feriados nacionais; por exemplo, o dia 20/01/2006, um feriado local na Cidade 
do Rio de Janeiro, teve uma alterafao no TMA nacional provocada pela alterafao 
no TMA na cidade mas, ao ser tratado como um dia normal pelo modelo, teve seu 
TMA mal estimado. Este problema tambem se observa nos feriados prolongados 
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e nos dias que se situam entre um feriado e o fim de semana, em que o TMA foge 
do padrao. 

Um problema ainda mais grave ocorreu com o dia 25/12/2005. Seu erro de previsao 
foi muito alto porque o dia em questao apresentava dois efeitos redutores de TMA, que 
foram somados pelo modelo de natureza aditiva: ser domingo e ser feriado. Mas, na 
pratica, o efeito de um dia ser feriado certamente e reduzido se ele cair em um domingo 
ou sabado; e o modelo nao levou isso em considera^o, ate porque so e possivel consi- 
derar um dia como sendo feriado ou nao, e nao como uma especie de feriado moderado. 

Na verdade, esses problemas aumentaram os erros de previsao, mas isso nao 
quer dizer que o processo futuro de previsao tera que incorrer nesses erros de grande 
magnitude quando ocorrerem dias com comportamentos diferentes. Na pratica, quando 
os analistas de previsao estiverem utilizando o modelo aqui proposto, eles poderao 
manipular os valores da variavel dummy feriado para os dias diferentes do padrao. 

Assim, por exemplo, o dia 24/12 pode receber manualmente o valor 0,7 para a sua 
variavel feriado, em vez de ser tratado como um dia normal; o dia 20/01 pode receber 
o valor 0,4 da mesma forma; um feriado que caia em um domingo pode ter o seu valor 
da variavel feriado reduzido de 1 para 0,6. 

Depois dessa manipulacao, basta considerar as informa 9 oes do setor de faturamento 
da Telemar sobre os ciclos de contas e usar a Equaijao 1 apresentada anteriormente 
para prever o TMA para o produto 103, a cada dia. Uma vez gerados os outputs da 
equa 5 ao, os analistas podem incorporar os impactos esperados pelos eventos especiais, 
como implanta 9 ao de novos servi 90 s, mudan 9 as no roteiro de atendimento etc. e, desta 
forma, reduzir substancialmente os erros de previsao. 

EXERCICIOS 

1 . A planilha Exercicio 1 relaciona diversas caracterfsticas de 151 automoveis, assim 
como seus pre 90 s de mercado (em US$). 

a. Construa um modelo de Regressao Multipla para explicar o comportamento da 
variavel pre 90 em fun 9 ao das variaveis relacionadas ao desempenho do carro. 

b. O modelo ficou bom? Todas as variaveis sao bons “explicadores lineares” do 
comportamento do pre 90 , a 10% de significance? 

c. Retire alguma(s) variavel(is) que nao esteja(m) explicando muito bem o pre 90 
e refa 9 a a regressao. Avalie o novo modelo. 

2. A companhia que vende leite deseja conhecer uma possivel rela 9 ao entre as quanti- 
dades vendidas e o valor investido em propaganda no mes. A planilha Exercicio 2 
apresenta os valores relativos aos ultimos 14 meses. Sabe-se que a empresa en- 
frentou uma greve do quinto ao setimo mes. 

a. Ignore a informa 9 ao da greve e rode uma Regressao Simples para explicar a 
quantidade vendida em fun 9 ao do investimento em propaganda. 

b. Refa 9 a a Regressao Simples, eliminando os meses de greve. Avalie o resultado. 

c. Trate a greve como uma variavel dummy e fa 9 a uma Regressao Multipla. 
Quanto seria vendido em um mes normal, em que fossem investidos R$ 1.000 
em propaganda? 
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3. A planilha Exercfcio 3 mostra a demanda de energia e a temperatura (em °F) para 
diversos dias e perfodos do dia em determinada regiao. 

a. Construa um modelo de Regressao Multipla para prever a demanda horaria 
de energia eletrica em fungao do perfodo do dia (madrugada, manha, tarde ou 
noite), do dia da semana e da Temperatura (°F). Avalie o modelo. 

b. Quanto seria a demanda das 14 as 15 horas de uma quarta-feira, com tempera¬ 
tura media de 80°F? 

c. E de 20 as 21 horas de um sabado, com temperatura media de 15°F? 
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Objetivos de aprendizagem do apendice 

• Aprender a prever os valores futuros de uma grandeza com base apenas no seu 
complexo comportamento historico. 

• Conhecer os componentes mais comuns das series temporais. 

• Saber isolar cada um desses componentes e entender a influencia deles no compor¬ 
tamento da grandeza estudada. 

• Aprender a estimar o comportamento futuro de cada um desses componentes. 

• Conseguir agrupar todos esses componentes novamente para prever a grandeza 
estudada. 

• Entender a influencia do acaso no comportamento da serie temporal e utilizar esse 
conhecimento para construir cenarios pessimistas e otimistas de previsao. 




Muitas grandezas apresentam um comportamento bastante irregular ao longo do tempo. 
Por exemplo, a venda de determinado produto ao longo de um mes em determinada 
filial de uma loja; ou a temperatura ao longo do dia em um determinado bairro; ou a 
intensidade pluviometrica ao longo do ano em uma determinada cidade; entre outras. 

Ao observarmos esse comportamento irregular (atraves de graficos, por exemplo), 
podemos imaginar a dificuldade em encontrar um padrao de comportamento para a 
grandeza em voga. Sem identificar esse padrao de comportamento, fica muito dificil, 
por exemplo, fazer um exercfcio de previsao acerca do futuro da grandeza. 

Imaginemos que o Grafico 1 representa o consumo de energia eletrica no Brasil ao 
longo das semanas. Esse comportamento irregular pode parecer desprovido de padrao 
mas, como veremos mais adiante, nao e exatamente o que ocorre. 

O metodo apresentado a partir de agora e conhecido como Decomposi^ao Classica 
e tem como objetivo entender o comportamento (decompondo-o) das grandezas ao 
longo do tempo, de modo a auxiliar na previsao para os valores futuros das grandezas. 

A irregularidade no comportamento das grandezas normalmente e resultado de 
diversos componentes. No exemplo do consumo de energia eletrica (Grafico 1), esta 
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grandeza varia ao longo do tempo em fungao da sazonalidade, da tendencia, do ciclo e 
de componentes aleatorios. Estes componentes aparecem em parte dos comportamentos 
irregulares apresentados pelas grandezas. Nesses casos, a Decomposigao Classica e o 
metodo de previsao mais adequado, ja que e nesses elementos que as grandezas sao 
decompostas neste metodo. A seguir, esses componentes sao explicados. 

SAZONALIDADE (S) 

Explica o comportamento cfclico (ao longo de um ano ou intervalo menor) das gran¬ 
dezas, muitas vezes sendo resultado de suas variagoes em fungao das estagoes do ano 
(por exemplo: venda de sorvetes, consumo de energia eletrica). Nestes casos, o padrao 
de sazonalidade e mensal ou trimestral ao longo do ano. A venda de sorvetes no Brasil, 
por exemplo, e bem grande no primeiro trimestre de cada ano, diminui no segundo, 
atinge seu menor valor no terceiro e volta a aumentar no quarto. 

Mas a sazonalidade pode se manifestar de forma diferente da mensal/semanal e/ 
ou nao ser fruto das variagoes climaticas. Exemplos: 

• as vendas semanais de uma empresa ao longo do mes - que podem ser maiores nas 
ultimas semanas e menores nas primeiras em razao do cumprimento de cotas por 
partes dos vendedores; 

• consumo diario de bebida alcoolica nos bares ao longo da semana - que e maior nos 
ultimos dias da semana e menor nos primeiros por causa do movimento nos bares; 

• a intensidade de trafego por hora ao longo do dia - que apresenta picos nos horarios 
de “rush” - de manha cedo e no fim da tarde. 

TENDENCIA (T) 

As grandezas, em boa parte das vezes, apresentam uma tendencia de comportamento 
em relagao a variavel tempo. Assim, nesses casos, elas revelam um comportamento que 
tende a ser maior (na maioria das vezes) ou menor a medida que o tempo passa. Esta 
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tendencia, normalmente, e consequencia do aumento da populacao e/ou da evolu 9 ao 
tecnologica/cientffica, mas tambem pode ser fruto de outras causas, como o aumento 
gradativo da temperatura na Terra ou ate mesmo de mudanpas no comportamento 
cultural das sociedades. 

Dessa forma: 

• a receita de vendas de uma empresa tende a crescer por causa do aumento da 
populacao; 

• o percentual de obitos em uma cidade tende a diminuir por causa da evolufao 
cientffica; 

• o percentual de mulheres trabalhando fora de casa tende a aumentar por causa de 
modifica 9 oes culturais. 

E, por causa disso, e de se esperar, a priori, que a receita de vendas de uma empresa 
tenha sido maior em 2001 do que em 1982; e que o percentual de obitos em uma cidade 
tenha sido menor em 2010 do que em 1950. 

CICLO (C) 

Parte do comportamento irregular das grandezas pode ser atribufda ao ciclo, que 
se manifesta na forma de flutua 9 oes semelhantes a ondas em torno das tendencias, 
que tern dura 9 ao maior do que um ano, mas que nao tern a exata periodicidade que 
caracteriza o padrao sazonal. A identifica 9 ao das epocas em que ocorrem os picos e 
depressoes destas ondas sugere seu relacionamento com a flutua 9 ao geral da atividade 
economica - o chamado ciclo economico. 

Portanto, mesmo que uma grandeza nao revele sazonalidade nem tendencia quando ob- 
servada a sua varia 9 ao em rela 9 ao ao tempo, ainda assim pode apresentar um comportamento 
cfclico nao aleatorio, que pode ser visualizado, entendido e decomposto de modo a ser usado 
num futuro exercfcio de previsao. O montante per capita de dinheiro gasto em restaurantes 
em Sao Paulo, por exemplo, nao apresenta, a principio, nenhum padrao sazonal e nenhuma 
tendencia de crescimento ou redu 9 ao ao longo do tempo; mas seu comportamento pode ser 
explicado pela varia 9 ao do m'vel economico e ter sua varia 9 ao atribufda exclusivamente 
ao ciclo (alem de a componentes aleatorios). Se for esperada uma atividade economica de 
intensidade em 2006 maior do que em 1999 e menor do que em 2003, ha de se prever que 
a grandeza em questao seja maior em 2003 e menor em 1999 (com um valor intermediario 
para 2006), independente do mes do ano. 

COMPONENTES ALEATORIOS (A) 

Se expurgarmos os efeitos sazonais, cfclicos e de tendencia de uma serie que revele 
o comportamento de uma grandeza ao longo do tempo, o restante da sua varia 9 ao e 
denominado flutua 9 ao residual e pode ser considerado como proveniente de compo¬ 
nentes aleatorios. Esta flutua 9 ao residual tambem pode ser entendida como oriunda 
de situa 9 oes indistingufveis, ja que a sazonalidade, a tendencia e o ciclo foram consi- 
derados como situa 9 oes distingufveis. 

Como dito anteriormente (e sera detalhadamente explicado a partir de agora, com 
o auxflio do exemplo da venda de lubrificantes no Brasil e da planilha Deeomposiijao 
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classica.xls), o metodo de Decomposigao Classica isola os efeitos (sazonalidade, 
tendencia, ciclo e componentes aleatorios) que geram a irregularidade do comporta- 
mento das series temporais, preve seu comportamento (de cada efeito, isoladamente) 
no futuro e, novamente, agrupa esses efeitos de forma a obter a melhor previsao para 
o futuro da serie temporal. 

A serie temporal pode, entao, ser representada da seguinte forma: 

R=SxTxCxA 

em que Reo valor real da grandeza e os fatores sao os efeitos mencionados ante- 
riormente 

O METODO DA DECOMPOSigAO CLASSICA 

Ao abrirmos a planilha Decomposigao classica.xls, podemos observar a serie de dados 
historicos que representa a venda de lubrificantes no Brasil, por trimestre, entre 1983 e 
1999. O primeiro passo para a elaboragao de uma previsao baseada em dados passados 
e fazer um grafico da serie historica, de modo a melhor perceber o seu comporta¬ 
mento. Fazemos isso pressionando o botao Executar passo da planilha. 

Ao observarmos o grafico, podemos verificar uma flutuagao muito grande dos 
dados, alem de uma tendencia relativamente nftida de crescimento. Para decompor- 
mos a serie nos seus quatro efeitos ja mencionados - tendencia, ciclo, componentes 
aleatorios e sazonalidade -, vamos comegar por este ultimo. 

Os dados sao apresentados em trimestres e parece que e nessa unidade de tempo 
que a sazonalidade se manifesta. Se compararmos as vendas totais do ano de 1984 
com as vendas totais dos anos de 1985, podemos entender que essa comparagao esta 
livre de sazonalidade, porque ambos os anos tem um primeiro trimestre, um segundo 
trimestre, um terceiro trimestre e um quarto trimestre, eliminando as diferengas devidas 
a esse efeito. Na verdade, considerando o intervalo de um ano, boa parte do efeito 
aleatorio e anulada, porque os componentes positivos e negativos tendem a se anular 
ao longo de um ano. 

Entretanto, e interessante continuarmos trabalhando na mesma grandeza dos dados 
originais - a venda trimestral. Para tal, nao lidaremos com a venda anual, mas, sim, 
com a media movel de quatro perfodos (trimestres). Assim, a primeira media movel 
sera a media de vendas do primeiro ao quarto trimestres; a segunda media movel sera 
a media de vendas do segundo ao quinto trimestres; e assim por diante. Podemos 
observar que cada uma dessas medias moveis contem um primeiro trimestre, um 
segundo trimestre, um terceiro trimestre e um quarto trimestre, anulando o efeito da 
sazonalidade. 

Nosso proximo passo e, entao, calcular essas medias moveis, pressionando no¬ 
vamente o botao Executar passo da planilha (esse procedimento devera ser repetido 
sempre que estiver sendo indicado o proximo passo). A formula das medias moveis 
pode ser observada na planilha. A rigor, a primeira media movel se situa entre o 
segundo e o terceiro trimestres (ponto medio do intervalo), a segunda entre o terceiro 
e o quarto trimestres, e assim por diante. Para efeito de apresentagao, elas foram 
colocadas “meia linha” acima na planilha, ou seja, a primeira media movel na segunda 
linha, a segunda media movel na terceira linha, e assim por diante. 
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Para fazer esse ajuste, e necessario avai^ar para o passo seguinte, que consiste em 
calcular as medias moveis centradas. A media movel centrada corresponde a media 
entre duas medias moveis consecutivas e se situa no ponto medio entre elas (sua formula 
pode ser observada). Assim, a primeira media movel centrada se situa na terceira linha 
(ja que consiste na media das duas primeiras medias moveis - situadas entre a segunda 
e a terceira linhas, e entre a terceira e a quarta linhas, respectivamente); a segunda 
media movel centrada na quarta linha; e assim por diante. 

O proximo passo e plotar as medias moveis centradas no grafico, de modo a 
observar como seu comportamento esta livre da sazonalidade (e de boa parte dos 
componentes aleatorios tambem). Observamos que a flutuafao da serie diminuiu bas- 
tante e podemos concluir que a sua varia 9 ao e fruto da tendencia e do ciclo economico. 
Esta conclusao pode ser mais bem entendida - ao consideramos que a media movel 
consiste nos dados originais (R), sem o efeito da sazonalidade e dos componentes 
aleatorios - atraves da formula: 

, SxTxCxA _ _ 

Media Movel =-= T x C 

Sx A 

Como se trata de um exercfcio de decompos^ao, esta na hora de decompor a media 
movel (T X C) em tendencia (T) e ciclo (C). Para tal, devemos executar o proximo 
passo, que e determinar a tendencia. Isso e feito a partir da op 9 ao (do Excel) de 
adicionar uma linha de tendencia a serie da media movel centrada. A reta plotada e 
crescente (revelando uma tendencia positiva, provocada, provavelmente, pelo aumento 
da popula 5 ao brasileira e do seu poder aquisitivo no perfodo da serie) e pode ser re- 
presentada pela equa^ao y = 61,324x - 118.654. Esta reta explica, aproximadamente, 
72% dos dados (R 2 = 0,7238: que pode ser considerado um bom nfvel de explica^ao) 
que compoem a serie das medias moveis centradas. 

O proximo passo e quantificar os valores dessa reta de tendencia, de modo que 
possamos ter dados para cada um dos trimestres do perfodo considerado. Essa quanti- 
fica§ao e realizada por meio da equa^ao da reta, anteriormente apresentada, que pode 
ser observada nas formulas da planilha. 

Essa quantifica 9 ao e importante como prepara 9 ao para o proximo passo, a deter- 
mina 9 ao dos residuos ciclicos. Como sabemos, a media movel consiste nos efeitos 
da tendencia e do ciclo. Portanto, se eliminarmos o efeito da tendencia (T), teremos 
apenas os residuos ciclicos (C): 

Media Movel _ Tx C 

T T 

Essa rela 9 ao tambem pode ser observada nas formulas da planilha. Esses residuos 
ciclicos consistem, na verdade, de flutua 9 oes em torno da tendencia, resultantes da 
maior ou menor intensidade do ciclo economico. Assim, o primeiro trimestre de 1986 
viveu uma situa 9 ao de baixa intensidade do ciclo economico (88,7%), reduzindo em 
11%, aproximadamente, o nfvel de vendas projetado pela tendencia para esse perfodo. 
Analogamente, o segundo trimestre de 1989 experimentou uma forte intensidade 
economica (109,3%), aumentando em 9%, aproximadamente, o nfvel de vendas 
projetado pela tendencia para este perfodo. 
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Essas flutuagoes podem ser mais bem evidenciadas depois de executado o passo 
seguinte: elaboragao do grafico dos residuos ciclicos. Nesse grafico, a tendencia esta 
representada pela linha horizontal de 100% e os residuos podem ser observados como 
flutuagoes em torno dessa linha. Assim, pode-se ver que o pico economico ocorreu 
no terceiro trimestre de 1965 (118,7%) e o menor nlvel do ciclo aconteceu no terceiro 
trimestre de 1983 (77,8%). 

Com a media movel, reunimos os efeitos da tendencia (T) e do ciclo economico 
(C), conforme ja constatado. Se dividirmos os valores originais (R) pelos valores da 
media movel, ficaremos apenas com o efeito da sazonalidade e dos componentes 
aleatorios, ja que: 


R SxTxCxA 0 , 

_—_= S x A 

Media Movel T x C 

E e esse o proximo passo: a avaliagao do efeito da sazonalidade e dos componentes 

aleatorios. As formulas podem ser observadas na planilha. 

Assim como fizemos anteriormente, ao separarmos os efeitos da tendencia dos efei¬ 
tos do ciclo, nos o faremos com o efeito da sazonalidade e dos componentes aleatorios. 
Comegaremos calculando os Indices sazonais (sao quatro: um para cada trimestre), 
que representam variagoes relativas entre os trimestres. A media de todos os Indices 
sazonais deve ter sempre o valor 1, ja que, ao considerarmos o perfodo de sazonalidade 
(neste caso, o ano inteiro), o efeito sazonal deve sumir (conforme ja mencionado). 

Ao agruparmos todos os valores (de S X A) por trimestre e tirarmos sua media, 
obteremos os indices sazonais, pois os componentes aleatorios tendem a se anular ao 
juntarmos 15, 20 valores, e o resultado e proveniente apenas de efeitos sazonais. Na 
verdade, e conveniente, no calculo dessas medias, eliminar os valores extremos (o 
menor e o maior) dos grupos de modo a se evitar considerar outliers (originados por 
fatores como greves e incendios), ja que um desvio muito grande pode distorcer uma 
media adequada. 

Assim, o indice sazonal do primeiro trimestre consiste na media dos valores (de 
S x A) dos primeiros trimestres (com excegao do menor e do maior). A formula e 
analoga para os outros trimestres e pode ser observada na planilha apos a execugao do 
passo seguinte (calculo dos indices sazonais). Esses indices se repetem, ja que existe 
apenas um indice sazonal para o primeiro trimestre, um indice sazonal para o segundo 
trimestre, e assim por diante. 

Como a media utilizada e modificada (em fungao da exclusao dos valores ex¬ 
tremos), e necessario executar o proximo passo e ajustar os indices sazonais. Este 
ajuste consiste em dividir cada indice sazonal pela media dos quatro indices (conforme 
formulas na planilha) de modo a fazer com que sua media seja igual ao valor 1. Pelo 
fato de as medias nao terem contido os valores extremos, a media dos indices sazonais 
ficou, nesse caso, um pouco maior do que 1 (1,001), devendo os indices ser ajustados 
(reduzidos) por esse fator. 

Nesse momento, podemos observar que os maiores indices sazonais ocorrem no 
segundo e no quarto trimestres, respectivamente. Esse comportamento pode ser enten- 
dido atraves do fato de as distribuidoras de petroleo trabalharem com cotas semestrais 
de vendas, deslocando os picos de venda para os finais dos semestres. 
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Analogamente a separagao dos efeitos do ciclo dos efeitos da tendencia, a separagao 
da sazonalidade dos componentes aleatoria e mais bem visualizada graficamente. O 
passo seguinte consiste em elaborar os graficos dos componentes aleatorios, um para 
cada trimestre. Assim, no grafico do primeiro trimestre, a reta horizontal representa 
o indice sazonal do primeiro trimestre e a flutuagao dos valores (de S x A) em torno 
dela e explicada exclusivamente por componentes aleatorios. Para os outros trimes- 
tres, a ideia e a mesma. 

Naturalmente, o proximo passo e calcular os componentes aleatorios, na forma de 
uma divisao dos valores (de S x A) pelos indices sazonais (S), conforme as formulas 
na planilha e a seguir: 


R 1 SxA A 

_X_—_— 

Media Movel S S 

Esses componentes aleatorios devem ser mais bem observados em uma distribuigao 
acumulada de probabilidade (cuja elaboragao consiste no passo seguinte), j a que serao 
utilizados na confecgao de uma previsao probabilistica de vendas, conforme veremos 
mais adiante. Assim, ao analisarmos alguns pontos notaveis desse grafico, podemos 
concluir que: existem 47% de probabilidade de o componente aleatorio ser menor do 
que 100%, ou seja, de o valor real futuro de venda ser inferior ao da previsao pontual 1 
(desconsiderando os componentes aleatorios); e que existem 5% de probabilidade de 
o componente aleatorio ser maior do que 109%. 

Para fazermos uma previsao de vendas para o futuro, imaginaremos que o padrao 
sazonal se mantera no periodo previsto. Assim, o proximo passo e repetir os indices 
sazonais para os trimestres futuros. 

E de se esperar, tambem, que a tendencia se manifeste - no periodo previsto - da 
mesma forma que se manifestou no passado. Desta forma, o passo seguinte consiste 
na extrapolagao da tendencia futura. Podemos observar que os valores da tendencia 
sao resultado da mesma formula que compos os valores passados. 

O proximo passo (preenchimento dos valores futuros do ciclo economico) requer 
um pouco de informagoes, bom senso e imaginagao por parte do previsor, ja que nao 
ha motivos para crer que essa variavel apresente, no futuro, o mesmo padrao de corn- 
portamento do passado. Desta forma, o previsor pode consultar especialistas da area 
economica, ou adotar um valor unico que represente a intensidade com a qual o ciclo 
economico se manifestara no periodo previsto (este valor pode ser a media dos valores his- 
toricos ou um outro valor que o previsor julgue ser mais representativo do comportamento 
futuro). No nosso exemplo, o previsor reuniu informagoes que o levaram a considerar que: 

• o ano de 2000 apresentaria um nivel economico um pouco mais intenso (88%) que 
o de 1999; 

• esse nivel cairia para 83% em 2001; 

• esse nivel subiria para 95% em 2002. 


1 A previsao de vendas pontual ignora os componentes aleatorios, e realizada compondo-se os efeitos de 
sazonalidade, tendencia e ciclo e sera elaborada a partir de agora, para o periodo de 2000 a 2002 (tres anos 
ou 12 trimestres). 
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A composigao das vendas futuras e executada no passo seguinte, simplesmente 
agrupando (compondo) os tres fatores: T x C x S, formula que pode ser observada na 
planilha. Essa e a previsao pontual, que ignora o efeito dos componentes aleatorios. 

A previsao probabilfstica de vendas considera esse efeito (e reproduz da mesma 
forma que na pontual o comportamento dos outros tres efeitos) e sera ilustrada por 
tres cenarios: 

• Pessimista, confeccionado no passo seguinte, em cor vermelha, para o qual o 
componente aleatorio escolhido foi o segundo menor (para retirar o valor extremo) 
entre os dados historicos (84,9%). Ao observarmos esse valor na distribuigao 
acumulada de probabilidade, podemos perceber que vendas reais inferiores a esses 
valores previstos tem menos de 3% de probabilidade de ocorrer. 

• Otimista, confeccionado no proximo passo, em cor verde, para o qual o componente 
aleatorio escolhido foi o segundo maior (para retirar o valor extremo) entre os 
dados historicos (110,3%). Ao observarmos esse valor na distribuigao acumulada 
de probabilidade, podemos perceber que vendas reais superiores a esses valores 
previstos tem, aproximadamente, 5% de probabilidade de ocorrer. 

• Meio-termo, confeccionado no passo seguinte, em cor amarela, para o qual o 
componente aleatorio escolhido foi o valor 104,2%. Ao observarmos esse valor 
na distribuigao acumulada de probabilidade, podemos perceber que vendas reais 
inferiores a esses valores previstos tem, aproximadamente, 72% de probabilidade 
de ocorrer. 

O ultimo passo consiste em plotar esses cenarios (juntamente com a previsao 
pontual) no grafico de vendas, de modo a comparar os posslveis comportamentos 
futuros dessa grandeza com o seu comportamento historico. O previsor pode e deve, 
inclusive, observar entre que faixas de valores as vendas se situarao no futuro, e com 
que probabilidade de ocorrencia. 

EXERCICIOS 

1 . A planilha Exereicio 1 relaciona as vendas mensais de todo o varejo dos Estados 
Unidos (em bilhoes de dolares) de janeiro de 1984 a julho de 1993. Faga a previsao 
do comportamento dessa grandeza para os dois anos seguintes. 

2. A planilha Exereicio 2 relaciona as vendas mensais de uma cadeia de lojas de 
departamento de janeiro de 1993 a dezembro de 1995 no Brasil. Faga a previsao 
do comportamento dessa grandeza para os dois anos seguintes. 

3. A planilha Exereicio 3 relaciona as vendas mensais (em unidades) de uma cadeia 
de lojas de artigos femininos de janeiro de 1986 a dezembro de 1992 no Brasil. 
Faga a previsao do comportamento dessa grandeza para os dois anos seguintes. 

4. A planilha Exereicio 4 relaciona as vendas mensais (em 1 .000m 3 ) de oleo diesel no 
Brasil de janeiro de 1984 a setembro de 1994. Faga a previsao do comportamento 
dessa grandeza para os dois anos seguintes. 
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5. A planilha Exercicio 5 relaciona a producao mensal (em 1.000 litros) de cerveja na 
Australia de janeiro de 1956 a agosto de 1995. Fa§a a previsao do comportamento 
dessa grandeza para os dois anos seguintes. 
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